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

CTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGG  >  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)
                                                      gttCGCCTGAGGGATCTGCGCGCGCCTGCCGGAGGCGGCGTGAACGCCTTATCCGGCTTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                                                                                                                                                                            <  2:468862/237‑1 (MQ=255)
                                                      gtgcccacgagggccttggggaggACGGCCGGGTGGGGGCTGACCGCCGTATCCGGCCTTCGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGGGTGAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                                                                                                                                                                            <  2:890410/216‑1 (MQ=11)
                                                           cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                                >  1:267682/1‑199 (MQ=38)
                                                           cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                                <  2:267682/199‑1 (MQ=38)
                                                           cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCATCAGGCATCTGCGCACGACTGCCGGATgcggctgg                                                                                                                                                                       >  1:265211/1‑236 (MQ=16)
                                                                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCCGGCACCTGCGCACGACTGCCGGATGCGGCCGGAACGCCTTATCCGCCCCACGGGTGGCGCGAGCATTTTTTGGCCTGTGCAAATCTGTGACCGTCGCACCACGCACCTTCGCACACCTGCCgcgatagggggcaccgcccttcc                                                                                                                                                            >  2:676369/1‑217 (MQ=2)
                                                                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACACGGTAGCGTGGCATAAGGCATCTGGGCACGACGCCCGGGAGCGGGCGAAACGCCTTaa                                                                                                                                                             >  2:956119/1‑237 (MQ=11)
                                                                      tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTTCGTACGACGTCCGGATGCGGGGTAAACGCCTTAt                                                                                                                                                             >  1:1162204/1‑238 (MQ=255)
                                                                             gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                                      <  1:1228931/238‑1 (MQ=38)
                                                                               cTGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACCCCTTATCCGGCCTACGGGGTGCGCGGGAATTTTTAAGCCCGATAAGAACCGGTTGCGGCCCATCAGGCATATTGCACCCACCTCCGGATGCGGCAGAAACGCCTTaacgcgtaaccg                                                                                                                                                   >  2:1252399/1‑228 (MQ=11)
                                                                                 gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                 <  1:894817/239‑1 (MQ=39)
                                                                                         cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTa                                                                                                                                                                                                                                      <  1:180996/147‑1 (MQ=255)
                                                                                         cggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGAGAATTTGTa                                                                                                                                                                                                                                      >  2:180996/1‑147 (MQ=255)
                                                                                            cgTGAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGGGAATTTGTAAGCCCGAGAAAACGCGTTAGCGTGCCATCCGGCATCTGCGCACCACTGCGCGATGCGGCGTAAACGCCTTTACCGGCCCACGAGGGGCGCgggca                                                                                                                                     >  2:443768/1‑237 (MQ=11)
                                                                                                         tCCGGCCCACGGATGGGTGGGGAATTGGTAGGCCTGATAAGACGCTTTCGCGTCGCATCAGGCATCTGCGCACGCCTGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGa                                                                                                                        <  1:676816/239‑1 (MQ=14)
                                                                                                             gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGAAAGCCCCTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGccgggtaca                                                                                                                    >  1:1256555/1‑232 (MQ=255)
                                                                                                             gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGCCGCGAGCATTTGTAGGCCTGATAAGACGCGTTCTCGTCGCCTCAGGCAACTGCGCACCACTGCCGGGTGCGGCGTAAACGCCCTATCCGGCCCACCGATGGCGCGGGAATTTgcatcccggcgaag                                                                                                                    >  2:1128675/1‑226 (MQ=14)
                                                                                                               cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGGATGGCGCGGGAATTTGTACGCATGATAAGACGCGTTCGCAGCGCATCAGGCATCTGGCCACGACTGCCCGGTGCGGCGTGAACGCCCTCTTCGGCCTCCGGGGAGCGGGGGAActtttcgggcgggagagag                                                                                                                  >  2:822895/1‑223 (MQ=14)
                                                                                                               cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTCGCGGCGCATCAGGCCTCTGGGCACGCCTGCCGGATCCGGCGAAAACCCCCTACGCAGCCCACCGATGGCCGCGGGAAatttggggcagataaaaa                                                                                                                   >  2:668495/1‑224 (MQ=14)
                                                                                                                       ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCAc                                                                                                                                                                                          >  1:1102430/1‑160 (MQ=40)
                                                                                                                       ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCAc                                                                                                                                                                                          <  2:1102430/160‑1 (MQ=40)
                                                                                                                          gggggAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGGATCAGGCATCGGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                        <  1:483787/236‑1 (MQ=255)
                                                                                                                                        gCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGCCGCGGGAATTTGTAGGCCTGTTAAAACGCGTTAGCGGCGCCTCAGGCATCt                                                                                          >  1:846650/1‑238 (MQ=255)
                                                                                                                                             aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATTTGCGCACGACTGCCGGATgc                                                                                                                                                                             >  1:732627/1‑151 (MQ=18)
                                                                                                                                             aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgc                                                                                                                                                                             <  2:732627/151‑1 (MQ=33)
                                                                                                                                                  aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGGTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACCGCCGGATGCGGCGGAAACTCCTTATCCGGCCTACGGGTGGCCGCGGGAATTGTAGGACCGATTAGACGCGTTCGCGTCCAATCAGGCACCTGCCAACGaacg                                                                               >  2:704861/1‑236 (MQ=11)
                                                                                                                                                  aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGGCATTTGTAGGCCTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCACACGaat                                                                                >  1:1175219/1‑236 (MQ=255)
                                                                                                                                                               gCATCGGGCACCTGCGCACGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTTCGGGTGGCGCGAGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                                  <  1:253905/239‑1 (MQ=255)
                                                                                                                                                               ccATCGGGCATCTGCGCACGGCTTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                                  <  1:956119/238‑1 (MQ=255)
                                                                                                                                                                cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGCGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTAGCATCAGGCATCTGCGCACGCCTGCCGGATGCGGCGAAAACGCCTTACCCGGCCTACGGCTGGCGCGGAGATTTGTAGGCATGATAAGAAGCGTTAGCGTCGGATACAGCATCTGccaccgcccgacggaggggataaa                                                                 >  2:253905/1‑216 (MQ=11)
                                                                                                                                                                         tcccagcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCGGCACACGACTGCCGGATGCGATAAACGCCTTGt                                                         >  1:932481/6‑238 (MQ=255)
                                                                                                                                                                          cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGTCCTGATAAGACGCGTTAGCGTTGCAACAGGCATCTGCACCCAACAGCCGGAGGCGACAAACGCCTTGTcc                                                       >  2:1255865/1‑239 (MQ=14)
                                                                                                                                                                          cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                        <  1:668495/238‑1 (MQ=255)
                                                                                                                                                                                aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATCGCGCCGGCATTTGTCAGCCTGATAAGACGCGGTAACGTCCACTCCGGCCACTCCCCCCCCCCGCCGGATGCCATAAAAGTATTTTCCgggcc                                                  >  2:620853/1‑235 (MQ=11)
                                                                                                                                                                                     tGCCGGTGGGGGAGTGAACGTCTTATCCGGCGTCGGGGTGGCGCGAGAATTTGTAGGGCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGGGGCGGAAACGCCTTATCCGGCCTACGGATGGCGGGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                            <  2:190341/239‑1 (MQ=14)
                                                                                                                                                                                     tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGCCTTTTCCGGCCTACAtt                                             >  1:1018734/1‑238 (MQ=255)
                                                                                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTc                                           >  1:1155726/1‑238 (MQ=255)
                                                                                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGAACGACTGCCGGATGCGGGGTAACAGCCTTATCCGGCCTACGGATCGGCCTGGGATTAGTAGGCCTGAGAAGATCCGGTACCGCCGCACACGGACTCTGCACACGACCGCCCGAGGCGATAAAAGTCTTGTCCGGcccaaaatac                                           >  2:304471/1‑230 (MQ=14)
                                                                                                                                                                                                 tgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCTACTGCGCACGAGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                 <  1:676369/237‑1 (MQ=255)
                                                                                                                                                                                                 gcggAACCCCTTATCCGGCCTACGGGTGGCCGGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                 <  1:552341/235‑1 (MQ=255)
                                                                                                                                                                                                   tGAACGCCTTAACCGGGCCACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                               <  1:443768/238‑1 (MQ=255)
                                                                                                                                                                                                      cccctttttCCGGTCTGGGGGTGGGGGGGGAATTTGTAGGCCTGATAAGAGGGGGTAGGGGCGCATCAGGCATCTGCGACCGACTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                            <  2:1224759/230‑1 (MQ=17)
                                                                                                                                                                                                      aCGGCTTAGCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                            <  1:822895/238‑1 (MQ=255)
                                                                                                                                                                                                      aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGTCATCTGCACACGCCTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTTGCCCGTAGGCAGTCAt                           >  1:1231537/1‑239 (MQ=255)
                                                                                                                                                                                                          cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                       <  1:1255865/239‑1 (MQ=255)
                                                                                                                                                                                                           tttcccGGGCTTCGGGTGGCGCGGGAATTTGGAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                       <  1:1251375/234‑1 (MQ=255)
                                                                                                                                                                                                                  gCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGCCTGCCGGATGCGATAAACTTCTTGTCCTGCCTACATTTCGCCCGTAGGCGGTCATTAAATAGTTTTg               >  1:908779/1‑239 (MQ=255)
                                                                                                                                                                                                                    cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                        >  1:1128952/1‑196 (MQ=255)
                                                                                                                                                                                                                    cTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                        <  2:1128952/196‑1 (MQ=255)
                                                                                                                                                                                                                      aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTAGTCCTGATAAGACGCGTTATCGTCGCATCAGCCATCTGCACACGACGGCCGGATGCCATAAACCTCTTGTCCGGCCTACCTTTCCCCCGTAgtccacccttaaatatttttaatt            >  2:340280/1‑215 (MQ=14)
                                                                                                                                                                                                                        gggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGTTAAACGTCTTGTCCGGCCTCAATTTCGCCCGTAGGCAGTCCTTAAATGGTTCTGTTTAcc         >  1:724576/1‑238 (MQ=255)
                                                                                                                                                                                                                        gggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCAGCAGGCATCTGCGAACGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg         <  1:340280/239‑1 (MQ=255)
                                                                                                                                                                                                                         ggTGGCGGGGGAATTTGGAGGCCTGAGAAGGCGGGTTAGGTGCCGATCAGGCATCTGCGCACGTCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg         <  2:265211/238‑1 (MQ=18)
                                                                                                                                                                                                                            ggggggagaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCgg      <  1:1252399/233‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc     <  1:582353/239‑1 (MQ=255)
                                                                                                                                                                                                                             gggcgGGAATTTGTAGGCCCGATAAGGCGGGGTAGCGTCGCATCAGGGATCGGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc     <  1:1343364/236‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGCTAAGACGCGTTAGCGTCGCATCCGGCATCTGCACCCCACTGCCGGATGGGATAAACGGCTTTTCCGGCTTCAATTTCGCCCGTAGGCATGCATTAAATCGTTTTGTTTatcgg       >  2:403104/1‑232 (MQ=17)
                                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCCCACGACTGCCGGATGCGGTAAACGTCTTGTCCGGCCTCCATTTGGCCCGTAGGAGGCCATTAAATAGTTCTGATTAtgcg       >  1:760855/1‑236 (MQ=17)
                                                                                                                                                                                                                                ggagaATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTgg  <  1:56128/237‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
CTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGG  >  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: