Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,1670GT27.5% 154.8 / 16.7 70L33F (TTG→TTTynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (10/40);  new base T (19/0);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 4.33e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTG  >  NC_000913/1657933‑1658404
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
aCCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTTCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGt                                                                                                                                                                                                                                                 <  2:42089/239‑1 (MQ=255)
 ccTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTcccgt                                                                                                                                                                                                                                                 >  1:2056791/1‑236 (MQ=255)
 ccTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGc                                                                                                                                                                                                                                                     >  2:2428178/1‑235 (MQ=255)
      ttaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTcttttcccttttcc                                                                                                                                                                                                                                            >  1:1556733/1‑238 (MQ=255)
              aCCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAAATTTTATATAAAGATAAAGAACAGGGAGTGAGTTATGTCAAAAAATGAACGAATGGTGGGCATCAGCCGCAGAAAACTCGTTAAATTTACCGCGATAGGTTATCTGGCGCTGGCGGCAGGCGGGTTTATTTTTGCCGTTTCCCTTgccaaa                                                                                                                                                                                                                                    >  2:1653539/1‑235 (MQ=255)
                   tCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTTCGCAATgcagc                                                                                                                                                                                                                              >  1:2175190/1‑239 (MQ=255)
                   tCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagc                                                                                                                                                                                                                              <  2:2299844/239‑1 (MQ=255)
                    cTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTACCTCCGCAAAgcagca                                                                                                                                                                                                                             >  2:547133/1‑239 (MQ=255)
                      aTTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  2:1404388/239‑1 (MQ=255)
                            atGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTTCGTTTTCCCCGTGCGATTCCGCAACACCggg                                                                                                                                                                                                                    >  2:420614/1‑238 (MQ=255)
                                 ttAATAATTCTTAACCCAAAATGGGTATATTCTCTCTATTTTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATATAGAACAGGGAGTGAGTTATGTACAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGGTCTCTGGCGCTGGCTGCAGGCGGCCTTTTCTTTCCGTTTACCCCGCGCACAGCATCAGCAGCAGtgtag                                                                                                                                                                                                                >  2:2091666/1‑234 (MQ=255)
                                     taatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGc                                                                                                                                                                                                           <  1:2452751/239‑1 (MQ=255)
                                         tCTTAACCCAAACTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGCTTACCCTTCGGCATTCAGCAGCAGCGGGACAACAGGCCcgc                                                                                                                                                                                                       >  2:672740/1‑239 (MQ=255)
                                                                 ccTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTgg                                                                                                                                                                                        <  1:1919845/230‑1 (MQ=255)
                                                                     tATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATACAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTCCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  1:2176100/192‑1 (MQ=255)
                                                                            aGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAATGGTCTGGGGTGGCTGGTCCGCCa                                                                                                                                                                       >  1:2032888/1‑235 (MQ=255)
                                                                             gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  1:1544207/184‑1 (MQ=255)
                                                                             gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           >  2:1544207/1‑184 (MQ=255)
                                                                               gcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAAAAGGGAGGGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTCAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTTGCGTTTTCCCTGTGGCAAGCAGCAGAAGCgggccaa                                                                                                                                                                                                                >  2:2176100/1‑187 (MQ=255)
                                                                                 gcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGt                                                                                                                                                               <  1:420614/238‑1 (MQ=255)
                                                                                           cAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGc                                                                                                                                                   <  1:2428178/239‑1 (MQ=255)
                                                                                             atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGt                                                                                                                                                                                                                                                 >  2:402335/1‑146 (MQ=255)
                                                                                             atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGGTTTACCTTCGCCATGCCACCGCAGCCGGGACAAACGACCGGCGAAAAAGGGGCTGGGGGGCGCGGTTCGGTAAACTGTGGTAGcgcacgcgc                                                                                                                                                    >  2:1574805/1‑230 (MQ=255)
                                                                                             atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGt                                                                                                                                                                                                                                                 <  1:402335/146‑1 (MQ=255)
                                                                                              tatataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACt                                                                                                                                                <  1:547133/239‑1 (MQ=255)
                                                                                                     aCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTa                                                                                                                                        <  1:1653539/239‑1 (MQ=255)
                                                                                                      cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTAc                                                                                                                                       <  1:1469231/239‑1 (MQ=255)
                                                                                                      cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTAc                                                                                                                                       <  1:1421499/239‑1 (MQ=255)
                                                                                                         ttATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGGTTACCCTTCGCCATGCATCAGCAGCGGGGCAACAAGCCCCCGGAAAAAAGGGCGTGGGGGGCTGTTCCGTCAACAGTGGGAGGCGCGTGTGCAGTTctctaacat                                                                                                                                      >  2:1714545/1‑231 (MQ=255)
                                                                                                             ataACGAAAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTaa                                                                                                                               <  1:1672083/239‑1 (MQ=255)
                                                                                                                 cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGcagcagcag                                                                                                                                                                                                                         >  1:1299974/1‑150 (MQ=255)
                                                                                                                 cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGcagcagcag                                                                                                                                                                                                                         <  2:1299974/150‑1 (MQ=255)
                                                                                                                 cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTTCCGCTTACCCCTCGCCAAGCAGCAGCAGCGGTAAAAACGGCCCGCGGAAAAGGGGGCTGGGGTGCCTGTTTCGTTTACCGGGGTAGGCGCTGTGTACTCTGTCTCTATTTTTAAGag                                                                                                                           >  2:2568681/1‑238 (MQ=255)
                                                                                                                       aGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                                     <  1:2182105/239‑1 (MQ=255)
                                                                                                                              ggAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCCGCGCAATGCAGCAGCAGCGGCACCACAAGCCCGCGGAAAAGTGGTCAGGGGGGGCCGTTCCGTGAACTATGGGAGCCGCCTTGCACTTTGTCTACCGGGTAAAGATAATGacatga                                                                                                                 >  2:1230672/1‑231 (MQ=255)
                                                                                                                               gagtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                                     <  1:1092478/231‑1 (MQ=255)
                                                                                                                               gagtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGAcc                                                                                                               <  1:2258841/237‑1 (MQ=255)
                                                                                                                                   gagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTggg                                                                                                           <  1:672740/237‑1 (MQ=255)
                                                                                                                                    gattATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGCACCCTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTgg                                                                                                        <  2:2056791/237‑1 (MQ=255)
                                                                                                                                    agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCcgcg                                                                                                                                                                                                      <  1:2388980/149‑1 (MQ=255)
                                                                                                                                      ttATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTgg                                                                                                                                                                                        <  1:724195/161‑1 (MQ=255)
                                                                                                                                        aTGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           >  1:1319955/1‑125 (MQ=255)
                                                                                                                                        aTGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  2:1319955/125‑1 (MQ=255)
                                                                                                                                         tGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCCGCAGCGGGACAACAAGTCCCCGGAAAAGTGGTCTGGGGGGGCCGTTCCGTCTACAGCGGTAGCCCCTGTTGCCCTCGTCTTACATTTAAAGTTACTGATGTGACCTggg                                                                                                           >  2:1092478/1‑231 (MQ=255)
                                                                                                                                          gTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGGAAAAAGGGTCTGGGGGGCCTTGTCCGGCACCGGGGGGAGCCCCCGGGGACTTCGGCTTCACGGTTAAAGATATTAAGGGGCCTTGGGGGaaaact                                                                                                  >  2:808677/1‑236 (MQ=255)
                                                                                                                                          gTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTg                                                                                                  <  1:1674831/239‑1 (MQ=255)
                                                                                                                                              aaaaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGGTTTCCCCTCCCCAAGCAGCCGCCGCGGGACAACAGGGCCCCGGAAAAATGGt                                                                                                                                                                                            >  2:2388980/1‑149 (MQ=255)
                                                                                                                                                    gAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCgg                                                                                        <  1:808677/239‑1 (MQ=255)
                                                                                                                                                      aCGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTgggg                                                                                                                                                                                      <  1:775690/147‑1 (MQ=255)
                                                                                                                                                             gtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGGAAAAGGGGGCTGGGGTGCCTGTTCCCTCAACTGTGGTAGCCGCTGTGCACTTCGTCGTTCAGTTTAAGATAATGAAGTGGCCCGGGGGGAAAAAGAGAAAAACGGCGGCGAGGa                                                                              >  2:373337/1‑239 (MQ=255)
                                                                                                                                                                ggCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTTCCGTTTAACCTTCTCAATGCAGCAGCAAGGGGACAAAAGGGCCGCGGAAAAGTGGTCCGGGGGGGCTGGTcc                                                                                                                                                                           >  2:775690/1‑147 (MQ=255)
                                                                                                                                                                          gCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCCGTTACCCCCCCCAATTCAGCCGCCCCGGTacaa                                                                                                                                                                                                                >  1:1180714/1‑101 (MQ=255)
                                                                                                                                                                                     gTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTaa                                                                                                                               <  1:2076496/167‑1 (MQ=255)
                                                                                                                                                                                          aTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGt                                                                                                                                                                              <  1:2530132/118‑1 (MQ=255)
                                                                                                                                                                                              aCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAGTGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTgcg                                             <  1:2227779/239‑1 (MQ=255)
                                                                                                                                                                                               ccGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTTCCTGCGCAAATCCGCAGCAGCGGGACAACAGGGCCGCGAAAAAATGGGCTGGGGGGCCCTGTCCGGTAACTGGGGGAGCCGCCGGGCACTTCGTCTTAATGTTTAAAATAATgta                                                                                                                     >  2:2076496/1‑165 (MQ=255)
                                                                                                                                                                                                 gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTaa                                                                                                                               <  1:1809806/155‑1 (MQ=255)
                                                                                                                                                                                                    aTAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGGTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGGCTGGGGGGGCTGTTCCGTCaaat                                                                                                                                                                    >  2:2530132/1‑116 (MQ=255)
                                                                                                                                                                                                       ggTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGtt                                                                                                                                 <  1:1329466/147‑1 (MQ=255)
                                                                                                                                                                                                         ttCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGGTTTACCCGCCCCATTCAGCAGCAGCGGTGCAACAAGCCCCCCGAAAAGGGGTCTGGGGGGCCTGTTCCGTCAAATGGGGTAGCCGCCGTGCCCTTCGGCTGCATGTTAAAGATGATGAAATGACCTGGGTGGTAACTGACCAAACCGGCAGCCGACGGTAATGCAACCAACCGGAGCGCGCCTGTGTGTACGGTCGCttc                                   >  2:146597/1‑236 (MQ=255)
                                                                                                                                                                                                           ctctGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGGAAAAGTGGGCTGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCCGTGCACTTCGTCTCCATGTTAAAGATAATGaa                                                                                                                     >  2:1809806/1‑155 (MQ=255)
                                                                                                                                                                                                              tGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgc                             <  1:373337/239‑1 (MQ=255)
                                                                                                                                                                                                                     gtctgCAGGCGTGTTGTCTTGGCCGTTTACCCTGGGCAATGCAGCAGCAGAGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt                       <  2:376755/236‑1 (MQ=255)
                                                                                                                                                                                                                        tgCAGTCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTa                   <  2:2175190/239‑1 (MQ=255)
                                                                                                                                                                                                                        gtccggCGTTTTGTCTTTGCCGTTTACCCTGCGCAATGCAGCGGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTa                   <  2:1457113/235‑1 (MQ=255)
                                                                                                                                                                                                                          cAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGtt                                                                                                                                 <  1:2406794/128‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCATCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGGCCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCACTCCATCCGCCGGCGTATTAATCATCCCGATc       <  1:2091666/239‑1 (MQ=255)
                                                                                                                                                                                                                                        cTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGc     <  1:376280/239‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAg                                                                                                                                                             <  1:958034/88‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAg                                                                                                                                                             >  2:958034/1‑88 (MQ=255)
                                                                                                                                                                                                                                           tGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTg  >  1:877590/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
ACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTG  >  NC_000913/1657933‑1658404

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