Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 374,683 G→C 5.6% P272P (CCG→CCC) ‡ 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,6830GC 5.6% 316.0 / 8.2 90P272P (CCG→CCC) ‡mhpE4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I
Reads supporting (aligned to +/- strand):  ref base G (42/43);  new base C (2/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

CCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCG  >  NC_000913/374536‑374828
                                                                                                                                                   |                                                                                                                                                 
ccTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCgg                                                                                                                                                  >  4:377838/1‑149 (MQ=255)
  tCGCGGAAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATAGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTa                                                                                                                                                <  4:193161/149‑1 (MQ=255)
    gcgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACg                                                                                                                                              >  4:817372/1‑149 (MQ=255)
     cgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                             <  1:250956/149‑1 (MQ=255)
     cgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                             <  2:297752/149‑1 (MQ=255)
            gggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                      <  1:71992/149‑1 (MQ=255)
              gcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcgc                                                                                                                                    >  1:93705/1‑149 (MQ=255)
              gcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcgc                                                                                                                                    >  2:16749/1‑149 (MQ=255)
                  gggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaaa                                                                                                                                <  4:547262/149‑1 (MQ=255)
                   ggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGTATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAc                                                                                                                               <  1:759560/149‑1 (MQ=255)
                      gcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                            >  2:665904/1‑149 (MQ=255)
                      gcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                            <  2:701409/149‑1 (MQ=255)
                       cAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctg                                                                                                                           <  2:93705/149‑1 (MQ=255)
                         ggTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggc                                                                                                                         >  4:223300/1‑149 (MQ=255)
                              cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctgg                                                                                                                    <  3:129518/149‑1 (MQ=255)
                               gCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGAGTCGTACCGGTCGATcct                                                                                                                                    >  3:855437/1‑130 (MQ=255)
                               gCACCGCTGGAAGTGTTTATCGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTggg                                                                                                                   <  2:529971/149‑1 (MQ=255)
                               gCACCGCTGGAAGTGATTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATcct                                                                                                                                    <  4:855437/132‑3 (MQ=255)
                                cACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                             <  2:40358/122‑1 (MQ=255)
                                cACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                             >  1:40358/1‑122 (MQ=255)
                                 aCCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATc                                                                                                                                      <  2:965128/128‑1 (MQ=255)
                                 aCCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATc                                                                                                                                      >  1:965128/1‑128 (MQ=255)
                                   cGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATcct                                                                                                                                    <  4:293024/128‑3 (MQ=255)
                                   cGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATcct                                                                                                                                    >  3:293024/1‑126 (MQ=255)
                                      tGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCt                                                                                                            >  1:48759/1‑149 (MQ=255)
                                       ggAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTg                                                                                                           <  3:240388/149‑1 (MQ=255)
                                         aaGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGgt                                                                                                         <  2:916030/149‑1 (MQ=255)
                                         aaGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGgt                                                                                                         <  1:377023/149‑1 (MQ=255)
                                            tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGtt                                                                                                      >  1:85230/1‑149 (MQ=255)
                                            tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGtt                                                                                                      <  2:928703/149‑1 (MQ=255)
                                              tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                    >  1:952437/1‑149 (MQ=255)
                                              tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                    >  3:450455/1‑149 (MQ=255)
                                              tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                    >  1:101788/1‑149 (MQ=255)
                                              tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                    >  4:157194/1‑149 (MQ=255)
                                              tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                    <  1:14043/149‑1 (MQ=255)
                                               ttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTAc                                                                                                   <  4:853273/149‑1 (MQ=255)
                                                   tGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGa                                                                                               <  1:321516/149‑1 (MQ=255)
                                                     ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGc                                                                                             <  1:446270/149‑1 (MQ=255)
                                                     ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGc                                                                                             <  4:243036/149‑1 (MQ=255)
                                                       gccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCtt                                                                                           >  1:317916/1‑149 (MQ=255)
                                                            ggATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGc                                                                                      <  3:292936/149‑1 (MQ=255)
                                                            ggATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGc                                                                                      >  4:557028/1‑149 (MQ=255)
                                                             gATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCg                                                                                     <  4:102173/149‑1 (MQ=255)
                                                                aaaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                             <  1:143570/90‑1 (MQ=255)
                                                                aaaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                             >  2:143570/1‑90 (MQ=255)
                                                                aaaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCa                                                                                  >  4:290566/1‑149 (MQ=255)
                                                                 aaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCAc                                                                                 >  1:277803/1‑149 (MQ=255)
                                                                 aaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCAc                                                                                 <  3:816975/149‑1 (MQ=255)
                                                                 aaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCAc                                                                                 >  4:195918/1‑149 (MQ=255)
                                                                    tGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtgt                                                                              <  1:69004/149‑1 (MQ=255)
                                                                      ggCTGTCAGCAGGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGa                                                                            >  4:672723/1‑149 (MQ=255)
                                                                         tGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctgg                                                                                                                    <  1:157824/106‑1 (MQ=255)
                                                                         tGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctgg                                                                                                                    >  2:157824/1‑106 (MQ=255)
                                                                         tGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAc                                                                         <  4:63457/149‑1 (MQ=255)
                                                                         tGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAc                                                                         <  2:774215/149‑1 (MQ=255)
                                                                                 tGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAcggcggcgg                                                                 >  3:925993/1‑149 (MQ=255)
                                                                                 tGGGACCGATCTCTATGCGTAAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAcggcggcgg                                                                 >  3:5603/1‑149 (MQ=255)
                                                                                       cGATCTCGATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGtt                                                           >  2:928680/1‑149 (MQ=255)
                                                                                         aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                         >  4:210005/1‑149 (MQ=255)
                                                                                           ctctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATgg                                                       <  1:72518/149‑1 (MQ=255)
                                                                                              tATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTAc                                                                                                   <  2:401256/102‑1 (MQ=255)
                                                                                              tATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTAc                                                                                                   >  1:401256/1‑102 (MQ=255)
                                                                                              tATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCtt                                                    >  1:45181/1‑149 (MQ=255)
                                                                                              tATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCtt                                                    >  3:713239/1‑149 (MQ=255)
                                                                                               aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTa                                                   >  1:448753/1‑149 (MQ=255)
                                                                                               aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTa                                                   >  2:508799/1‑149 (MQ=255)
                                                                                                   gTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTg                                               >  1:687194/1‑149 (MQ=255)
                                                                                                    ttAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGc                                              >  3:704553/1‑149 (MQ=255)
                                                                                                       aTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATgg                                                       <  1:161020/137‑1 (MQ=255)
                                                                                                       aTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATgg                                                       >  2:161020/1‑137 (MQ=255)
                                                                                                       aTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGt                                           <  3:512376/149‑1 (MQ=255)
                                                                                                          gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGa                                        <  3:186033/149‑1 (MQ=255)
                                                                                                            tGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGata                                      <  2:112178/149‑1 (MQ=255)
                                                                                                               cgccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCGGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtc                                   <  3:695652/149‑1 (MQ=255)
                                                                                                                 ccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                 <  3:210005/149‑1 (MQ=255)
                                                                                                                 ccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                 <  1:135797/149‑1 (MQ=255)
                                                                                                                 ccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                 >  1:1037622/1‑149 (MQ=255)
                                                                                                                 ccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                 <  2:906488/149‑1 (MQ=255)
                                                                                                                 ccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                 <  3:287292/149‑1 (MQ=255)
                                                                                                                  cgacgaCCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCg                                <  3:445819/149‑1 (MQ=255)
                                                                                                                   gacgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGt                               >  1:699773/1‑149 (MQ=255)
                                                                                                                         cTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCt                         <  1:859896/149‑1 (MQ=255)
                                                                                                                         cTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCt                         >  3:220621/1‑149 (MQ=255)
                                                                                                                           ggTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTgg                       <  3:377838/149‑1 (MQ=255)
                                                                                                                                gTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTCGAGCTGGGCaaa                  >  2:657712/1‑149 (MQ=255)
                                                                                                                                 tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAAc                 >  2:305974/1‑149 (MQ=255)
                                                                                                                                 tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAAc                 >  4:59367/1‑149 (MQ=255)
                                                                                                                                         aGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATg         >  1:536664/1‑149 (MQ=255)
                                                                                                                                         aGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATg         >  4:785539/1‑149 (MQ=255)
                                                                                                                                            aTCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGtt      <  1:665904/149‑1 (MQ=255)
                                                                                                                                                aCCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTggcg  >  4:1034227/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCG  >  NC_000913/374536‑374828

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