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,9410CT97.1% 94.3 / ‑4.1 36T58I (ACT→ATT) yahHputative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (14/20);  minor base G (0/1);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCA  >  NC_000913/339649‑340226
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                             
cGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGa                                                                                                                                                                                                                                                                                        <  2:651750/300‑1 (MQ=255)
   gCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATgc                                                                                                                                                                                                                                                                                     <  1:483069/300‑1 (MQ=255)
        ggCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGt                                                                                                                                                                                                                                                                                >  1:373108/1‑300 (MQ=255)
                     cgtgCATTGTACATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                                                                                                                                   <  1:479102/297‑1 (MQ=255)
                      gggCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTNTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                                                                                                                                  <  1:592663/300‑1 (MQ=255)
                             gtgCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAANAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTNAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGACTGATAAGAGTAGCGCATCAGGGATTTGTAAACATTGCCGGATGTGGCGTGAACGACTTATTCGGCCTata                                                                                                                                                                                                                                                           >  2:682562/1‑298 (MQ=18)
                               gCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAg                                                                                                                                                                                                                                                         <  1:322652/300‑1 (MQ=255)
                                      ccTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGAGCAAGAATCTGTAGGCCTGACACGCGTAGCGCATCAGGCATTTGTCGCCATTGCCGGATACGGCGTAAATGCATTATCCGGCCTACGAATGGcgcg                                                                                                                                                                                                                                                  >  2:182980/1‑300 (MQ=18)
                                              cTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAttgt                                                                                                                                                                                                                                          <  1:668764/300‑4 (MQ=255)
                                                        aaaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATa                                                                                                                                                                                                                               <  1:6640/300‑1 (MQ=255)
                                                               ttCTTGGCTGGTGCGAACGTTACGGCGCCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGTGTAGCGCATCAGGNATTTGTCACCATTGCCGGATGCGGCGNGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATAAGGCATTTGTCACCATCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGAGAGAATTGTAGGCCTGATAAGTGTAg                                                                                                                                                                                                                        >  2:570312/1‑300 (MQ=14)
                                                                  ttGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGNAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                     <  1:617827/300‑1 (MQ=255)
                                                                  ttGGCTGGTGCGAACGTTACGGCGGCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAgcgc                                                                                                                                                                                                                     <  1:379918/300‑1 (MQ=255)
                                                                         gTGCGAACTTTACGGCTTCTGACTTACCTGTTCGTGCCGGAGGCTGCGTGGACGCCTTATCCTTCCTACGACTGGCGCAAGAATCTCTAGGCCTGATGAGCTTAGCACATCATGTATTTTTTGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCAGGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGACATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                                                                                                                                              <  2:186446/300‑1 (MQ=11)
                                                                         gTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                                                                                                                                              <  1:316443/300‑1 (MQ=255)
                                                                         gTGCGAACGTTACGGCGTCTGACCTACACGTTCATGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGTCTGATAAGCGTAGCGCATCAGGCATTTGTCACCACTGCCGGATGCGGCGTGAACGCCTTATCTGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCACTAGGCATTTGTCACCACTGCCGACTGCGGCGGGAACGCCTTATCCGGCCTACGAGCGGCGGGAGACTTGTAGGCCTGGTACACGTAGCGCATCAGGc                                                                                                                                                                                                              >  2:38405/1‑300 (MQ=14)
                                                                                          tCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGTCCTGATAAGCTTAGCGCATCAGGCATTTGTCACCANGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGTCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAgcgcctgtctcttatacacatctcgcag                                                                                                                                                                                             <  1:581213/300‑25 (MQ=11)
                                                                                             gACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGTGAGGACTGTAGGCCTGACGAGCGTAATGCATCAGGCagtctggcgtcggtcataac                                                                                                                                                                                          >  2:163048/1‑281 (MQ=11)
                                                                                                                                                                                    gcgtttgacctacatgttcgtgctggATGTGGCGTGAATACCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGt                                                                                                    <  2:164910/276‑1 (MQ=11)
                                                                                                                                                                                          gacctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTGGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAc                                                                                              <  1:556639/286‑1 (MQ=11)
                                                                                                                                                                                           acctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                             >  1:514043/14‑300 (MQ=11)
                                                                                                                                                                                                              aTGCTGCGTGAACGCCTTATCCGGCCTACGAATGGCACAAGAATCTGTAGGCCTNATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATANGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                          <  2:787266/300‑1 (MQ=14)
                                                                                                                                                                                                              aTGCGGCGTGAACGCCTTGTCCGGCCTACGAAGGGCACAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                          <  2:387057/300‑1 (MQ=14)
                                                                                                                                                                                                              aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                          <  1:8012/300‑1 (MQ=255)
                                                                                                                                                                                                                  ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcctgtctcttatacacatctcgcag                                                                      <  1:771577/300‑24 (MQ=11)
                                                                                                                                                                                                                     gTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaac                                                                   <  1:538928/300‑1 (MQ=255)
                                                                                                                                                                                                                     gTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaaa                                                                   >  1:682943/1‑299 (MQ=255)
                                                                                                                                                                                                                               tATCCTGCCTACGAATGGCGAAAGAATCTGTAGGCCTGATAAGCGTAACGCATCAGGCATTGTTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAag                                                         <  1:209049/300‑1 (MQ=17)
                                                                                                                                                                                                                                      ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGACATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgt                                                  >  1:733297/1‑300 (MQ=255)
                                                                                                                                                                                                                                         aCGAATGGCGTAAGAGTCTGTAGGCCTGATAAGCGTAGCACAGCAGGCATTTGTCATCATTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGACTGGATCAAAAATCTTTGGGCTTGATAAGCTTAGCGCATCAGGCATTTGTCCCCATTGCCGGATTCGGCGTGAACGCCTCATCCGTCTTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGt                                               <  2:435377/300‑1 (MQ=255)
                                                                                                                                                                                                                                                gcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAACCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGcc                                        >  2:511569/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          cTGTAGGCCTGATAAGTGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCATTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTCCGAGTGGCGCGAGGATTGTCGGCCTGATGAGCTTAGCGCATCAGGCAGTCTGGCGTTGTTCATAAACCCATCACCCCCCGTAGCGGACATAAACACATGGAAGAGCGTGTGGTTGTTGCCATtggtggtg                              >  2:544990/1‑298 (MQ=14)
                                                                                                                                                                                                                                                               ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGc                         >  1:660515/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 ccTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAGCGCCTTATCCGGCATACGAGTGGCGCGAGAATTGTCGGCCTGATAAGCGTAGCGCATCAGGCACTCTGGCGTTGGTCATAACCCAATCACCCTCTGTAGCGGACATAGAAAAACGAAAGAGCTTGTGGTCGTTGCCATTGGTGACAACAGCAt                       >  2:337197/1‑300 (MQ=17)
                                                                                                                                                                                                                                                                         gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGGATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGAGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAg               >  2:769886/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                 cATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGCGGAGCGAGAATTGTAGGCCTGCTAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAACCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCCACAGCATTATCCAAGATacagcc       >  2:53431/1‑295 (MQ=18)
                                                                                                                                                                                                                                                                                      agcagttgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCATGCATTTGTCACCATTGCCGGATGCGGCATGAACGCCTTATCCGGCCTACAAGTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCCGGCGTTGGTCATAACTCCATCACCCTCTGTAGCGGACAGAACAACATGAAAGAGTTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgccagcca  <  2:158909/291‑1 (MQ=18)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                             
CGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCA  >  NC_000913/339649‑340226

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: