Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 374,688 G→C 5.3% R274P (CGA→CCA) ‡ mhpE → 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,6880GC 5.3% 464.3 / 13.0 132R274P (CGA→CCA) ‡mhpE4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I
Reads supporting (aligned to +/- strand):  ref base G (63/62);  new base C (3/4);  total (66/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 134 total aligned reads displayed.

AAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGC  >  NC_000913/374420‑374959
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
aaaGCAGGGCTGAAACCAGAAACGGCAACTGGCAGGCACGCTCACCATAACCGGAGTCTTGGCGTGGCGAACGCTATGGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAAGGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGa                                                                                                                                                                                                                                                                            <  4:255094/275‑1 (MQ=255)
                   aaaCGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGGGTCGGTCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                         >  1:122638/1‑275 (MQ=255)
                      cGCAAACTGGCATGCACGCTCACCAGAACCTGAGTCTTGGCGGGGCGAAATCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTggg                                                                                                                                                                                                                                                       <  1:155836/274‑1 (MQ=255)
                                          cACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGGCGATCGCGAAACGCTGGGTCGGGGATCCGCCGGGGTTTACTCGAg                                                                                                                                                                                                                                  >  4:338042/1‑275 (MQ=255)
                                                    tGAGTCTTGGCGTGGCGACCTCCATCGGCGCGGTGGAAGGTGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGt                                                                                                                                                                                                                        <  2:208706/275‑1 (MQ=255)
                                                                                  cggTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTTAAACGGCGGCGGctcggttgggc                                                                                                                                                                                          >  4:319211/1‑268 (MQ=255)
                                                                                    gTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCCCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGACCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCtt                                                                                                                                                                                        >  2:508704/1‑275 (MQ=255)
                                                                                            gggCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTCAGTGCgg                                                                                                                                                                                >  2:278909/1‑275 (MQ=255)
                                                                                            gggCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCg                                                                                                                                                                                 <  2:243235/274‑1 (MQ=255)
                                                                                              gCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTCATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCCAGCTTCCTGCGTCACTGTGAAACGGCGGCGGGGCGTTATGGCTTAAGTGCGTTg                                                                                                                                                                              >  1:161012/1‑275 (MQ=255)
                                                                                                        aaTCGAGGCCAGCCTCGCGGGAAGGGGCGCGGTCGCAGTTAACGAACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTATATGCGTTAATGTATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCg                                                                                                                                                                    <  2:308269/275‑1 (MQ=255)
                                                                                                                             aaTGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTc                                                                                                                                                                                                                                                                              <  6:494358/148‑1 (MQ=255)
                                                                                                                               tGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGa                                                                                                                                                                                                                                                                            <  5:844385/148‑1 (MQ=255)
                                                                                                                               tGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGAACGGNACGAGTCGATCGCGCAACGATGGCGCTGGGGAGCGCGGGTGTTTAGTCGAGCGTACTGCGTCACTGTGAACCGGNGGCGGGGCGGTATGGACTAGCTGCGGTGGGNNATCTGGTGGAGTCGGGGCaagggaggggg                                                                                                                                             >  3:78937/1‑266 (MQ=255)
                                                                                                                                   ggcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTg                                                                                                                                        <  2:73777/274‑1 (MQ=255)
                                                                                                                                       ggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                          >  8:579550/1‑142 (MQ=255)
                                                                                                                                       ggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                          <  7:579550/142‑1 (MQ=255)
                                                                                                                                            aGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                                <  2:338043/274‑1 (MQ=255)
                                                                                                                                            aGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                                <  2:437860/274‑1 (MQ=255)
                                                                                                                                               tAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcg                                                                                                                                                                                                                                                            >  7:421167/1‑148 (MQ=255)
                                                                                                                                                  cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaa                                                                                                                                                                                                                                                                     >  7:237212/1‑136 (MQ=255)
                                                                                                                                                  cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaa                                                                                                                                                                                                                                                                     <  8:237212/136‑1 (MQ=255)
                                                                                                                                                  cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                         >  8:592437/1‑148 (MQ=255)
                                                                                                                                                               gggtTTATGGCCGCCGCGGATAAACGGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGt                                                                                                            <  2:171450/273‑1 (MQ=255)
                                                                                                                                                                tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                              >  5:586863/1‑145 (MQ=255)
                                                                                                                                                                 gtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                         >  6:677795/1‑133 (MQ=255)
                                                                                                                                                                   ttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                                                                                                                                                        <  7:592437/148‑1 (MQ=255)
                                                                                                                                                                   ttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                                                                                                                                                        >  7:1125410/1‑148 (MQ=255)
                                                                                                                                                                   ttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGgcg                                                                                                        <  2:518269/275‑1 (MQ=255)
                                                                                                                                                                       tGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    <  8:421167/148‑1 (MQ=255)
                                                                                                                                                                       tGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    >  5:1089860/1‑148 (MQ=255)
                                                                                                                                                                         ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    >  6:566487/1‑146 (MQ=255)
                                                                                                                                                                         ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    <  5:566487/146‑1 (MQ=255)
                                                                                                                                                                            ccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                                                                                                                                                                >  4:132615/1‑115 (MQ=255)
                                                                                                                                                                            ccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCtt                                                                                                                                                                                                                               <  5:818994/148‑1 (MQ=255)
                                                                                                                                                                                ggATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcg                                                                                                                                                                                                                                                            <  2:132616/115‑1 (MQ=255)
                                                                                                                                                                                  aTAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGcaac                                                                                         <  2:404448/275‑1 (MQ=255)
                                                                                                                                                                                   tAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTTCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTg                                                                                                                                                                                                                           <  8:821914/145‑1 (MQ=255)
                                                                                                                                                                                   tAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTg                                                                                                                                                                                                                           >  7:821914/1‑145 (MQ=255)
                                                                                                                                                                                          ggCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtgtg                                                                                                                                                                                                                 <  8:1125410/148‑1 (MQ=255)
                                                                                                                                                                                          ggCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGAATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGGTCTGCGCAACAACAaata                                                                                 >  2:100793/1‑275 (MQ=255)
                                                                                                                                                                                               gcagcaTGGGACCGATCTCTATGCGGTAAGGGAGGGTGCCGACGACCTGGGGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGGGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatga                                                                             <  4:287361/274‑1 (MQ=255)
                                                                                                                                                                                                  gcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAg                                                                                                                                                                                                                                                                                >  8:885494/1‑77 (MQ=255)
                                                                                                                                                                                                  gcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAg                                                                                                                                                                                                                                                                                <  7:885494/77‑1 (MQ=255)
                                                                                                                                                                                                  gcaGGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGGCCTGGTGCGGCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGAc                                                                         <  1:78937/275‑1 (MQ=255)
                                                                                                                                                                                                      gggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                         >  8:232966/1‑96 (MQ=255)
                                                                                                                                                                                                      gggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                         <  7:232966/96‑1 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACg                                                                                                                                                                                                                                                                                  <  7:786645/70‑1 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACg                                                                                                                                                                                                                                                                                  >  8:786645/1‑70 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcg                                                                                                                                                                                                                                                                         <  6:1100924/79‑1 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcg                                                                                                                                                                                                                                                                         >  5:1100924/1‑79 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCt                                                                                                                                                                                                                                                >  8:1137631/1‑104 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCt                                                                                                                                                                                                                                                <  7:1137631/104‑1 (MQ=255)
                                                                                                                                                                                                         aCCTATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTg                                                                                                           >  2:269535/1‑234 (MQ=255)
                                                                                                                                                                                                          ccGATCTCTATGCGTTAATGGAGGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCgag                                                                  <  2:433861/274‑1 (MQ=255)
                                                                                                                                                                                                              tctcTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                          >  5:978093/1‑71 (MQ=37)
                                                                                                                                                                                                              tctcTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                          <  6:978093/71‑1 (MQ=38)
                                                                                                                                                                                                               ctctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                                                                                                                                                                <  5:181461/80‑1 (MQ=255)
                                                                                                                                                                                                                 ctATGCGTTAATGGATGCCGCCGACGCCCGGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCa                                                          <  4:100792/275‑1 (MQ=255)
                                                                                                                                                                                                                  tATGCGTTAATGGAGGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCa                                                          <  2:319212/274‑1 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctgg                                                                                                                                                                                                                                                              >  8:887028/1‑78 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctgg                                                                                                                                                                                                                                                              <  7:887028/78‑1 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAg                                                                                                                                                                                                                                  >  8:651764/1‑106 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAg                                                                                                                                                                                                                                  <  7:651764/106‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtg                                                                                                                                                                                                                   <  5:928766/119‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtg                                                                                                                                                                                                                   >  6:928766/1‑119 (MQ=255)
                                                                                                                                                                                                                          aaTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGt                                                                                                                                                                                                                                                                               >  7:164148/1‑54 (MQ=255)
                                                                                                                                                                                                                          aaTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGt                                                                                                                                                                                                                                                                               <  8:164148/54‑1 (MQ=255)
                                                                                                                                                                                                                           aTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    <  6:1098468/96‑1 (MQ=255)
                                                                                                                                                                                                                                 gccgccGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGata                                                                                                                                                                          >  5:1024389/1‑148 (MQ=255)
                                                                                                                                                                                                                                  ccgccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    >  7:283312/1‑89 (MQ=255)
                                                                                                                                                                                                                                  ccgccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                    <  8:283312/89‑1 (MQ=255)
                                                                                                                                                                                                                                           acCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGa                                                                                                                                                                <  6:1024389/148‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGaaa                                                                                                                                                                                                              <  8:328480/99‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGaaa                                                                                                                                                                                                              >  7:328480/1‑99 (MQ=255)
                                                                                                                                                                                                                                               ggTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTTCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATGCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGCTATGATCGTTTACGTGCCGCTGGCTCTGCGCAACAACAAATCATGCTTACTCCCCAGAACGTGCATTTTGTTTGCCCATCCCAgacccaacccc                             >  4:341489/1‑264 (MQ=255)
                                                                                                                                                                                                                                                    gTCCGTGACAGGATCGACCGGTACGAGTCGAGCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATgc                       <  1:52817/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      ccGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCgg                                                                                                                                                                                >  7:147513/1‑121 (MQ=255)
                                                                                                                                                                                                                                                      ccGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCgg                                                                                                                                                                                <  8:147513/121‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                                                                                                                                                             <  7:921793/105‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTgg                                                                                                                                                          <  6:229661/138‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTgg                                                                                                                                                          >  5:229661/1‑138 (MQ=255)
                                                                                                                                                                                                                                                               gATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                              >  5:592096/1‑50 (MQ=255)
                                                                                                                                                                                                                                                               gATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                              <  6:592096/50‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACt                                                                                                                                                                                                                                      <  7:817670/56‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACt                                                                                                                                                                                                                                      >  8:817670/1‑56 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                                                                                                                                                             >  5:775160/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                                                                                                                                                             <  6:775160/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCCGGGCAAACGCCGGATGGTTgg                                                                                                                                       <  6:884385/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aCCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGCCTGCCGGGAACGTGCCGTTTGTATGCCCATCCCAGCACGACTGCCGGAAGCGGCGGTAACGCCTTAt       >  2:444361/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                     cgggTGCGAGTGGGTCGGGAAAGGGGGGGGCTGGGGTAGGCTGGGGGTTACGCGGGGTTCTTGAGGCAATGGGAAGCGGAGGCGGGGCGTTATGGCGTAAGGGCGGTGGGTATTCTCGTTGGGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGGTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAAGAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTAt       <  4:286551/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTGCGGGTCGTGCGCGAAACGCTGGGGCTGGGATACGCTGGGGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTGAGGGTTTAAGGGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGGGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATc      <  4:440910/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGc                                                                                                                                                                                                                                 <  8:311311/57‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGc                                                                                                                                                                                                                                 >  7:311311/1‑57 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGtt                                                                                                                                         >  7:146588/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGtt                                                                                                                                         <  8:146588/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTa                                                                                                                                                                                       >  7:1098081/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                         tACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTa                                                                                                                                                                                       <  8:1098081/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGCGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGCTGCGGCGTGAACGCCTTACCCGGc  >  2:7448/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                           cGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTg                                                                                                           >  2:567692/1‑168 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
AAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGC  >  NC_000913/374420‑374959

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