Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3870AG37.3% 71.0 / 15.6 59M43V (ATG→GTG) rhsDprotein RhsD
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base G (22/0);  total (40/19)
Fisher's exact test for biased strand distribution p-value = 2.11e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCT  >  NC_000913/523143‑523637
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                
tgaGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGACTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGAGGTCCCCATGTCCAGGGTTTGGCAGGTGAAAGAATTGGCGCGCCCACCGGCTGGGCGGGCGCGGGGTGT‑GCGGGGGGGGTGAcct                                                                                                                                                                                                                                                            >  2:446145/1‑249 (MQ=255)
           aaggttATGTGAAAAAGATATAAATACAATAGCTGGTCTTGGGTGTCATGTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCGCTGTATATGAGCGGAAAACCAGCGGCGCGTCAGGTATATGTGACTCAGTAGGGCGGTCCCAGTGTCGAGGGTGCGGGAGTTGTAGGAGTTGTCGCCCCGACCGCCGTGGCGGGCTCGGTTTGT‑CCGGACGCTATGGCTG‑C‑GGG‑CAAccc                                                                                                                                                                                                                                               <  2:277884/246‑1 (MQ=255)
             gTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGGGGCGGGCTCGGGGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCggg                                                                                                                                                                                                                                            >  1:93520/1‑250 (MQ=255)
               ttATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCAGGGTTCGGCCGGGGTAAAAATTGGCGCGCCCACCGGCGTGGCGGGGTCGG‑GTGTGCCGGGGGGGGTGGCTT‑C‑GGG‑GAAccccgca                                                                                                                                                                                                                                           >  1:125259/1‑246 (MQ=255)
                    tgAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCCGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGGTAGAACTGGCGCGCCCCCCGGCGTGGCGGGCTCGGTGTGT‑CCCGGGGGGGTGGCTT‑C‑GGG‑CCACCCGGTaaaaccc                                                                                                                                                                                                                                     >  1:230452/1‑247 (MQ=255)
                       aaaaTATATAACTACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGT‑CCGGGGGGGGGGACTT‑C‑GGG‑GAACCCGGGAAAACCgccg                                                                                                                                                                                                                                  >  1:417839/1‑249 (MQ=255)
                       aaaaTATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCgctg                                                                                                                                                                                                                                  <  2:6161/251‑1 (MQ=255)
                                aaaTACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGGGGGCTCGGGGTGT‑CCGGGGGGGGTGACTT‑C‑GGG‑CAACCCGGGAAATCCGCTG‑CCggggggg                                                                                                                                                                                                                        >  1:117769/1‑249 (MQ=255)
                                         aGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGTTCGGCAGGGGTAAGAAATGGCGCGCCCACCGGCGGGGCGGGCTCGGGGTGT‑CCGGGCGGGGTGGCTT‑C‑GGG‑GAAACCGGGAAAACCGCCG‑CTggggggggaggggcgg                                                                                                                                                                                                               >  1:549739/1‑240 (MQ=255)
                                             ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAAGCCGCTG‑C‑GGGGGTGTAAGGTGCTGCCCg                                                                                                                                                                                                          <  2:69901/251‑1 (MQ=255)
                                                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGcc                                                                                                                                                                                                            <  2:85656/245‑1 (MQ=255)
                                                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGGTCGGGGGGG‑CCGGGCGGGGGGAATT‑C‑GGG‑CAACCCGGGAAATCCCCTG‑CTGGGGGGG‑GAGGGGCTGCCCgaacg                                                                                                                                                                                                      >  1:85656/1‑247 (MQ=255)
                                                     gtgtCATTTTATTTTTCTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGAc                                                                                                                                                                                                   <  2:149528/250‑1 (MQ=255)
                                                      tgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGGTCGGCCGGTGTAAGAATTGGGGCGCCCACCGGCGTGGCGGGCTCGGGGGGT‑CCCGGGGGGGGGACTT‑C‑CGG‑GAAACCGGGAAATCCCCTG‑CTGGGGGCG‑AAGGGGGTGCCCcgcggaacgg                                                                                                                                                                                                 >  1:329689/1‑245 (MQ=255)
                                                       gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCGGGGGGT‑GCGGGCGGGGTGACTT‑C‑GGG‑GAACCCGGGAAAACAGCTG‑CCGGGGGGG‑GAGGAGAGGCCCGGGGAGACggg                                                                                                                                                                                                >  1:119104/1‑250 (MQ=255)
                                                           ttttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGTGGGT‑CCGGGGGGGGTGAATT‑C‑GGG‑GCACCCGGGAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCgagaagggactt                                                                                                                                                                                            >  1:105718/1‑243 (MQ=255)
                                                               attttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGGGGGGGGCTCGGGGTGT‑GCGGGCGGGGGGGCTT‑T‑GGG‑GCACCCGGGAAAACCGCTG‑CTGGGGGGG‑AAGGGGCTGCCCggggggggcggcctttggg                                                                                                                                                                                        >  1:501755/1‑234 (MQ=255)
                                                                   ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCCGGGTGT‑GCGGGGGGGGTGGCTT‑CCGGC‑CAACCCGG‑AAATCCGCTG‑CTGGGGGCG‑GAGGTGCTGGCCGGCggggccggccatgcgccgcc                                                                                                                                                                                    >  1:293976/1‑232 (MQ=255)
                                                                       ttgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGT‑CCGGGCGGGGTGGCTT‑C‑GGG‑CAACCCGGGAAATCCGCTG‑CTGGGGGGG‑GAGGGGCTGCCCGGCGGGGCGGGCCCTGCGC‑Tgccccgc                                                                                                                                                                               >  1:475671/1‑248 (MQ=255)
                                                                        tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCGGTGTGG‑CCGGGGGGGGTTGCTT‑C‑GGG‑CAACCCGGCAAATTCCATC‑CTTGGGGGG‑CAGGTGCTGCCCGGCGCGAGGGACCCTGCGC‑Tgccccgcc                                                                                                                                                                              >  1:529123/1‑251 (MQ=255)
                                                                         gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑Tgcccggccc                                                                                                                                                                             <  2:321910/251‑1 (MQ=255)
                                                                           tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGGCCAGGGGTCGGCCGGTGTGAGAATTGGCGGGCCCCCCGGCGTGGCGTGCCCGGGGGGT‑CCGGGCGGGGTGGCTT‑C‑TGG‑CAACCCGGGAAATCCCCTG‑CCGGGGGCG‑GAGGGGGTGCCCGGCGGGGCGGGCCTTGCTC‑Tgcccccccccc                                                                                                                                                                           >  2:545842/1‑244 (MQ=255)
                                                                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGTGTGT‑CCGGGCGGGGGGGCTT‑C‑GGG‑GAACCCGGGAAAACCGCCG‑CTTGGGGCG‑GAGGGGCTGCCCGGCGAGACGGGCCTTGCGCCTGCCCG‑Cccgccgccggttctt                                                                                                                                                                >  1:250182/1‑245 (MQ=255)
                                                                                                   aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                                                                                                                                                   <  1:201140/251‑1 (MQ=255)
                                                                                                   aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGCGCGGGGTGT‑CCGGGCGGGGTGACTT‑C‑GGG‑GAACCCCGGAAAACCGCCG‑CTGGGGGGG‑GAGGGGCTGCCCGGCGGGGCGGACCTTTGGC‑TGCCCGGCCCCCTGcccgttaatctctccccgcccc                                                                                                                                                   >  1:151627/1‑231 (MQ=255)
                                                                                                    aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGGCCCGGGTTCGGCCGGGGGAAGAATTGGGGCGGCCACCGGCGTGGGGGGCTCGGGGGGT‑TCGGGGGGGGGGGCTT‑C‑GGG‑GAACCCGGGAAATACGCTG‑GTGGGGGGG‑CGGGGGGGGCCCGGCGCGACGGGGCTTGCGC‑CGCCccgcccgcgcccgggccttctctcccccacccc                                                                                                                                                  >  2:12182/1‑226 (MQ=255)
                                                                                                               aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTAcc                                                                                                                                       <  1:162676/251‑1 (MQ=255)
                                                                                                                cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGCTCCGGGTGT‑TCGGGCGGGGTGACTT‑C‑GGG‑CAACCCGGGAAAACCGCTG‑CTGGGGGGG‑GAGGGGCTGCCCGGCGAGACGGACCTTTCGC‑TGCCCCGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAACAGCTACCg                                                                                                                                      >  1:217816/1‑251 (MQ=255)
                                                                                                                  ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCg                                                                                                                                      <  1:339058/249‑1 (MQ=255)
                                                                                                                      atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCCGGCGGGGTGACTT‑C‑TGG‑CAACCCCGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGGCCTTGCGG‑CGCCCCGCCCCCCGCCCTTCCTTCTCTCCCCCCCCTACagcagcccccggcccg                                                                                                                                 >  1:9906/1‑240 (MQ=255)
                                                                                                                               aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGGCCAGGGTTTCGCCGGTGTAAGAATTTGCGCGCCCCCCGGCGGCGCGGGCTCCGGGGGT‑TCTGGCGGGGTGGATC‑C‑GGG‑CAACCCGGTTAATCCCCTG‑CTGGGGGGG‑AAGGCGGTGCCCGGCGAGACCGACCTTCCGC‑TCCCCGCCCCGCTGCCGTTTCTTCTCTCCCGCACCCACCGCAGCGACCGGACGAAGaacccgcc                                                                                                                       >  1:42879/1‑244 (MQ=255)
                                                                                                                               aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGTTTGGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGGGTGCCCGGTGTGT‑CCGGGCGGGGTGGCTTTC‑GGG‑‑AAACCGGGCAATCCCCTGGCGGGGGGCG‑AGGGGGC‑GCCCGGCGGGGCGGACCCTTCGC‑CGCCCGGCCCCCGGGCGTTCCTTTTCTCCCGCACCCAAAGGCGCCTCCGGACGAAGACGcccgc                                                                                                                       >  1:45380/1‑248 (MQ=255)
                                                                                                                                  ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTACATTCTCTCCCGCACCTACAGCAGCTATCGGACGAAGACGCCTGCAcc                                                                                                                    <  2:42467/251‑1 (MQ=255)
                                                                                                                                                   agaTATGACTCAGTATGGCGGTCCCATTGGCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGcc                                                                                                   <  2:417839/251‑1 (MQ=255)
                                                                                                                                                    gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGG‑CCGGGCGGGATGACTT‑C‑GGGGAAACCCGG‑AAAACCGCCG‑GTGGGGGGG‑AAAGGGCTGCCCGGCGGGACGGACCTTTCGC‑TGCCCCGCCCGCCGCCGGTCCTTCTCTCCCCCACCCACAGCCGCCACCGGGCGGAAACGGCTGCACCGGTGGGGGGTTTTGGccc                                                                                                  >  1:562776/1‑251 (MQ=255)
                                                                                                                                                               gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGGTTTCGGCCCCGGCTGGaa                                                                                         >  1:217429/1‑249 (MQ=255)
                                                                                                                                                                      ggTCCCATTGGCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTc                                                                                <  2:349315/251‑1 (MQ=255)
                                                                                                                                                                               gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACCGCAGCTACCGGCCGACGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTGTTATATCCg                                                                       >  2:137010/1‑251 (MQ=255)
                                                                                                                                                                                     ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTAc                                                                  >  1:549230/1‑250 (MQ=255)
                                                                                                                                                                                     ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGGTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTAc                                                                  >  2:221522/1‑250 (MQ=255)
                                                                                                                                                                                      gTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGCGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTAc                                                                  >  1:235970/1‑249 (MQ=255)
                                                                                                                                                                                                  tAAGACTTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGGTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGAc                                                    >  1:166301/1‑251 (MQ=255)
                                                                                                                                                                                                   aaGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACg                                                   <  1:459314/251‑1 (MQ=255)
                                                                                                                                                                                                          ggCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGAt                                            >  1:522780/1‑251 (MQ=255)
                                                                                                                                                                                                          ggCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTGCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGAt                                            <  2:230452/251‑1 (MQ=255)
                                                                                                                                                                                                             gcgcCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACCGCTACGTGATGACGGACTGATACt                                         >  1:13349/1‑251 (MQ=255)
                                                                                                                                                                                                               gcCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTc                                        >  1:386805/1‑250 (MQ=255)
                                                                                                                                                                                                                 ccACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAAc                                     >  1:298841/1‑251 (MQ=255)
                                                                                                                                                                                                                      ggcgtggcgtGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGCGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGCTCATTCTCTCCCGCCCCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGAcc                                 >  1:575463/1‑249 (MQ=255)
                                                                                                                                                                                                                         gtggcgtgCTCGGTGTTT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGCAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcgg                             >  2:575481/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggggtgCGCGGTGTGT‑CCGTGCGGGATGATGT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGTCGGACCTTGCGC‑TGCCCGGCCCGCGGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCTGACGAAGACGCCTGCACCGGTGGGCGTTTGCGGCGCCGGCTGGGAAGCGCCCTCTGATAGCCGCTTACAGCTACGTGATGATGGACTGATACTCAACTATAAcggcg                           <  2:306154/248‑1 (MQ=255)
                                                                                                                                                                                                                           ggcgtgCTCGGTGTGG‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACCACAAcggcg                           >  2:464938/1‑251 (MQ=255)
                                                                                                                                                                                                                            gcgtgCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcggcgg                          <  2:93520/251‑1 (MQ=255)
                                                                                                                                                                                                                              gtgCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACggcgggc                        <  1:126899/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGGCTGATACTCAACGACAACggcggggcg                      >  2:417771/1‑248 (MQ=255)
                                                                                                                                                                                                                                     gtgtgt‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCAt                 <  2:125259/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gt‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACggcggggggggggattac             >  1:496707/1‑239 (MQ=255)
                                                                                                                                                                                                                                            ccGGGCGGGACGACTG‑C‑GGG‑CAACCCGGTAAATCCCCTG‑CTGGGGTCT‑AAGGTGCTGCCCGGCAAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCGACCGGAGGAAGACTCCTGCACCGGTGGGCGTTTTCGGCCCCGGGTGTAAGGTGCCTTCTGATTCCCTCTTACAGCTACGTGATGACGGACGGATGCTCAACGACaagtgcgggttgatcagtcac             <  2:533014/249‑20 (MQ=255)
                                                                                                                                                                                                                                                  gggATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcggcggctggggcagttccattgt        >  2:793/1‑230 (MQ=255)
                                                                                                                                                                                                                                                     aTGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCTCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATTCTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCgct  <  2:533099/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                
TGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTT‑C‑GGG‑CAACCCGGTAAATCCGCTG‑CTGGGGGCG‑AAGGTGCTGCCCGGCGAGACGGACCTTGCGC‑TGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCT  >  NC_000913/523143‑523637

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