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.

TGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCT  >  NC_000913/523143‑523637
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                
tgaGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGACTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGAGGTCCCCATGTCCAGGGTTTGGCAGGTGAAAGAATTGGCGCGCCCACCGGCTGGGCGGGCGCGGGGTGTGCGGGGGGGGTGAcct                                                                                                                                                                                                                                                            >  2:446145/1‑249 (MQ=255)
           aaggttATGTGAAAAAGATATAAATACAATAGCTGGTCTTGGGTGTCATGTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCGCTGTATATGAGCGGAAAACCAGCGGCGCGTCAGGTATATGTGACTCAGTAGGGCGGTCCCAGTGTCGAGGGTGCGGGAGTTGTAGGAGTTGTCGCCCCGACCGCCGTGGCGGGCTCGGTTTGTCCGGACGCTATGGCTGCGGGCAAccc                                                                                                                                                                                                                                               <  2:277884/246‑1 (MQ=255)
             gTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGGGGCGGGCTCGGGGTGTCCGGGCGGGATGACTTCGGGCAACCCggg                                                                                                                                                                                                                                            >  1:93520/1‑250 (MQ=255)
               ttATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCAGGGTTCGGCCGGGGTAAAAATTGGCGCGCCCACCGGCGTGGCGGGGTCGGGTGTGCCGGGGGGGGTGGCTTCGGGGAAccccgca                                                                                                                                                                                                                                           >  1:125259/1‑246 (MQ=255)
                    tgAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCCGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGGTAGAACTGGCGCGCCCCCCGGCGTGGCGGGCTCGGTGTGTCCCGGGGGGGTGGCTTCGGGCCACCCGGTaaaaccc                                                                                                                                                                                                                                     >  1:230452/1‑247 (MQ=255)
                       aaaaTATATAACTACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCGGGGGGGGGGACTTCGGGGAACCCGGGAAAACCgccg                                                                                                                                                                                                                                  >  1:417839/1‑249 (MQ=255)
                       aaaaTATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctg                                                                                                                                                                                                                                  <  2:6161/251‑1 (MQ=255)
                                aaaTACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGGGGGCTCGGGGTGTCCGGGGGGGGTGACTTCGGGCAACCCGGGAAATCCGCTGCCggggggg                                                                                                                                                                                                                        >  1:117769/1‑249 (MQ=255)
                                         aGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGTTCGGCAGGGGTAAGAAATGGCGCGCCCACCGGCGGGGCGGGCTCGGGGTGTCCGGGCGGGGTGGCTTCGGGGAAACCGGGAAAACCGCCGCTggggggggaggggcgg                                                                                                                                                                                                               >  1:549739/1‑240 (MQ=255)
                                             ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAGCCGCTGCGGGGGTGTAAGGTGCTGCCCg                                                                                                                                                                                                          <  2:69901/251‑1 (MQ=255)
                                                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGcc                                                                                                                                                                                                            <  2:85656/245‑1 (MQ=255)
                                                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGGTCGGGGGGGCCGGGCGGGGGGAATTCGGGCAACCCGGGAAATCCCCTGCTGGGGGGGGAGGGGCTGCCCgaacg                                                                                                                                                                                                      >  1:85656/1‑247 (MQ=255)
                                                     gtgtCATTTTATTTTTCTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGAc                                                                                                                                                                                                   <  2:149528/250‑1 (MQ=255)
                                                      tgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGGTCGGCCGGTGTAAGAATTGGGGCGCCCACCGGCGTGGCGGGCTCGGGGGGTCCCGGGGGGGGGACTTCCGGGAAACCGGGAAATCCCCTGCTGGGGGCGAAGGGGGTGCCCcgcggaacgg                                                                                                                                                                                                 >  1:329689/1‑245 (MQ=255)
                                                       gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCGGGGGGTGCGGGCGGGGTGACTTCGGGGAACCCGGGAAAACAGCTGCCGGGGGGGGAGGAGAGGCCCGGGGAGACggg                                                                                                                                                                                                >  1:119104/1‑250 (MQ=255)
                                                           ttttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGTGGGTCCGGGGGGGGTGAATTCGGGGCACCCGGGAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgagaagggactt                                                                                                                                                                                            >  1:105718/1‑243 (MQ=255)
                                                               attttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGGGGGGGGCTCGGGGTGTGCGGGCGGGGGGGCTTTGGGGCACCCGGGAAAACCGCTGCTGGGGGGGAAGGGGCTGCCCggggggggcggcctttggg                                                                                                                                                                                        >  1:501755/1‑234 (MQ=255)
                                                                   ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCCGGGTGTGCGGGGGGGGTGGCTTCCGGCCAACCCGGAAATCCGCTGCTGGGGGCGGAGGTGCTGGCCGGCggggccggccatgcgccgcc                                                                                                                                                                                    >  1:293976/1‑232 (MQ=255)
                                                                       ttgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGGCTTCGGGCAACCCGGGAAATCCGCTGCTGGGGGGGGAGGGGCTGCCCGGCGGGGCGGGCCCTGCGCTgccccgc                                                                                                                                                                               >  1:475671/1‑248 (MQ=255)
                                                                        tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGGGGGCTCGGTGTGGCCGGGGGGGGTTGCTTCGGGCAACCCGGCAAATTCCATCCTTGGGGGGCAGGTGCTGCCCGGCGCGAGGGACCCTGCGCTgccccgcc                                                                                                                                                                              >  1:529123/1‑251 (MQ=255)
                                                                         gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                                                                                                                                                                             <  2:321910/251‑1 (MQ=255)
                                                                           tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGGCCAGGGGTCGGCCGGTGTGAGAATTGGCGGGCCCCCCGGCGTGGCGTGCCCGGGGGGTCCGGGCGGGGTGGCTTCTGGCAACCCGGGAAATCCCCTGCCGGGGGCGGAGGGGGTGCCCGGCGGGGCGGGCCTTGCTCTgcccccccccc                                                                                                                                                                           >  2:545842/1‑244 (MQ=255)
                                                                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGTGTGTCCGGGCGGGGGGGCTTCGGGGAACCCGGGAAAACCGCCGCTTGGGGCGGAGGGGCTGCCCGGCGAGACGGGCCTTGCGCCTGCCCGCccgccgccggttctt                                                                                                                                                                >  1:250182/1‑245 (MQ=255)
                                                                                                   aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                                                                                                                                                   <  1:201140/251‑1 (MQ=255)
                                                                                                   aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGGGGCGGGCGCGGGGTGTCCGGGCGGGGTGACTTCGGGGAACCCCGGAAAACCGCCGCTGGGGGGGGAGGGGCTGCCCGGCGGGGCGGACCTTTGGCTGCCCGGCCCCCTGcccgttaatctctccccgcccc                                                                                                                                                   >  1:151627/1‑231 (MQ=255)
                                                                                                    aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGGCCCGGGTTCGGCCGGGGGAAGAATTGGGGCGGCCACCGGCGTGGGGGGCTCGGGGGGTTCGGGGGGGGGGGCTTCGGGGAACCCGGGAAATACGCTGGTGGGGGGGCGGGGGGGGCCCGGCGCGACGGGGCTTGCGCCGCCccgcccgcgcccgggccttctctcccccacccc                                                                                                                                                  >  2:12182/1‑226 (MQ=255)
                                                                                                               aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTAcc                                                                                                                                       <  1:162676/251‑1 (MQ=255)
                                                                                                                cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGCTCCGGGTGTTCGGGCGGGGTGACTTCGGGCAACCCGGGAAAACCGCTGCTGGGGGGGGAGGGGCTGCCCGGCGAGACGGACCTTTCGCTGCCCCGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAACAGCTACCg                                                                                                                                      >  1:217816/1‑251 (MQ=255)
                                                                                                                  ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCg                                                                                                                                      <  1:339058/249‑1 (MQ=255)
                                                                                                                      atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCCGGCGGGGTGACTTCTGGCAACCCCGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGGCCTTGCGGCGCCCCGCCCCCCGCCCTTCCTTCTCTCCCCCCCCTACagcagcccccggcccg                                                                                                                                 >  1:9906/1‑240 (MQ=255)
                                                                                                                               aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGGCCAGGGTTTCGCCGGTGTAAGAATTTGCGCGCCCCCCGGCGGCGCGGGCTCCGGGGGTTCTGGCGGGGTGGATCCGGGCAACCCGGTTAATCCCCTGCTGGGGGGGAAGGCGGTGCCCGGCGAGACCGACCTTCCGCTCCCCGCCCCGCTGCCGTTTCTTCTCTCCCGCACCCACCGCAGCGACCGGACGAAGaacccgcc                                                                                                                       >  1:42879/1‑244 (MQ=255)
                                                                                                                               aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGTTTGGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGGGTGCCCGGTGTGTCCGGGCGGGGTGGCTTTCGGGAAACCGGGCAATCCCCTGGCGGGGGGCGAGGGGGCGCCCGGCGGGGCGGACCCTTCGCCGCCCGGCCCCCGGGCGTTCCTTTTCTCCCGCACCCAAAGGCGCCTCCGGACGAAGACGcccgc                                                                                                                       >  1:45380/1‑248 (MQ=255)
                                                                                                                                  ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTACATTCTCTCCCGCACCTACAGCAGCTATCGGACGAAGACGCCTGCAcc                                                                                                                    <  2:42467/251‑1 (MQ=255)
                                                                                                                                                   agaTATGACTCAGTATGGCGGTCCCATTGGCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGcc                                                                                                   <  2:417839/251‑1 (MQ=255)
                                                                                                                                                    gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGGCCGGGCGGGATGACTTCGGGGAAACCCGGAAAACCGCCGGTGGGGGGGAAAGGGCTGCCCGGCGGGACGGACCTTTCGCTGCCCCGCCCGCCGCCGGTCCTTCTCTCCCCCACCCACAGCCGCCACCGGGCGGAAACGGCTGCACCGGTGGGGGGTTTTGGccc                                                                                                  >  1:562776/1‑251 (MQ=255)
                                                                                                                                                               gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGGTTTCGGCCCCGGCTGGaa                                                                                         >  1:217429/1‑249 (MQ=255)
                                                                                                                                                                      ggTCCCATTGGCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTc                                                                                <  2:349315/251‑1 (MQ=255)
                                                                                                                                                                               gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACCGCAGCTACCGGCCGACGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTGTTATATCCg                                                                       >  2:137010/1‑251 (MQ=255)
                                                                                                                                                                                     ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTAc                                                                  >  1:549230/1‑250 (MQ=255)
                                                                                                                                                                                     ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGGTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTAc                                                                  >  2:221522/1‑250 (MQ=255)
                                                                                                                                                                                      gTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGCGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTAc                                                                  >  1:235970/1‑249 (MQ=255)
                                                                                                                                                                                                  tAAGACTTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGGTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGAc                                                    >  1:166301/1‑251 (MQ=255)
                                                                                                                                                                                                   aaGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACg                                                   <  1:459314/251‑1 (MQ=255)
                                                                                                                                                                                                          ggCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGAt                                            >  1:522780/1‑251 (MQ=255)
                                                                                                                                                                                                          ggCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTGCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGAt                                            <  2:230452/251‑1 (MQ=255)
                                                                                                                                                                                                             gcgcCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACCGCTACGTGATGACGGACTGATACt                                         >  1:13349/1‑251 (MQ=255)
                                                                                                                                                                                                               gcCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTc                                        >  1:386805/1‑250 (MQ=255)
                                                                                                                                                                                                                 ccACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAAc                                     >  1:298841/1‑251 (MQ=255)
                                                                                                                                                                                                                      ggcgtggcgtGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGCGACGGACCTTGCGCTGCCCGGCCCGCTGCCGCTCATTCTCTCCCGCCCCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGAcc                                 >  1:575463/1‑249 (MQ=255)
                                                                                                                                                                                                                         gtggcgtgCTCGGTGTTTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGCAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcgg                             >  2:575481/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggggtgCGCGGTGTGTCCGTGCGGGATGATGTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGTCGGACCTTGCGCTGCCCGGCCCGCGGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCTGACGAAGACGCCTGCACCGGTGGGCGTTTGCGGCGCCGGCTGGGAAGCGCCCTCTGATAGCCGCTTACAGCTACGTGATGATGGACTGATACTCAACTATAAcggcg                           <  2:306154/248‑1 (MQ=255)
                                                                                                                                                                                                                           ggcgtgCTCGGTGTGGCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACCACAAcggcg                           >  2:464938/1‑251 (MQ=255)
                                                                                                                                                                                                                            gcgtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcggcgg                          <  2:93520/251‑1 (MQ=255)
                                                                                                                                                                                                                              gtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACggcgggc                        <  1:126899/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGGCTGATACTCAACGACAACggcggggcg                      >  2:417771/1‑248 (MQ=255)
                                                                                                                                                                                                                                     gtgtgtCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCAt                 <  2:125259/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gtCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACggcggggggggggattac             >  1:496707/1‑239 (MQ=255)
                                                                                                                                                                                                                                            ccGGGCGGGACGACTGCGGGCAACCCGGTAAATCCCCTGCTGGGGTCTAAGGTGCTGCCCGGCAAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCGACCGGAGGAAGACTCCTGCACCGGTGGGCGTTTTCGGCCCCGGGTGTAAGGTGCCTTCTGATTCCCTCTTACAGCTACGTGATGACGGACGGATGCTCAACGACaagtgcgggttgatcagtcac             <  2:533014/249‑20 (MQ=255)
                                                                                                                                                                                                                                                  gggATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcggcggctggggcagttccattgt        >  2:793/1‑230 (MQ=255)
                                                                                                                                                                                                                                                     aTGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCTCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATTCTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCgct  <  2:533099/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                
TGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCT  >  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: