Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,176 A→G 18.5% intergenic (+295/‑283) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,1760AG18.5% 123.1 / 2.8 64intergenic (+295/‑283)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (27/23);  new base G (9/5);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGG  >  NC_000913/374950‑375411
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
cTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGGCGCGTTAGCGTCGCATGAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGc                                                                                                                                                                                                                                   <  1:1277498/239‑1 (MQ=255)
          cgaagggtGGCGGGAGAATTTGTAGGGTTGATAAGGAGCGTTAGGGTCGCATCAGGGATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                                                                                                                                          <  1:368389/231‑1 (MQ=11)
          cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTTTCGGCCTGATaaac                                                                                                                                                                                                                          >  1:651517/1‑236 (MQ=255)
          cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCCTATCCGGCCTACCGGGGGGGCGGGAATTTTTAGGCCTGTTAAGAc                                                                                                                                                                                                                         >  1:710885/1‑239 (MQ=255)
             cGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTTTAGGCCTGAtcaaacc                                                                                                                                                                                                                        >  1:190341/1‑231 (MQ=255)
                          aaTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATCGCGCGGGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATTAGGCATCTGCGACCGCCTGCCGGATGCGGGGCGACCGCCTTCTGCGCGCTCCGCCTCGCGCGGGAATTTTTAGGCCAGATagatagcgccggcgcgggagg                                                                                                                                                                                                         >  2:710188/1‑219 (MQ=11)
                                                  ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGGGGGAATTTGTAGGCCTGAGAAGACGAGTTAGCGGCGCATCAGGCATTTGCGCACGAATGCCGGATGGCGGGTGAACGCCTCATCCGGCCTACGGTGAGCGCGGGAGTTTGTAGGCCGGTAAAAAAGCGTTACAGGCGCATAAgcgaatgcgcaagccatggcgga                                                                                                                                                                                 >  2:1269098/1‑217 (MQ=2)
                                                      gttCGCCTGAGGGATCTGCGCGCGCCTGCCGGAGGCGGCGTGAACGCCTTATCCGGCTTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATgcg                                                                                                                                                                            <  2:468862/237‑1 (MQ=255)
                                                      gtgcccacgagggccttggggaggACGGCCGGGTGGGGGCTGACCGCCGTATCCGGCCTTCGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGGGTGAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATgcg                                                                                                                                                                            <  2:890410/216‑1 (MQ=11)
                                                           cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                                >  1:267682/1‑199 (MQ=38)
                                                           cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                                <  2:267682/199‑1 (MQ=38)
                                                           cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCATCAGGCATCTGCGCA‑CGACTGCCGGATgcggctgg                                                                                                                                                                       >  1:265211/1‑236 (MQ=16)
                                                                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCCGGCACCTGCGCACGACTGCCGGATGCGGCCGGAACGCCTTATCCGCCCCACGGGTGGCGCGAGCATTTTTTGGCCTGTGCAAATCTGTGACCGTCGCACCACGCACCTTCGCA‑CACCTGCCgcgatagggggcaccgcccttcc                                                                                                                                                            >  2:676369/1‑217 (MQ=2)
                                                                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACACGGTAGCGTGGCATAAGGCATCTGGGCA‑CGACGCCCGGGAGCGGGCGAAACGCCTTaa                                                                                                                                                             >  2:956119/1‑237 (MQ=11)
                                                                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTTCGTA‑CGACGTCCGGATGCGGGGTAAACGCCTTAt                                                                                                                                                             >  1:1162204/1‑238 (MQ=255)
                                                                             gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                                      <  1:1228931/238‑1 (MQ=38)
                                                                               cTGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACCCCTTATCCGGCCTACGGGGTGCGCGGGAATTTTTAAGCCCGATAAGAACCGGTTGCGGCCCATCAGGCATAT‑TGCACCCACCTCCGGATGCGGCAGAAACGCCTTaacgcgtaaccg                                                                                                                                                   >  2:1252399/1‑228 (MQ=11)
                                                                                 gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                 <  1:894817/239‑1 (MQ=39)
                                                                                         cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTa                                                                                                                                                                                                                                      <  1:180996/147‑1 (MQ=255)
                                                                                         cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGAGAATTTGTa                                                                                                                                                                                                                                      >  2:180996/1‑147 (MQ=255)
                                                                                            cgTGAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGGGAATTTGTAAGCCCGAGAAAACGCGTTAGCGTGCCATCCGGCATCTGCGCA‑CCACTGCGCGATGCGGCGTAAACGCCTTTACCGGCCCACGAGGGG‑CGCgggca                                                                                                                                     >  2:443768/1‑237 (MQ=11)
                                                                                                         tCCGGCCCACGGATGGGTGGGGAATTGGTAGGCCTGATAAGACGCTTTCGCGTCGCATCAGGCATCTGCGCACGCCTGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGa                                                                                                                        <  1:676816/239‑1 (MQ=14)
                                                                                                             gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGGAAAGCCCCTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGccgggtaca                                                                                                                    >  1:1256555/1‑232 (MQ=255)
                                                                                                             gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGCCGCGAGCATTTGTAGGCCTGATAAGACGCGTTCTCGTCGCCTCAGGCAACTGCGCA‑CCACTGCCGGGTGCGGCGTAAACGCCCTATCCGGCCCACCGATGG‑CGCGGGAATTTgcatcccggcgaag                                                                                                                    >  2:1128675/1‑226 (MQ=14)
                                                                                                               cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGGATGGCGCGGGAATTTGTACGCATGATAAGACGCGTTCGCAGCGCATCAGGCATCTGGCCA‑CGACTGCCCGGTGCGGCGTGAACGCCCTCTTCGGCCTCCGGGGAG‑CGGGGGAActtttcgggcgggagagag                                                                                                                  >  2:822895/1‑223 (MQ=14)
                                                                                                               cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTCGCGGCGCATCAGGCCTCTGGGCA‑CGCCTGCCGGATCCGGCGAAAACCCCCTACGCAGCCCACCGATGGCCGCGGGAAatttggggcagataaaaa                                                                                                                   >  2:668495/1‑224 (MQ=14)
                                                                                                                       ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑c                                                                                                                                                                                          >  1:1102430/1‑160 (MQ=40)
                                                                                                                       ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑c                                                                                                                                                                                          <  2:1102430/160‑1 (MQ=40)
                                                                                                                          gggggAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGGATCAGGCATCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                        <  1:483787/236‑1 (MQ=255)
                                                                                                                                        gCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGC‑CGCGGGAATTTGTAGGCCTGTTAAAACGCGTTAGCGGCGCCTCAGGCATCt                                                                                          >  1:846650/1‑238 (MQ=255)
                                                                                                                                             aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATTTGCGCA‑CGACTGCCGGATgc                                                                                                                                                                             >  1:732627/1‑151 (MQ=18)
                                                                                                                                             aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATgc                                                                                                                                                                             <  2:732627/151‑1 (MQ=33)
                                                                                                                                                  aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGGTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACCGCCGGATGCGGCGGAAACTCCTTATCCGGCCTACGGGTGGCCGCGGGAA‑TTGTAGGACCGATTAGACGCGTTCGCGTCCAATCAGGCACCTGCCAACGaacg                                                                               >  2:704861/1‑236 (MQ=11)
                                                                                                                                                  aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGCATTTGTAGGCCTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCACACGaat                                                                                >  1:1175219/1‑236 (MQ=255)
                                                                                                                                                               gCATCGGGCACCTGCGCACGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTTCGGGTGGCGCGAGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                                  <  1:253905/239‑1 (MQ=255)
                                                                                                                                                               ccATCGGGCATCTGCGCACGGCTTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                                  <  1:956119/238‑1 (MQ=255)
                                                                                                                                                                cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGCGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTAGCATCAGGCATCTGCGCA‑CGCCTGCCGGATGCGGCGAAAACGCCTTACCCGGCCTACGGCTGG‑CGCGGAGATTTGTAGGCATGATAAGAAGCGTTAGCGTCGGATACAGCATCTGccaccgcccgacggaggggataaa                                                                 >  2:253905/1‑216 (MQ=11)
                                                                                                                                                                         tcccagcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCGGCACACGACTGCCGGATGCGATAAACGCCTTGt                                                         >  1:932481/6‑238 (MQ=255)
                                                                                                                                                                          cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGTCCTGATAAGACGCGTTAGCGTTGCAACAGGCATCTGCACCCAACAGCCGGAGGCGACAAACGCCTTGTcc                                                       >  2:1255865/1‑239 (MQ=14)
                                                                                                                                                                          cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                        <  1:668495/238‑1 (MQ=255)
                                                                                                                                                                                aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATCG‑CGCCGGCATTTGTCAGCCTGATAAGACGCGGTAACGTCCACTCCGGCCACTCCCCCCCCCCGCCGGATGCCATAAAAGTATTTTCCgggcc                                                  >  2:620853/1‑235 (MQ=11)
                                                                                                                                                                                     tGCCGGTGGGGGAGTGAACGTCTTATCCGGCGTCGGGGTGGCGCGAGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGGGGCGGAAACGCCTTATCCGGCCTACGGATGG‑CGGGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                            <  2:190341/239‑1 (MQ=14)
                                                                                                                                                                                     tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGCCTTTTCCGGCCTACAtt                                             >  1:1018734/1‑238 (MQ=255)
                                                                                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTc                                           >  1:1155726/1‑238 (MQ=255)
                                                                                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGAA‑CGACTGCCGGATGCGGGGTAACAGCCTTATCCGGCCTACGGATCG‑GCCTGGGATTAGTAGGCCTGAGAAGATCCGGTACCGCCGCACACGGACTCTGCACACGACCGCCCGAGGCGATAAAAGTCTTGTCCGGcccaaaatac                                           >  2:304471/1‑230 (MQ=14)
                                                                                                                                                                                                 tgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCTACTGCGCA‑CGAGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                 <  1:676369/237‑1 (MQ=255)
                                                                                                                                                                                                 gcggAACCCCTTATCCGGCCTACGGGTGGCCGGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                 <  1:552341/235‑1 (MQ=255)
                                                                                                                                                                                                   tGAACGCCTTAACCGGGCCACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                               <  1:443768/238‑1 (MQ=255)
                                                                                                                                                                                                      cccctttttCCGGTCTGGGGGTGGGGGGGGAATTTGTAGGCCTGATAAGAGGGGGTAGGGGCGCATCAGGCATCTGCGAC‑CGACTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                            <  2:1224759/230‑1 (MQ=17)
                                                                                                                                                                                                      aCGGCTTAGCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                            <  1:822895/238‑1 (MQ=255)
                                                                                                                                                                                                      aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGTCATCTGCACACGCCTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTTGCCCGTAGGCAGTCAt                           >  1:1231537/1‑239 (MQ=255)
                                                                                                                                                                                                          cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                       <  1:1255865/239‑1 (MQ=255)
                                                                                                                                                                                                           tttcccGGGCTTCGGGTGGCGCGGGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                       <  1:1251375/234‑1 (MQ=255)
                                                                                                                                                                                                                  gCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGCCTGCCGGATGCGATAAACTTCTTGTCCTGCCTACATTTCGCCCGTAGGCGGTCATTAAATAGTTTTg               >  1:908779/1‑239 (MQ=255)
                                                                                                                                                                                                                    cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                        >  1:1128952/1‑196 (MQ=255)
                                                                                                                                                                                                                    cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                        <  2:1128952/196‑1 (MQ=255)
                                                                                                                                                                                                                      aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGG‑CGCGGGAATTTGTAGTCCTGATAAGACGCGTTATCGTCGCATCAGCCATCTGCACACGACGGCCGGATGCCATAAACCTCTTGTCCGGCCTACCTTTCCCCCGTAgtccacccttaaatatttttaatt            >  2:340280/1‑215 (MQ=14)
                                                                                                                                                                                                                        gggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGTTAAACGTCTTGTCCGGCCTCAATTTCGCCCGTAGGCAGTCCTTAAATGGTTCTGTTTAcc         >  1:724576/1‑238 (MQ=255)
                                                                                                                                                                                                                        gggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCAGCAGGCATCTGCGAA‑CGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg         <  1:340280/239‑1 (MQ=255)
                                                                                                                                                                                                                         ggTGGCGGGGGAATTTGGAGGCCTGAGAAGGCGGGTTAGGTGCCGATCAGGCATCTGCGCA‑CGTCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg         <  2:265211/238‑1 (MQ=18)
                                                                                                                                                                                                                            ggggggagaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCgg      <  1:1252399/233‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc     <  1:582353/239‑1 (MQ=255)
                                                                                                                                                                                                                             gggcgGGAATTTGTAGGCCCGATAAGGCGGGGTAGCGTCGCATCAGGGATCGGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc     <  1:1343364/236‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGCTAAGACGCGTTAGCGTCGCATCCGGCATCTGCACCCCACTGCCGGATGGGATAAACGGCTTTTCCGGCTTCAATTTCGCCCGTAGGCATGCATTAAATCGTTTTGTTTatcgg       >  2:403104/1‑232 (MQ=17)
                                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCCCACGACTGCCGGATGCGGTAAACGTCTTGTCCGGCCTCCATTTGGCCCGTAGGAGGCCATTAAATAGTTCTGATTAtgcg       >  1:760855/1‑236 (MQ=17)
                                                                                                                                                                                                                                ggagaATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg  <  1:56128/237‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
CTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCA‑CGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGG‑CGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGG  >  NC_000913/374950‑375411

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