Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913728,0240TA22.4% 136.8 / 12.7 67N237Y (AAT→TAT) kdpAK(+) transporting P‑type ATPase subunit KdpA
Reads supporting (aligned to +/- strand):  ref base T (33/19);  new base A (0/15);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 8.37e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGATGCTGCTGTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAA  >  NC_000913/727776‑728270
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
atTGATGCTGCTGTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTaa                                                                                                                                                                                                                                                       <  1:130188/251‑1 (MQ=255)
  tGATGCTGCTGTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑GTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTaa                                                                                                                                                                                                                                                       >  2:497418/1‑249 (MQ=255)
   gATGCTGCTGTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGTCAG‑CATCTGCACGAGGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAaag                                                                                                                                                                                                                                                     <  1:184299/250‑1 (MQ=255)
    aTGCTGATGTCCGAGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCAAGCCTCAGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGCCAAAGCACAGCGCCGTTGGGATCAAGAAGC‑TCGCCAG‑CACCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAaaga                                                                                                                                                                                                                                                    >  2:401421/1‑250 (MQ=255)
         gctgTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGcc                                                                                                                                                                                                                                                <  2:501032/249‑1 (MQ=255)
           tgTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTCGCGG‑CGATCGCACATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGC‑TCGTCAG‑CATCTGCACGAAGTTGGTCAGTGAGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCAccc                                                                                                                                                                                                                                            >  1:127161/1‑250 (MQ=255)
            gTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCAccgc                                                                                                                                                                                                                                           >  1:164473/1‑251 (MQ=255)
                       gTGCCAGCAGATGCGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGACGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCCGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTACAGAAGCCACCGCCGttcgaaccg                                                                                                                                                                                                                                >  2:70493/1‑244 (MQ=255)
                        tGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCg                                                                                                                                                                                                                               >  1:390852/1‑250 (MQ=255)
                               agATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCtt                                                                                                                                                                                                                       >  1:78967/1‑251 (MQ=255)
                               agATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCAAGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGCCGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAgaatatt                                                                                                                                                                                                                       >  2:270446/1‑245 (MQ=255)
                                       aTTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGACGA‑TCGCCAG‑CATCTGCACGAAGTTGGGCCGTGCGGTTGGG‑TTTTCCAA‑CGG‑GTGCCACGAGTAGGCATTAAAGAAGCCACAGCCGTTAGTAACTGA‑CATCTTGATCACtt                                                                                                                                                                                                               >  2:238365/1‑251 (MQ=255)
                                       aTTACCCTGAACTTCGCCCCACATCACCCCGCCTACGCAGATGACAAAAA‑TCACTGACACCGCCCACAGCAACAGGCGCCCCGGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CTTCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTAAAA‑CGG‑ATGCGACGAGTGGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCtt                                                                                                                                                                                                               <  1:291518/251‑1 (MQ=255)
                                        ttACCCTGACCTTCTCCCCACACCCCCCCGCCTACGCAGATGAAAAAAA‑TCACTGACATCGCCCACAGCAACAGGCGCCCCTGGCGG‑CGATCGCCCATCCCTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATAAAGAAGA‑TCGCCAG‑CATCTGCACGAATTGGGTCAGTGCGGTTGGT‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCCCCGCCGTTAGTA‑CCGAGCATCTTGATCGCtt                                                                                                                                                                                                               <  1:260258/250‑1 (MQ=255)
                                          ccccTGACCTTCTGCCCACACCCCCCCGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCGGGCGG‑CGATCGCCCATCACTCCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGTTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCt                                                                                                                                                                                                            <  1:287740/250‑1 (MQ=255)
                                          ccccTGAACTTCTCCCCACATCCCCACGCCAACGCAGATGAAAAAAA‑CCACTGACACCCCCCACAGCAACATGCGCCCCGGGCGG‑CGACCCCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATAAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTGGGTCAGTGGGGTTGGG‑TTTCCAAA‑CGG‑ATGGGACGGGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCt                                                                                                                                                                                                            <  1:383316/250‑1 (MQ=255)
                                           cccTGCACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCCAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑CTGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTg                                                                                                                                                                                                           >  1:328778/1‑251 (MQ=255)
                                             cTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTg                                                                                                                                                                                                           >  1:499917/1‑249 (MQ=255)
                                                aaCTTTTGCCCACATCCCCACCCCTACCCAGATGACAAAAA‑CCCCTGACCCCCACCACCGCACCATGC‑CCCCCGGCGGGCGGCCGCCCATCCCTTCACCAAAGGAAAAGACCAGCCCCTTTGGGATCAAGAAAA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGT‑TTTTCAAC‑CGG‑ATGCGACGAGTTGGCATAAAAAAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAg                                                                                                                                                                                                      <  1:112796/251‑1 (MQ=255)
                                                 aCTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGc                                                                                                                                                                                                     >  2:141404/1‑251 (MQ=255)
                                                     cTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTAc                                                                                                                                                                                                  >  1:20658/1‑250 (MQ=255)
                                                        cccccaTCACCACGCCTACGCAGAGTACAAAAA‑TCACTAACATCGCCCACAGCAACATGCCCCCCTGGCGG‑CGATCCCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTGGGGATCAAGAAGA‑TCGCTGG‑CTTCTGCACGAAGTGGGTCAGTGCGGTTGGT‑TTTCCAAA‑CGG‑ATGCGACGAGTTGGCATAAAAGAACCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGc                                                                                                                                                                                              <  1:504565/247‑1 (MQ=255)
                                                         ccactcccccacGCCTACGCAGA‑GACAAAAAATCACTGACACGGCCCAAAGAAACAGGCGCCCGGGGCGG‑GGCTCGCCCACCACCTAACCAAAGGAAAAGCCCGGGGCGGTTGGGATAAGGAAGA‑TCTCCTG‑CATTGGCGAGAATTGGGTAAGTGGGGTGGGG‑TTTTAAAA‑GGG‑ATGCGACGGGTTGGCATAAAAGAAGCATACGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTAAAggac                                                                                                                                                                                             <  2:354573/243‑3 (MQ=255)
                                                             aTCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCAt                                                                                                                                                                                         >  2:56559/1‑251 (MQ=255)
                                                                      ccTACACAGATGAAAAAAA‑TCACGGACCTCCCCCACAGCAACAGGCCCCCCTGGCGG‑CGCCCGCCCATCACTCCCCAAAAGGAAAAGCACAGCCCCTTTGGAATCAAGAAGA‑TCGCCAG‑CTTCGGCACGAATTTGGTCAGTGCGGTTGGT‑TTTTCAAA‑CGG‑ATGCGACGATTTGGAATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAg                                                                                                                                                                                <  1:353046/251‑1 (MQ=255)
                                                                       cTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGc                                                                                                                                                                               >  1:234341/1‑251 (MQ=255)
                                                                          cGCAGATGAAAAAAA‑TCCCTGACACCGCCCACAGCAACATGCCCCCCTGGCGG‑CGACCGCCCATCACTTCACCAAAGGAAAAGCACAGCGCCGTGGGGATCAAGAAGA‑TCGCAAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTAAAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGt                                                                                                                                                                            <  1:193665/251‑1 (MQ=255)
                                                                           tcAGATGACAAAAA‑TCACTGACATCGCCTACAGCAACATGCGCCCCTGGAGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACGGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTACGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCAACTTGATCGCTTCCTGAGAACCTACAGGCCCCATGGGAAACAGCTGtt                                                                                                                                                                           >  2:76312/2‑251 (MQ=255)
                                                                             aGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTgc                                                                                                                                                                         >  1:257109/1‑251 (MQ=255)
                                                                               ataaCAAAAA‑TCCCTGACCCCGCCCCCCGCAACAGGCCCCCCGGCCGG‑CGATCCCCCATCACTTCACCAAAGTCAAAGCACAGCGCCGTTGGGATCAAAAAGA‑CCGCCAG‑CATCTGCACGAAGTTGGTCAGGGGGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTGGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTgcgc                                                                                                                                                                       <  1:181391/248‑1 (MQ=255)
                                                                                      aaa‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGCGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAAGCCCCATGGGTAACAGCTGTTGCGCTCCTTc                                                                                                                                                                 >  2:97675/1‑250 (MQ=255)
                                                                                            aCTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGt                                                                                                                                                           >  1:455334/1‑251 (MQ=255)
                                                                                                  accGCCCACACCAACATGCCCCCCGGGCGG‑CGACCGCCCATCACTCCCCCAAAGGCAAAGCACAGCGCGGTTGGGAACAAGAAGA‑TCGCCAG‑CATCTGCACGAGGTGGGTAAGTG‑GGTTGGGGTTTTAAAA‑CGG‑ATGCGACGATTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTcac                                                                                                                                                     <  1:378075/249‑1 (MQ=255)
                                                                                                    cGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTcac                                                                                                                                                     >  2:202106/1‑249 (MQ=255)
                                                                                                    cGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTcac                                                                                                                                                     >  1:346680/1‑249 (MQ=255)
                                                                                                       cccaaGCAACATGCGCCCCGGGCGG‑CGATCCCCCATCACTCCACAAAAGGCAAAGCACAGCGCCGTGGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGGGCGGTTGGG‑TTTCAAAA‑CGG‑ATGCGACGAGTTGGCATAAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTCCAACGGTATTCACAGCCt                                                                                                                                                <  1:241897/247‑1 (MQ=255)
                                                                                                          aaGCAACACGCGCCCCGGGCGG‑CTGCCGCCCATCTCATCACCAAGGGAAAAGACCACCCCCGTGGGGACCAAGAAGA‑TCACCAG‑CATCGGCAAGAAGTTTGTCAGGGCG‑TTGGG‑TTTGCAAAGCTG‑ATGCGCAGTGTTGGCTTAAAAGCCACCCCCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACTGTCCTCACAGCATGAt                                                                                                                                             <  2:178240/250‑1 (MQ=255)
                                                                                                           aGCAACAGGCGCCCCGGGCGG‑CGCCCGCCCATCACTACACAAAAGGCAAAGCACAGCGCCGTTGGGATAAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTGGGTCAGTGCGGTGGGG‑TTTTAAAA‑CGG‑ATGCGACGAGTTGGCATAAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATa                                                                                                                                            <  1:497418/251‑1 (MQ=255)
                                                                                                                 aTGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGc                                                                                                                                        >  1:113195/1‑249 (MQ=255)
                                                                                                                   gcgcCCCTGGCGG‑CGACCGCCCATCACTTCCCCAAAGGCAAACCACAGCGCCTTGGGGAT‑AAAAAGAACCGCAAG‑CATCTGCACGAGGTTGGTAGGGGCGGTTGGT‑TTTTAAAA‑GGG‑TTGCGACGAGTTGGCATTAAGGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGaa                                                                                                                                    <  1:477758/251‑1 (MQ=255)
                                                                                                                      ccccTGGCGG‑CGATCGCCAATCACTCCCCCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGt                                                                                                                                 <  1:267470/251‑1 (MQ=255)
                                                                                                                           ggcgg‑cGACCGCCCACCCCTCCACAAAAGGCAAAGCACGGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATAAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGc                                                                                                                            <  1:141404/251‑1 (MQ=255)
                                                                                                                             ggggctaTCGCCCATCCCTTACCCAAGGGAAAAGCAGGGCGCGGTGGGGATCAAGAAGA‑TCGCAAG‑CTTCTCCACAAAGTGTGTAGGTGCGTTTGGT‑TTTCCAAA‑CGG‑AGGCGACGGTTTGGCATAAAAGAGGCAACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCa                                                                                                                           <  1:56700/245‑1 (MQ=255)
                                                                                                                                  gATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAgcgc                                                                                                                       >  1:207543/1‑250 (MQ=255)
                                                                                                                                        ccATCATTTCACAAAAGCCAAAGCACAGCCCCGTTGGGACAAAGAAGA‑TCGCTCT‑CATCGGCACGAATTTGTTCAGTGCGGTGGGT‑TTTTAAAA‑GGG‑AGGCGACGAGTTGTCATTAAAGAGCCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCttgt                                                                                                                <  1:335384/251‑1 (MQ=255)
                                                                                                                                         cATCACTTCACAAAAGGCAAAGCACAGCGCCGTTGGAATAAAGAAGA‑TCGCAAG‑CATCGGCACGAAGTGGGTAAGGGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGATTGGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCttgtt                                                                                                               <  1:70493/251‑1 (MQ=255)
                                                                                                                                           cctctccccaaagggcaaagacgaGCGCCGTGGGGA‑CAAAAAGA‑CCGCCAGGCTTCGGCAGGAAGTTGGTCGGGGGGGTTGGT‑TTTCCAAA‑CGG‑GTGCGGCGGTTTGGAATAAAAGAACCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGa                                                                                                             <  2:346680/228‑1 (MQ=255)
                                                                                                                                            cACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGa                                                                                                             >  1:139274/1‑250 (MQ=255)
                                                                                                                                               cccacaaaaGGAAAAGCACAGCGCCTTGGGGATAAAGAAGA‑TCGCCAG‑CATTTGCGCGAAGTGGGGCAGGGCGGTTTGT‑TTTTCAAC‑CGG‑GTGCGACGAGTTGGCATAAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATaa                                                                                                         <  2:499917/245‑1 (MQ=255)
                                                                                                                                                 cACCAAAGACAAAGCACAGCGCCGTGGGGATAAAAAAGA‑GCGCCAG‑CATCTGCAAGAAGTGGGTCAGTGCGGTTGGT‑TTTTCAAA‑CGG‑ATGCGACGAGTGGGCATAAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATaaaa                                                                                                       <  1:420038/251‑1 (MQ=255)
                                                                                                                                                  cccAAAGGCAAAGCACGGCGCCGTTGGGACCAAGAAGT‑TCGCCGG‑CATCTGCAGGAAGTTGGTCGGGGCGGTTGGT‑TTTTCAAC‑GGG‑AGGCGACGAGTGGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATaaaaa                                                                                                      <  1:56559/250‑1 (MQ=255)
                                                                                                                                                      aaGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGACGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAAc                                                                                                    >  1:375718/1‑249 (MQ=255)
                                                                                                                                                         gCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGc                                                                                               >  1:185061/1‑251 (MQ=255)
                                                                                                                                                             ggCACAGCGCCGGTGGGTACAAGAAGA‑TCGCGAG‑CTTCTGCACGAGGTGGGTCGGTGGGGTTGGG‑TTAT‑AAA‑CGGGATGCGACGAGGTGACATAAAAGAGCCCCCCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATc                                                                                           <  2:209056/250‑1 (MQ=255)
                                                                                                                                                                 caGCGCCGTTGAGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAATAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACCGAGCAATcaac                                                                                        >  2:187370/1‑250 (MQ=255)
                                                                                                                                                                             gATCAAGAGGA‑TCGCGAG‑CATCTGCACGAAGTTGGTCAGTGAGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGATGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTCTTCACAGCCTGATCAGGCAGAAAGTTTTGCAGCGCACCTTGTTGACTAAAAAACAGTGCAATCAACAACGCCACAggg                                                                            >  2:273506/1‑250 (MQ=255)
                                                                                                                                                                              ataaaaaaga‑GCGCCAG‑CATCTGCGCAAGGTTGGTCAGTGCGGTTGGG‑TTTTAAAC‑CGG‑AGGCGCGGGGTTGGCATTAAAGAACCCTCCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGAc                                                                          <  1:389822/245‑1 (MQ=255)
                                                                                                                                                                                                      gaccgatgttggaaattggggttggt‑tTTTAAAA‑CGT‑ATGCGACGTGTTGACATAAAATACCCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGc                                                      <  1:148119/230‑1 (MQ=255)
                                                                                                                                                                                                           aaGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGc                                                 >  1:271412/1‑249 (MQ=255)
                                                                                                                                                                                                               tgttcggTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTGGGCATAAAAGAACCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATc                                            <  1:202106/244‑1 (MQ=255)
                                                                                                                                                                                                                ggTCAGTGCGGTTGGT‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGAATAAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCg                                           <  1:270446/250‑1 (MQ=255)
                                                                                                                                                                                                                             gggtttttCAAA‑CGG‑AGGCGAGGAGTTGACTTTAAAAAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTcc                              <  1:401421/247‑1 (MQ=255)
                                                                                                                                                                                                                               g‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCg                            >  2:213081/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aaa‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATcaa                                                                                         >  1:503889/1‑183 (MQ=255)
                                                                                                                                                                                                                                            g‑ATGCGACGAGTGGGAATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATg                <  1:388875/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aTGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGAATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGCGATCATGc               >  2:64315/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            cATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGCAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTaaa  >  2:421646/1‑250 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
ATTGATGCTGCTGTCCGTGCCCAGTGCCAGCAGATGAGGATTACCCTGAACTTCTGCCCACATCACCACGCCTACGCAGATGACAAAAA‑TCACTGACATCGCCCACAGCAACATGCGCCCCTGGCGG‑CGATCGCCCATCACTTCACCAAAGGCAAAGCACAGCGCCGTTGGGATCAAGAAGA‑TCGCCAG‑CATCTGCACGAAGTTGGTCAGTGCGGTTGGG‑TTTTCAAA‑CGG‑ATGCGACGAGTTGGCATTAAAGAAGCCACCGCCGTTAGTA‑CCGAGCATCTTGATCGCTTCCTGAGAAGCTACAGGCCCCATGGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAA  >  NC_000913/727776‑728270

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