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→TTT) ynfEputative 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.

ACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTG  >  NC_000913/1657933‑1658404
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
aCCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTTCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑Gt                                                                                                                                                                                                                                                 <  2:42089/239‑1 (MQ=255)
 ccTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTcccgt                                                                                                                                                                                                                                                 >  1:2056791/1‑236 (MQ=255)
 ccTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGc                                                                                                                                                                                                                                                     >  2:2428178/1‑235 (MQ=255)
      ttaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTcttttcc‑cttttcc                                                                                                                                                                                                                                            >  1:1556733/1‑238 (MQ=255)
              aCCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAAATTTTATATAAAGATAAAGAACAGGGAGTGAGTTATGTCAAAAAATGAACGAATGGTGGGCATCAGCCGCAGAAAACTCGTTAAATTTACCGCGATAGGTTATCTGGCGCTGGCGGCAGGCGGG‑TTTATTTTTGCC‑GTTTCCCTTgccaaa                                                                                                                                                                                                                                    >  2:1653539/1‑235 (MQ=255)
                   tCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑GTTTACCCTTCGCAATgcagc                                                                                                                                                                                                                              >  1:2175190/1‑239 (MQ=255)
                   tCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATgcagc                                                                                                                                                                                                                              <  2:2299844/239‑1 (MQ=255)
                    cTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑GTTTTACCTCCGCAAAgcagca                                                                                                                                                                                                                             >  2:547133/1‑239 (MQ=255)
                      aTTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  2:1404388/239‑1 (MQ=255)
                            atGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTTC‑GTTTTCCCCGTGCGATTCCGCAACAC‑Cggg                                                                                                                                                                                                                    >  2:420614/1‑238 (MQ=255)
                                 ttAATAATTCTTAACCCAAAATGGGTATATTCTCTCTATTTTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATATAGAACAGGGAGTGAGTTATGTACAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGGTCTCTGGCGCTGGCTGCAGGC‑GGCCTTTTCTTT‑CC‑GTTTACCCCGCGCACAGCATCAGCAG‑CAGtgtag                                                                                                                                                                                                                >  2:2091666/1‑234 (MQ=255)
                                     taatTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGc                                                                                                                                                                                                           <  1:2452751/239‑1 (MQ=255)
                                         tCTTAACCCAAACTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTGCC‑GCTTACCCTTCGGCATTCAGCAGCAG‑CGGGACAACAGGCCcgc                                                                                                                                                                                                       >  2:672740/1‑239 (MQ=255)
                                                                 ccTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTgg                                                                                                                                                                                        <  1:1919845/230‑1 (MQ=255)
                                                                     tATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATACAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTCCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  1:2176100/192‑1 (MQ=255)
                                                                            aGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTTCCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAATGGTCTGGGGTGG‑CTGGTCCGCCa                                                                                                                                                                       >  1:2032888/1‑235 (MQ=255)
                                                                             gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  1:1544207/184‑1 (MQ=255)
                                                                             gcgcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           >  2:1544207/1‑184 (MQ=255)
                                                                               gcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAAAAGGGAGGGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTCAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑ATTTTTTTTTGC‑GTTTTCCCTGTGGCAAGCAGCAGAAG‑Cgggccaa                                                                                                                                                                                                                >  2:2176100/1‑187 (MQ=255)
                                                                                 gcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGt                                                                                                                                                               <  1:420614/238‑1 (MQ=255)
                                                                                           cAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGc                                                                                                                                                   <  1:2428178/239‑1 (MQ=255)
                                                                                             atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTGCC‑Gt                                                                                                                                                                                                                                                 >  2:402335/1‑146 (MQ=255)
                                                                                             atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GGTTTACCTTCGCCATGCCACCGCAGCCGGGACAAACGACCGGCGAAAAAG‑GGGCTGGGGGGCGCGGTTCGGTAAACTGTGGTAGcgcacgcgc                                                                                                                                                    >  2:1574805/1‑230 (MQ=255)
                                                                                             atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑Gt                                                                                                                                                                                                                                                 <  1:402335/146‑1 (MQ=255)
                                                                                              tatataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACt                                                                                                                                                <  1:547133/239‑1 (MQ=255)
                                                                                                     aCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑Ta                                                                                                                                        <  1:1653539/239‑1 (MQ=255)
                                                                                                      cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TAc                                                                                                                                       <  1:1469231/239‑1 (MQ=255)
                                                                                                      cTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TAc                                                                                                                                       <  1:1421499/239‑1 (MQ=255)
                                                                                                         ttATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GGTTACCCTTCGCCATGCATCAGCAG‑CGGGGCAACAAGCCCCCGGAAAAAAGGGCGTGGGGGG‑CTGTTCCGTCAACAGTGGGAGGCGCGTGTGCAGTTctctaacat                                                                                                                                      >  2:1714545/1‑231 (MQ=255)
                                                                                                             ataACGAAAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTaa                                                                                                                               <  1:1672083/239‑1 (MQ=255)
                                                                                                                 cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GTTTACCCTGCGCAATGcagcagcag                                                                                                                                                                                                                         >  1:1299974/1‑150 (MQ=255)
                                                                                                                 cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGcagcagcag                                                                                                                                                                                                                         <  2:1299974/150‑1 (MQ=255)
                                                                                                                 cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑ATTTTTTTTTCC‑GCTTACCCCTCGCCAAGCAGCAGCAG‑CGGTAAAAACGGCCCGCGGAAAAGGGGGCTGGGGTGC‑CTGTTTCGTTTACCGGGGTAGGCGC‑TGTGTACTCTGTC‑TCTA‑TTTTTAAGag                                                                                                                           >  2:2568681/1‑238 (MQ=255)
                                                                                                                       aGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGaa                                                                                                                     <  1:2182105/239‑1 (MQ=255)
                                                                                                                              ggAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GTTTACCCCGCGCAATGCAGCAGCAG‑CGGCACCACAAGCCCGCGGAAAAGTGGTCAGGGGGGG‑CCGTTCCGTGAACTATGGGAGCCGC‑CTTGCACTTTGTC‑TACC‑GGGTAAAGATAATGacatga                                                                                                                 >  2:1230672/1‑231 (MQ=255)
                                                                                                                               gagtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGaa                                                                                                                     <  1:1092478/231‑1 (MQ=255)
                                                                                                                               gagtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGAcc                                                                                                               <  1:2258841/237‑1 (MQ=255)
                                                                                                                                   gagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTggg                                                                                                           <  1:672740/237‑1 (MQ=255)
                                                                                                                                    gattATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGCACCCTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTgg                                                                                                        <  2:2056791/237‑1 (MQ=255)
                                                                                                                                    agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCcgcg                                                                                                                                                                                                      <  1:2388980/149‑1 (MQ=255)
                                                                                                                                      ttATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTgg                                                                                                                                                                                        <  1:724195/161‑1 (MQ=255)
                                                                                                                                        aTGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           >  1:1319955/1‑125 (MQ=255)
                                                                                                                                        aTGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATgcagcagc                                                                                                                                                                                                                           <  2:1319955/125‑1 (MQ=255)
                                                                                                                                         tGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GTTTACCCTGCGCAATGCAGCCGCAG‑CGGGACAACAAGTCCCCGGAAAAGTGGTCTGGGGGGG‑CCGTTCCGTCTACAGCGGTAGCCCC‑TGTTGCCCTCGTCTTACA‑T‑TTAAAGTTACTGATGTGACCTggg                                                                                                           >  2:1092478/1‑231 (MQ=255)
                                                                                                                                          gTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGGACAACAGGCCCGCGGAAAAAGGGTCTGGGGGGC‑CTTGTCCGGCACCGGGGGGAGCCCC‑CGGGGACTTCGGC‑TTCACGGTTAAAGAT‑ATTAAGGGGCCTTGGGGGaaaact                                                                                                  >  2:808677/1‑236 (MQ=255)
                                                                                                                                          gTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTg                                                                                                  <  1:1674831/239‑1 (MQ=255)
                                                                                                                                              aaaaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTGCCGGTTTCCCCTCCCCAA‑GCAGCCGCCG‑CGGGACAACAGGGCCCCGGAAAAATGGt                                                                                                                                                                                            >  2:2388980/1‑149 (MQ=255)
                                                                                                                                                    gAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCgg                                                                                        <  1:808677/239‑1 (MQ=255)
                                                                                                                                                      aCGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTgggg                                                                                                                                                                                      <  1:775690/147‑1 (MQ=255)
                                                                                                                                                             gtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGGAAAAGGGGGCTGGGGTGC‑CTGTTCCCTCAACTGTGGTAGCCGC‑TGTGCACTTCGTCGTTCA‑‑GTTTAAGATAATGAAGTGGCCCGGGGGGAAAAAGAGAAAAACGGCGG‑CGAGGa                                                                              >  2:373337/1‑239 (MQ=255)
                                                                                                                                                                ggCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑ATTTTTTTTTCC‑GTTTAACCTTCTCAATGCAGCAGCAA‑GGGGACAAAAGGGCCGCGGAAAAGTGGTCCGGGGGGG‑CTGGTcc                                                                                                                                                                           >  2:775690/1‑147 (MQ=255)
                                                                                                                                                                          gCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑CGTTACCCCCCCCAATTCAGCCGCCC‑CGGTacaa                                                                                                                                                                                                                >  1:1180714/1‑101 (MQ=255)
                                                                                                                                                                                     gTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTaa                                                                                                                               <  1:2076496/167‑1 (MQ=255)
                                                                                                                                                                                          aTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGt                                                                                                                                                                              <  1:2530132/118‑1 (MQ=255)
                                                                                                                                                                                              aCCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAGTGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTgcg                                             <  1:2227779/239‑1 (MQ=255)
                                                                                                                                                                                               ccGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTTCC‑GTTTTTCCTGCGCAAATCCGCAGCAG‑CGGGACAACAGGGCCGCGAAAAAATGGGCTGGGGGGC‑CCTGTCCGGTAACTGGGGGAGCCGC‑CGGGCACTTCGTC‑TTAA‑TGTTTAAAATAATgta                                                                                                                     >  2:2076496/1‑165 (MQ=255)
                                                                                                                                                                                                 gcgATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTaa                                                                                                                               <  1:1809806/155‑1 (MQ=255)
                                                                                                                                                                                                    aTAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑GGTTACCCTGCGCAATGCAGCAGCAG‑CGGGACAACAGGCCCGCGAAAAAGTGGGCTGGGGGGG‑CTGTTCCGTCaaat                                                                                                                                                                    >  2:2530132/1‑116 (MQ=255)
                                                                                                                                                                                                       ggTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGtt                                                                                                                                 <  1:1329466/147‑1 (MQ=255)
                                                                                                                                                                                                         ttCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCCGGTTTACCC‑GCCCCATTCAGCAGCAG‑CGGTGCAACAAGCCCCCCGAAAAGGGGTCTGGGGGGC‑CTGTTCCGTCAAATGGGGTAGCCGC‑CGTGCCCTTCGGC‑TGCA‑TGTTAAAGATGATGAAATGACCTGGGTGGTAACTGACCAAACCGGCAGCCGACG‑GTAATGCAACCAACCGGAGCGCGCCTGTGTGTACGGTCGCttc                                   >  2:146597/1‑236 (MQ=255)
                                                                                                                                                                                                           ctctGGCGCTGGCTGCAGGC‑GG‑TTTTTTTTTTCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGGAAAAGTGGGCTGGGGGGC‑CTGTTCCGTCAACTGTGGGAGCCGC‑CGTGCACTTCGTC‑TCCA‑TGTTAAAGATAATGaa                                                                                                                     >  2:1809806/1‑155 (MQ=255)
                                                                                                                                                                                                              tGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccgc                             <  1:373337/239‑1 (MQ=255)
                                                                                                                                                                                                                     gtctgCAGGC‑GT‑GTTGTCTTGGCC‑GTTTACCCTGGGCAATGCAGCAGCAG‑AGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt                       <  2:376755/236‑1 (MQ=255)
                                                                                                                                                                                                                        tgCAGTC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTa                   <  2:2175190/239‑1 (MQ=255)
                                                                                                                                                                                                                        gtccggC‑GT‑TTTGTCTTTGCC‑GTTTACCCTGCGCAATGCAGCGGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTa                   <  2:1457113/235‑1 (MQ=255)
                                                                                                                                                                                                                          cAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGtt                                                                                                                                 <  1:2406794/128‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCATCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGGCCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCACTCCATCCGCCGGCGTATTAATCATCCCGATc       <  1:2091666/239‑1 (MQ=255)
                                                                                                                                                                                                                                        cTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGc     <  1:376280/239‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAg                                                                                                                                                             <  1:958034/88‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAg                                                                                                                                                             >  2:958034/1‑88 (MQ=255)
                                                                                                                                                                                                                                           tGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTg  >  1:877590/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
ACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGC‑GG‑TTTTTCTTTGCC‑GTTTACCCTGCGCAATGCAGCAGCAG‑CGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGC‑CTGTTCCGTCAACTGTGGTAGCCGC‑TGTGCACTTCGTC‑TACA‑TGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTG  >  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: ‑