Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,972,0900AC28.2% 135.9 / 12.1 79V201G (GTG→GGG) tarmethyl‑accepting chemotaxis protein Tar
Reads supporting (aligned to +/- strand):  ref base A (41/15);  new base C (0/22);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 6.61e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATATTACGACTGGTAGC  >  NC_000913/1971841‑1972340
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          
gccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATaccac                                                                                                                                                                                                                                                           <  2:207069/251‑1 (MQ=255)
   gcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACAACCGGCGATTTCGCGAATGTGAGCAATACTTTTTGCCAGCGGAGTAAGCAACATACGGCGAAGGCCGTACCACGCCCACAGCAGAATAAATaccacca                                                                                                                                                                                                                                                         >  2:498418/1‑250 (MQ=255)
    cAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAGTGCCGTACCACGCCACCAGCAGAATCAATaccaccacc                                                                                                                                                                                                                                                       <  1:319030/251‑1 (MQ=255)
              ggTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATaa                                                                                                                                                                                                                                             <  1:349782/251‑1 (MQ=255)
               gTACCGGCATAGTTGGCATCTGAACCTTCGCGGACATGAGTGGCGGTGTCAGACAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAAc                                                                                                                                                                                                                                            <  2:7546/251‑1 (MQ=255)
               gTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAAc                                                                                                                                                                                                                                            <  2:707197/251‑1 (MQ=255)
                 aCCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCg                                                                                                                                                                                                                                          >  1:697449/1‑251 (MQ=255)
                  ccGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGAAGAATCAATACCACCACCAGCGCGTTAACCGc                                                                                                                                                                                                                                         >  1:470489/1‑251 (MQ=255)
                   gggCATAAAGGGCATCTGAACCTTCGCGGACATGAGGGACGGTGTCAGTCAAAGAGCGTTGCATATGGTAAACGCTCTGCGCCAGGTCGCCCATTGCACTGCGCCCGTCAATGGTCAGGGTATTCGCCTGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCCCGCCACCAGCAGAAGCAATACCCCCACCAGCGCGATAACCGcc                                                                                                                                                                                                                                        <  2:298348/250‑1 (MQ=255)
                    ggaatagagggattCTGAACCTTCGGGGACAGTGGTGCCGGGGCCGGTAAAAGAGCGTTGCAGATGGTAACCCCCGTGCCCCAGGCCGCCCATTTCCCGGCGCCCTGAAATGGTCAGGGTATTCGCCAGGTAACCACCGGCGATTTCGCGAAGGGGAGAAATAATTTTTGCCAGCGGAGTAAGCAACATACGCCGAAGGCCGTCCCACGCCACCACCAGAATCAATACCCCCACCAGCGCGATAACCGCCa                                                                                                                                                                                                                                       <  2:431875/238‑1 (MQ=255)
                     gCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGcc                                                                                                                                                                                                                                        >  1:145511/1‑249 (MQ=255)
                       aTAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGtt                                                                                                                                                                                                                                    >  2:322765/1‑251 (MQ=255)
                         aGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCTTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGTTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAACCGCCAGTTGc                                                                                                                                                                                                                                  <  1:655607/251‑1 (MQ=255)
                            tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCAc                                                                                                                                                                                                                               >  1:150921/1‑251 (MQ=255)
                             ggCATGTGAACCTGCGCGGAAAGGGGTGACGGGGTGAGTCAAAGTGCGTTTCATATGTGAATCCCCGTGCGCCAGGCCGCCCTTCTCCCTGCCCCCGTAAATGGTCAGGGTATTCGCCAGGTTACCCCCGGCAATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGTAGTAAGCAACATACGGCGAATGCCGTCCCACGCCACCAGCAAAATCAAGCCCCCCACCAGCGCGATAACCGCCAGTTGCCACt                                                                                                                                                                                                                              <  2:525084/251‑1 (MQ=255)
                             ggCATCTGAACCTGCGAGGACATGAGTGACGGGGTCAGTCAAAGAGCGTTGCATAGGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTAAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGAAATAATTTTTGCCAGGGGAGTAAGCAACATACGGAGAAGGCCGTACCACGCCACCAGCAGAATCAATCCCCCCACCAGCGCGATAACCGCCAGTTGCCACt                                                                                                                                                                                                                              <  1:380674/251‑1 (MQ=255)
                                     aaCCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAAt                                                                                                                                                                                                                      <  1:294334/251‑1 (MQ=255)
                                              ggACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAatc                                                                                                                                                                                                              >  1:326493/1‑250 (MQ=255)
                                                    ggAGTGACGGTTTCAGTCAAAGAGCGTTGCATATGTGAATCGCCCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCAAGGTTACCACCGGCGATGTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATCCCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGc                                                                                                                                                                                                        <  1:498418/250‑1 (MQ=255)
                                                     gAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAACCGCTCTGCGCCAGGTCCCCCATTTCCCTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCCCCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGc                                                                                                                                                                                                        <  1:846749/250‑1 (MQ=255)
                                                       gTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTCATCATCTGCGtt                                                                                                                                                                                                     >  1:354692/1‑251 (MQ=255)
                                                             gtgtCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGGCAATGGTAAGGGTATTCGCCAGGTTACCACCGGCGCTTTAGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGAGATACCAGCCAGTTTCCACTGGGCAACTCGGTACTCATATGCGTTTTCAGt                                                                                                                                                                                               >  2:768971/1‑251 (MQ=255)
                                                              tgtCAGTCAAAGAGCGATGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGACCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAAAATACGGAGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCAAGTTGAAACTGGGCAACTCGGTAATCATCTGCGTTGTCAGt                                                                                                                                                                                               >  1:464578/1‑250 (MQ=255)
                                                               gtcagtaaAAGAGCGTTGCATACGTGAACCGCCCTGCGCAAGGGCGCCCAGTTCCCTGCGCCCGTAAATGGTCAGGGTATCCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGAAACATACGGCGAATGCCGTCCCCCGCCCCCAGCAGATTAAATACCCCCACCAGCGCGATACCCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGa                                                                                                                                                                                             <  1:471425/251‑1 (MQ=255)
                                                                     aaaaGAGCTGTGCATATGGGAACCGCTCTGCGCCAGGTCGCCATTTCCCCTCCGCCTGTCAATGGTCAGGTAATGCGCCAGGTTACCACCGGCGATTTCGCGAAGGGGAGCAATAATTTTTGCCAGGGGAGTAAGCAACATCGGGCAAATGCCGTCCCCCGCCACCAGCAGAATCAATACCACCACCAGCGTGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGatat                                                                                                                                                                                       <  2:294135/250‑1 (MQ=255)
                                                                      ccttttccgatttaatggaaaaCGCTCTGCGCCAGGCCGCCCTTTTCACGGCGCCCGACAAGGTGCAGTTTATTCGCCAGGTTCCCACCGGCGATTTCGCGAAGGGGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCCCCAGCAGAATCAATCCCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATc                                                                                                                                                                                      <  2:307547/232‑1 (MQ=255)
                                                                         gagCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATcgcg                                                                                                                                                                                   >  1:277021/1‑251 (MQ=255)
                                                                             gTTGCATATGTGAAACGCTTGGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGTAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATAc                                                                                                                                                                               <  1:752634/251‑1 (MQ=255)
                                                                                   tatGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGTTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATc                                                                                                                                                                                      >  1:307547/1‑238 (MQ=255)
                                                                                        gaaccGCTCTGCGCCAGTTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATCCCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCAc                                                                                                                                                                     <  1:501676/246‑1 (MQ=255)
                                                                                        gAAACGCTCTGCGCCAGGCCGCCCATTTCTCTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCCCCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCAct                                                                                                                                                                    <  1:7213/251‑1 (MQ=255)
                                                                                         aaaCGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCActg                                                                                                                                                                   >  1:299972/1‑251 (MQ=255)
                                                                                             gCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCActgctg                                                                                                                                                                >  1:546015/1‑250 (MQ=255)
                                                                                               tctGCGCCAGTTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTAAGCAATCATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAAACAATACCACCACCAGCTCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAAGCATCTGAGCTATCAGTGACGCTCTCGCGATAAAGTTTTTCACTGCTAAgg                                                                                                                                                             >  1:35566/1‑251 (MQ=255)
                                                                                                   gcgcCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGCTAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGAc                                                                                                                                                                                            >  2:158963/1‑217 (MQ=255)
                                                                                                        aGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATACCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCTCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTg                                                                                                                                                     >  1:307397/1‑251 (MQ=255)
                                                                                                         ggTCGCCCATTTAACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGGTACCACCGGCGATTTCGCGAATGTGAGAAATAATTTTTGCCAGCGGAGTAAGCAACATATGGCGAATGCCGTACCACGCCCACAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTTCCACTGGGCAAATCGGTAATCATCTGCGTTATAAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTg                                                                                                                                                     >  2:871706/1‑250 (MQ=255)
                                                                                                          gccGCCCATTTCACGGCCCCGGTCAATGGTCGGGGTATTCGCCAGGTTACCACCGGCATTTGCGCGAATGGGAGCAAAAATTTTTGCCGGCGGATAAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAg                                                                                                                                                   <  2:354166/249‑1 (MQ=255)
                                                                                                          gTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAAGGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAg                                                                                                                                                   <  1:50347/251‑1 (MQ=255)
                                                                                                              cccATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCaaa                                                                                                                                               <  1:292956/251‑1 (MQ=255)
                                                                                                               ccATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAAc                                                                                                                                              >  1:209870/1‑251 (MQ=255)
                                                                                                                     cccTGCGCCCTTCAATGGCAAGGGTATTCGCCAGGTTACCACCGGCGATTCCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATCCCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTc                                                                                                                                         <  1:687838/248‑1 (MQ=255)
                                                                                                                     cACTGCTCCCGCTAATGGTCAGGTTATTCGCCAGGTGTCCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATCCGGCAAAGGCCGTACCACTCCACCAGCAGAATCAATCCCCCCCCCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTc                                                                                                                                         <  1:841487/250‑1 (MQ=255)
                                                                                                                            cgcccgcaaatgggaggggatgtCGCCAGGTTACCACCGGCGATTTCGCGAATGGGAGCAAAAATTTTGCCCTGGGTAGAAAGAAACATCCGGCGAATGCCGTACCCCGCCCCCAGCAGAACCAATACCACCACCAGCGGGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCAt                                                                                                                                   <  2:398823/229‑1 (MQ=255)
                                                                                                                             cccGTCAATGGTAAGGGTATTCGCCAGGTTACCCCCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCAttg                                                                                                                                 >  1:287245/1‑251 (MQ=255)
                                                                                                                                        tCAGGGGATTCGCCAGGATACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACAACCAGAGCGATAACCGCCAGTTGCCACTAGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGTCGTACTACTCAAACGCTTCGCCCAttgcattttgcatt                                                                                                                      >  1:668558/1‑251 (MQ=255)
                                                                                                                                               tATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATAAATACCACCACCAGCGCGATAACCGCCAGTTGACACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGCTATAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTgg                                                                                                                >  2:178221/1‑251 (MQ=255)
                                                                                                                                                 ttCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGt                                                                                                              >  1:655985/1‑251 (MQ=255)
                                                                                                                                                  tCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAACACGGTAATCATCTGATTTGTCAGTGACGATATCGAGCTAACGTTTTTCCTTGCTGGGTGCGTACTGAGCAAACGCTTCTCCCATTGGCTTTTTCAATCCCTGGGtt                                                                                                             >  2:250174/1‑251 (MQ=255)
                                                                                                                                                      cAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATACCGTACCACGCCACCACCAGAATCAATACCACCACCAGCGAGATAACAACCAGATGCCACTGGGCAAATCGGTAATCATCTACGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCATCGCCCACTGCATTATGCATCTCCTGGGtt                                                                                                             >  1:58375/1‑247 (MQ=255)
                                                                                                                                                       aGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTg                                                                                                        >  1:208185/1‑251 (MQ=255)
                                                                                                                                                         gTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTg                                                                                                        >  1:239952/1‑249 (MQ=255)
                                                                                                                                                              cctccggcgttatcggaaatgtggtcAATAATTTTTCCAGGGGGAGTAAGCACCCTGCGGCAAATGTCGTACCCCGCCCCCAGCATAAACACTCCCCCCACCAGCGCGATAACCGCCAGTTGCCACTGTGCAAATCGGTAATCATCTGCGTTGTCAGGGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGaa                                                                                                  <  2:837696/226‑1 (MQ=255)
                                                                                                                                                                 cGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATa                                                                                               >  1:300563/1‑250 (MQ=255)
                                                                                                                                                                           cggaaaggtgAGCAATAATTTTTGCCAGGGGAGTAAGCAACATCCGGCAAATGCCGTCCCCCCCAACCAGCAGAATCAATACCACCACCAGCGCGATATCCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTa                                                                                     <  2:150921/244‑1 (MQ=255)
                                                                                                                                                                             cgAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCCCCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTAtt                                                                                   <  1:863984/251‑1 (MQ=255)
                                                                                                                                                                             cgAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTAtt                                                                                   >  2:47966/1‑251 (MQ=255)
                                                                                                                                                                              aaaggtgAGCAATAATTTGTACCAGCGAGGTAAGCAACATACGCGGAAGGCCGAACCACGCCCCAAACAGAATCAATCCCCCCCCCTGCGGGATACCCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTg                                                                                  <  2:715240/247‑1 (MQ=255)
                                                                                                                                                                                        aataatTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCt                                                                         >  2:805134/1‑250 (MQ=255)
                                                                                                                                                                                          aaatTTTTGCCAGCGGAGTAAGCACCATACGGCAAAGGCCGTACCACGCCACCAGCAGAATCAAGACCACCACCAGCTCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATTTAGa                                                                      <  1:463704/250‑1 (MQ=255)
                                                                                                                                                                                            atTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGAt                                                                     >  1:499242/1‑250 (MQ=255)
                                                                                                                                                                                             tttttGCCAGCGGCGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATa                                                                    >  1:256422/1‑250 (MQ=255)
                                                                                                                                                                                             tttttGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATa                                                                    >  1:279176/1‑250 (MQ=255)
                                                                                                                                                                                                ttGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGCGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGCAATAAGATCCAGTATTGCCATAATCTAGATaatc                                                                 >  2:656311/1‑250 (MQ=255)
                                                                                                                                                                                                    cAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATaatcaata                                                             >  2:130594/1‑249 (MQ=255)
                                                                                                                                                                                                     aGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATaatcaatc                                                             >  1:731428/1‑249 (MQ=255)
                                                                                                                                                                                                              aaGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGGATAAATACCAACACCAGCGCGATAACCACCACTAGCCACTGGGCAAATCGGTAATCATCCGCGTTGTCAGTGACGATATCGCGATCAAGTTTTTCACTGCTGATGACGTACTGAGCAAACGCATCGCCCATCGCATTTTGCATTCCCTGGGTTGGCTGAGCGGAATAAGCTCCAGAATTGCCCTAATCTAGATaatcaatacgttatgtt                                                    >  2:778155/1‑240 (MQ=255)
                                                                                                                                                                                                                 aaccctccGGCAAATGCCGTCCCACCCCACCCCGAAAATCAAGACCCCCCCCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACg                                                <  2:326493/244‑1 (MQ=255)
                                                                                                                                                                                                                     aTACGGCGAATGCTGCACCACGCCCCAGTCAGAATAAATCCCACCACCAGCGCGATCACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATATTCTAGATAATCAATCAGTTCAGTTAACGCtgt                                            <  1:441473/251‑1 (MQ=255)
                                                                                                                                                                                                                     aTACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCtgt                                            >  1:651943/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggCGAATGCCGTCCCACGCCTCCAGCAGAATCAATCCCACCCCCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTaa                                        <  1:231786/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTa                                         >  1:705637/1‑250 (MQ=255)
                                                                                                                                                                                                                          gagattgcctacccACGCCACCAGCAGATTCAATACCTCCACCAGGGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAAt                                       <  1:863260/239‑1 (MQ=255)
                                                                                                                                                                                                                          gCGACTGCCGTCCCCCCCCACCAGCAGAATCAATACCCCCCCCAGCGCTATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAAt                                       <  1:2889/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCCGTCCCACGCCCCCAGCAAAATCAATACCACCCCCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGttt                                  <  1:49148/250‑1 (MQ=255)
                                                                                                                                                                                                                                      ccacccccccaGCAGAATCAATACCACACCCATCGTGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATAtt                           <  1:306600/243‑1 (MQ=255)
                                                                                                                                                                                                                                      ccACGCCACCAGCAGAATCAATACCCCCCCCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATAtt                           <  1:805134/251‑1 (MQ=255)
                                                                                                                                                                                                                                             accaGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTtcatc                    >  2:734293/1‑251 (MQ=255)
                                                                                                                                                                                                                                                agcagAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGc                                                                                 >  1:236860/1‑187 (MQ=255)
                                                                                                                                                                                                                                                agcagAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAat                 <  1:672767/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       aacAATACCACCACCCGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATATTACg           <  2:145511/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGCTATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGATATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTACCGCTGTGTAATAGTTTTTATATTTTTCATCAATATTACg           >  2:304656/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             accaccaccAGCGCGATAAACGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACCGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATATTACGACTGGt     >  1:380306/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                accaccaGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATATTACGACTGGTAGc  >  2:807000/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          
GCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGCCGTACCACGCCACCAGCAGAATCAATACCACCACCAGCGCGATAACCGCCAGTTGCCACTGGGCAAATCGGTAATCATCTGCGTTGTCAGTGACGATATCGCGATACAGTTTTTCACTGCTGAGGGCGTACTGAGCAAACGCTTCGCCCATTGCATTTTGCATTCCCTGGGTTGGCTGAGCGAAATAAGCTCCAGTATTGCCATAATCTAGATAATCAATCAGTTCAGTTAACGCTGTGTAATAGTTTTTATATTTTTCATCAATATTACGACTGGTAGC  >  NC_000913/1971841‑1972340

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