Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 339,941 C→T T58I (ACT→ATT)  yahH → putative uncharacterized protein YahH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,9410CT89.7% 76.3 / NA 40T58I (ACT→ATT) yahHputative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (11/24);  minor base G (0/4);  total (12/28)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCC  >  NC_000913/339653‑340225
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                            
cGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATgcg                                                                                                                                                                                                                                                                                   >  1:674030/1‑300 (MQ=255)
    ggCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGAACTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGTATTTGTCACCACTGCCGGATGCGGCGt                                                                                                                                                                                                                                                                               >  2:302232/1‑300 (MQ=255)
                       ttGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                            >  1:699899/1‑300 (MQ=255)
                            cATCCAGCTTTTGCCTGCGTTGAGAGAGCTAGTTTTGGCTGGGTCGAATGTTACGGAGTCTGGCCTACATGTTCATGCTGGGTGCGGCGTGAACGTCTTATCCTGCCTACGAATGGTGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCATTAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAgc                                                                                                                                                                                                                                                       <  2:719326/300‑2 (MQ=14)
                                               tCGACAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGATCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGcc                                                                                                                                                                                                                                   <  1:242384/300‑1 (MQ=255)
                                                    agagctaTTCGTGGCTGGGGCGAACGTTACGGCGTCTGGCCTACATGTTCATGCNGGATGCGGCGTGAACGCCTTATCTGGCGTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATa                                                                                                                                                                                                                              <  2:625840/294‑1 (MQ=14)
                                                         cattttttattGGTGCGAACGTTACGGCGTCTGACTTGCATGTTCATGCCGGATTCGGCGTGAACGCCTTATCCGGCCTATGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTTACCATTGTCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGt                                                                                                                                                                                                                         <  2:453338/290‑1 (MQ=14)
                                                              ttggctagtgcgagcgttgctttgTCTGACCTACATGTTCAGGCCGGATGCGGCTTGAACACCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGAGAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGAGTGGTGCGAGAATTGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                    <  2:163703/277‑1 (MQ=11)
                                                              ttGGCTGGGGCGGACGTTATGGCGTCTGACCTACCTGTTTGTGCCGGATGCGGCATGAACGCCTTATTCGGCCTACGAGTGGCGCAAGAGGCTGTAGGCCTGATAAGCGTAGCGCATCATGCATTTGTTACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAACCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAGGCGTAgcgc                                                                                                                                                                                                                    <  2:177788/300‑1 (MQ=11)
                                                                         gAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGCTAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCAGCGTGAACGCCTTATCCGGCATACGCGTGGCGAGAGAATTGTAGGCCTGATAAGCGTAGAGCATCAggtagtc                                                                                                                                                                                                         >  2:638883/1‑295 (MQ=14)
                                                                              ttACGGCATCTGACCTACATGTTCATGCTGGAGGCGTCGTGAACGCCTTAGCCGGCCTACGAACGGCGCAAGAATCTGGAGGCCCGATAAGCGTAGCGCATCAGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGGATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcg                                                                                                                                                                                                    <  2:109706/300‑9 (MQ=11)
                                                                               tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCTCAAGAATCTGTAGGCGTGATAAGCTTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                   <  2:509667/300‑10 (MQ=11)
                                                                               tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCNTACGAATGGCGCAAGAATCTGTAGGCNTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                   <  2:569264/300‑10 (MQ=11)
                                                                               tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                   <  1:22707/300‑10 (MQ=255)
                                                                                    cGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCagcctggcgttggtc                                                                                                                                                                                              >  1:139529/1‑286 (MQ=11)
                                                                                                                                                                                      gacctaaatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCGGGCATTTGTCGCCATTGCCGGATGCGGTGTGAACGTCTTATCCGGCCTCCGAATGGCTAGAGAATTGTAAGCTTGATACGCGCCGCGCCCCAGGCAGTCTGGCGTTGGTCATagc                                                                                             >  2:19548/15‑298 (MQ=11)
                                                                                                                                                                                       acctgcgtgttcatGCCGGATGCGTCGTGAACACCTTATCCTGCCTGCGAATGCNGCAAGAATCTTTATGCCTGATAAGCGTTACGCATCGGGCATTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGAATGCGGCGTGAACGCCTTACCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                            <  2:635013/287‑1 (MQ=11)
                                                                                                                                                                                       acctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                            <  1:165262/287‑1 (MQ=11)
                                                                                                                                                                                       accgacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCAGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                            <  1:321387/287‑1 (MQ=11)
                                                                                                                                                                                        cctgcatgttcatGCCGGATGCGGCGTGAACGCCCTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCTTAGCACATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGGTGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                                           <  2:44894/288‑1 (MQ=11)
                                                                                                                                                                                                     gcctgATGCTGCGTGAACGCCTTATCCGGCCTACGAACGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCGGATAAGCGTGGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                              <  1:635894/296‑1 (MQ=14)
                                                                                                                                                                                                     gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATGAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                              <  1:798161/300‑1 (MQ=255)
                                                                                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                                                            <  1:264644/300‑1 (MQ=255)
                                                                                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGCGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                                                            <  1:266323/300‑1 (MQ=255)
                                                                                                                                                                                                          attcggcgTGAACGCCCTATCCGGCTTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                         <  2:434287/297‑1 (MQ=14)
                                                                                                                                                                                                          aTGCGGCGTGAACGCCTTATCTTGCCTACGAATTGCGCAAGCATCTGTAGGCCTAATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTACCGGATGCGGCGTGAACGCCTTATCCGGCCTACGATTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                         <  2:465861/300‑1 (MQ=14)
                                                                                                                                                                                                          aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAATAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                         <  1:560656/300‑1 (MQ=255)
                                                                                                                                                                                                              ggcgTGAACACCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCAACCTACAAATGGCGCAAGAATCTGTAGGCCTAATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGACATACGAATGGCGCGAGGGTTGTAGGCCTGATAAGCGCAGCGCACCAGACAGTCTGACGTCGGTCATAACCACCTCACCCTCTACAGCGGatat                                                                     >  2:303982/1‑297 (MQ=11)
                                                                                                                                                                                                                 gTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaac                                                                  <  1:453483/300‑1 (MQ=255)
                                                                                                                                                                                                                      cGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACAt                                                             >  1:151564/1‑300 (MQ=255)
                                                                                                                                                                                                                           gattcGGCCTACGAACGGCGCAAGAAGCTGTAGGCCTGATAAGCGTAGCGAAGCAGGCATTTGGCAACAGTGCCGGAGGCGGCGGGAACGCCTTATCCGGCCTACGAATNGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCTGACGAGTGGTGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAag                                                        <  1:600291/296‑1 (MQ=17)
                                                                                                                                                                                                                                 gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTCGGCCTGCTAAGCGTAGCGCATCAGGCCGTCTGGCGTTGGTCATAACCCCATCACCCCCTGTAGCGGACATAACAACATGAAAGAGCTtg                                                  >  2:237954/1‑300 (MQ=17)
                                                                                                                                                                                                                                            gcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGANTTGTCGGCCTGATACGCTTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGAACATAACAACATGAAAGAGCGTGTGGTCGTTGcc                                       >  2:727541/1‑300 (MQ=17)
                                                                                                                                                                                                                                             cacAAGNATCTGTATGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGNCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCa                                      <  2:695636/298‑1 (MQ=255)
                                                                                                                                                                                                                                              gcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCAt                                     >  1:734371/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     tctgtggacttgatgaGCGTGGCGCATCCTGCATTTGGCACTCGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGc                              <  2:358214/290‑1 (MQ=18)
                                                                                                                                                                                                                                                        gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAAc                           <  2:266863/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gCGTAGCGCATCAGGCATTTGTAACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCCGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGAGTGGCGCGAGAGTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATCTGGCGTTGGTCACCACCCAATCCACCTCTGTAGAGGAAACAACAACATGAAAGAGCCCGTGGTCGTTGACATTGGTGGCGACAGCATCATCAAAg              >  2:11973/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                       gTAGCGCATCAGGCGTTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGNCCTACGAGTGGCGCGAGAATTGTAGGNCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGAt            <  2:573573/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCNGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCNTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgccagcc  <  2:659969/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                            
CGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCC  >  NC_000913/339653‑340225

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