Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,089,761 G→A intergenic (‑30208/‑42342) unknown ← / → unknown tRNA‑Ser‑GGA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,089,7610GA100.0% 194.2 / NA 72intergenic (‑30208/‑42342)unknown/unknowntRNA‑Ser‑GGA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/40);  total (32/40)

CAGGTGGTTAATGAAGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATGCTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCGAAGAAGAAGTG  >  1/1089521‑1090003
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
cAGGTGGTTAATGAAGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATgg                                                                                                                                                                                                                                           <  1:894795/250‑1 (MQ=255)
   gtggtTAATGAAGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGcac                                                                                                                                                                                                                                        <  1:869420/250‑1 (MQ=255)
   gtggtTAATGAAGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGca                                                                                                                                                                                                                                         <  1:684635/249‑1 (MQ=255)
     ggtTAATGAAGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATg                                                                                                                                                                                                                                     >  1:345444/1‑251 (MQ=255)
             aaGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAg                                                                                                                                                                                                                              <  1:477497/250‑1 (MQ=255)
              aGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGca                                                                                                                                                                                                                                         <  1:76400/238‑1 (MQ=255)
               gTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGt                                                                                                                                                                                                                           >  1:14668/1‑251 (MQ=255)
                            gATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAAt                                                                                                                                                                                                              >  1:115227/1‑251 (MQ=255)
                            gATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAAt                                                                                                                                                                                                              >  1:499167/1‑251 (MQ=255)
                                          aTTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATggg                                                                                                                                                                                                   <  1:541178/248‑1 (MQ=255)
                                          aTTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAg                                                                                                                                                                                                 <  1:762328/250‑1 (MQ=255)
                                                           aaCTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGAt                                                                                                                                                                                                                       <  1:354817/211‑1 (MQ=255)
                                                             cTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGtt                                                                                                                                                                             >  1:391101/1‑251 (MQ=255)
                                                                  aTAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCg                                                                                                                                                                        <  1:776728/251‑1 (MQ=255)
                                                                           ccTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACatgat                                                                                                                                                               <  1:1308/251‑1 (MQ=255)
                                                                                  gCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGt                                                                                                                                                                                                                                                >  1:59820/1‑163 (MQ=255)
                                                                                   cGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGAtgggc                                                                                                                                                           <  1:610406/247‑1 (MQ=255)
                                                                                                          gTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTAAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAgg                                                                                                                                >  1:450171/1‑250 (MQ=255)
                                                                                                               gTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCa                                                                                                                                                                                                                                              >  1:306733/1‑136 (MQ=255)
                                                                                                                    ggATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACatgat                                                                                                                                                               <  1:834310/210‑1 (MQ=255)
                                                                                                                       tatTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGc                                                                                                                                                                                     <  1:560104/185‑1 (MQ=255)
                                                                                                                           ttGGTTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCg                                                                                                               <  1:388996/251‑1 (MQ=255)
                                                                                                                                tGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATg                                                                                                          <  1:163727/251‑1 (MQ=255)
                                                                                                                                         cGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 <  1:30307/251‑1 (MQ=255)
                                                                                                                                         cGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 <  1:143055/251‑1 (MQ=255)
                                                                                                                                           gAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACatgat                                                                                                                                                               <  1:226717/187‑1 (MQ=255)
                                                                                                                                           gAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGa                                                                                                <  1:871915/250‑1 (MQ=255)
                                                                                                                                               aaGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 <  1:749764/245‑1 (MQ=255)
                                                                                                                                                aGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGca                                                                                                                                                                                                                                         <  1:313585/108‑1 (MQ=255)
                                                                                                                                                aGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGt                                                                                          >  1:854172/1‑251 (MQ=255)
                                                                                                                                                aGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGt                                                                                          >  1:494215/1‑251 (MQ=255)
                                                                                                                                                 gAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGt                                                                                                                                                                                                                                                <  1:412536/100‑1 (MQ=255)
                                                                                                                                                 gAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGa                                                                                                <  1:80781/244‑1 (MQ=255)
                                                                                                                                                  aataaataGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGa                                                                                                            >  1:281174/1‑231 (MQ=255)
                                                                                                                                                        taGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCg                                                                                                               >  1:150058/1‑222 (MQ=255)
                                                                                                                                                        taGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 <  1:302500/236‑1 (MQ=255)
                                                                                                                                                              gggTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGtatat                                                                                                     <  1:175745/226‑1 (MQ=255)
                                                                                                                                                                gTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTc                                                                                                                                                                                                                                >  1:847041/1‑101 (MQ=255)
                                                                                                                                                                gTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGaaaaa                                                                          >  1:86766/1‑251 (MQ=255)
                                                                                                                                                                   ataAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAAcc                                                                                                                                                                         <  1:370602/153‑1 (MQ=255)
                                                                                                                                                                        gTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAg                                                                                                                         <  1:887119/196‑1 (MQ=255)
                                                                                                                                                                                 gATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACg                                                         >  1:623380/1‑251 (MQ=255)
                                                                                                                                                                                   ttGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 >  1:637235/1‑209 (MQ=255)
                                                                                                                                                                                          cGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCtt                                                <  1:314609/251‑1 (MQ=255)
                                                                                                                                                                                                aTTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGa                                                                                                >  1:431761/1‑197 (MQ=255)
                                                                                                                                                                                                aTTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCag                                          >  1:831778/1‑251 (MQ=255)
                                                                                                                                                                                                     gggggATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGTTTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGCGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAgaagaa                                     >  1:897099/1‑251 (MQ=255)
                                                                                                                                                                                                        ggATGGAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCgg                                                                                                                                     >  1:280257/1‑152 (MQ=255)
                                                                                                                                                                                                             gAATTAATTATACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGa                                                                                                <  1:660389/184‑1 (MQ=255)
                                                                                                                                                                                                                 taattaTACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAatat                          >  1:841388/1‑250 (MQ=255)
                                                                                                                                                                                                                 taattaTACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTAATGTATATAAAGATCCCGTTACCGGTGCGGaaaaa                                                                          <  1:489271/202‑1 (MQ=255)
                                                                                                                                                                                                                  aattaTACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGa                                                                                                                                                 >  1:659043/1‑130 (MQ=255)
                                                                                                                                                                                                                    ttaTACTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGG‑AAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTaa                      >  1:227441/1‑250 (MQ=255)
                                                                                                                                                                                                                     tataCTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAAc                     <  1:47011/251‑1 (MQ=255)
                                                                                                                                                                                                                       taCTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACtt                   <  1:548702/251‑1 (MQ=255)
                                                                                                                                                                                                                         cTTCATTGACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCtt                                                <  1:272291/220‑1 (MQ=255)
                                                                                                                                                                                                                                gACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGc                                                                                                                           <  1:351207/138‑1 (MQ=255)
                                                                                                                                                                                                                                gACACACATACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCgaag          >  1:566665/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cacacaTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGc                                                                                                                                       <  1:27742/124‑1 (MQ=255)
                                                                                                                                                                                                                                  cacacaTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 <  1:52122/162‑1 (MQ=255)
                                                                                                                                                                                                                                  cacacaTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGa                                                                                                >  1:553122/1‑163 (MQ=255)
                                                                                                                                                                                                                                  cacacaTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCgaagaa        >  1:763253/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cacaTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCagaag                                       <  1:592012/218‑1 (MQ=255)
                                                                                                                                                                                                                                      caTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGc                                                                                                                       <  1:404934/136‑1 (MQ=255)
                                                                                                                                                                                                                                      caTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACtt                   >  1:278716/1‑236 (MQ=255)
                                                                                                                                                                                                                                      caTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCGaagaagaag    >  1:893053/1‑251 (MQ=255)
                                                                                                                                                                                                                                      caTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTCGGCGATGAAGGCAAGTGCGGTGAGTATTGCAGGAAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCGaagaagaag    <  1:579623/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGa                                                                                                >  1:337303/1‑158 (MQ=255)
                                                                                                                                                                                                                                       aTACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCGAAGAAGAAGt   >  1:279909/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tACTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCGAAGAAGAAGTg  >  1:60606/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aCTGGATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAg                                                                                                 >  1:732400/1‑155 (MQ=255)
                                                                                                                                                                                                                                            ggATACTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCg                                                                                                                                                                        <  1:456228/81‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
CAGGTGGTTAATGAAGTTTGGGGGGAACGATTTAAGCAGCAGATTTCACAAATGTCGGGAACTTTGATAATTCGTCCTGACAGCGGTGAACCGGGCACGGTAATTTGTCGCGTACTGGATATTTTGGCTGAAAAGTTCGGAACAAGAATAAATAGTAAGGGTTATAAAGTATTGCCGGATTGCATTCGTGTCATTCAGGGGGATGGAATTAATTATACTTCATTGACACACATACTGGATGCTGTCATGGCACATGGATTCAGCGTGGATAATGTGAATTTCGGTATGGGAGGCGGCTTGTTGCAACAGGTTAACCGTGACATGATGGGCTGGGCGATGAAGGCAAGTGCGGTGAGTGTTGCAGGCAAATGGCGTGATGTATATAAAGATCCCGTTACCGGTGCGGAAAAACGTTCGAAAAAAGGACGGCTTGCCTTAGTCAGAAGAAATGGCGAATATCTAACTTTGCGCGAAGAAGAAGTG  >  1/1089521‑1090003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑