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

GCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGG  >  NC_000913/374535‑374829
                                                                                                                                                    |                                                                                                                                                  
gCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCg                                                                                                                                                    >  1:958466/1‑149 (MQ=255)
    cgcgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCTGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATTGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTAc                                                                                                                                                <  3:146187/149‑1 (MQ=255)
     gcgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACg                                                                                                                                               >  1:950981/1‑149 (MQ=255)
      cgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                              <  4:746009/149‑1 (MQ=255)
      cgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGa                                                                                                                                              <  4:604863/149‑1 (MQ=255)
       gggAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAg                                                                                                                                             <  3:1030879/149‑1 (MQ=255)
           aTGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGa                                                                                                                                         <  1:726540/149‑1 (MQ=255)
            tGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGAt                                                                                                                                        <  2:74301/149‑1 (MQ=255)
             gggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTGCGTTGCAGGATCGACCGGTACGAGTCGATc                                                                                                                                       >  1:716417/1‑149 (MQ=255)
             gggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATc                                                                                                                                       >  1:521572/1‑149 (MQ=255)
              ggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATcc                                                                                                                                      <  4:138452/149‑2 (MQ=255)
                            tAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                       >  2:86362/1‑134 (MQ=255)
                            tAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                       <  1:86362/134‑1 (MQ=255)
                             aaCGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgct                                                                                                                       >  3:1037312/1‑149 (MQ=255)
                               cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctgg                                                                                                                     <  2:850527/149‑1 (MQ=255)
                               cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctgg                                                                                                                     <  4:157801/149‑1 (MQ=255)
                                gCACCGCTGGAAGGGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGTTGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTggg                                                                                                                    >  4:214013/1‑149 (MQ=255)
                                 cACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGa                                                                                                                   <  2:521572/149‑1 (MQ=255)
                                       tGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCt                                                                                                             <  3:441280/149‑1 (MQ=255)
                                        ggAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTg                                                                                                            <  2:911821/149‑1 (MQ=255)
                                         gAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                           <  4:904568/149‑1 (MQ=255)
                                             tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGtt                                                                                                       <  2:11177/149‑1 (MQ=255)
                                               tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                     <  3:167421/149‑1 (MQ=255)
                                               tttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTa                                                                                                     >  3:745510/1‑149 (MQ=255)
                                                 tATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACt                                                                                                   >  4:533694/1‑149 (MQ=255)
                                                 tATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACt                                                                                                   <  3:956004/149‑1 (MQ=255)
                                                  aTTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTc                                                                                                  <  3:173615/149‑1 (MQ=255)
                                                  aTTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTc                                                                                                  <  2:854659/149‑1 (MQ=255)
                                                   ttGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                 >  1:92673/1‑149 (MQ=255)
                                                   ttGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                 >  2:17839/1‑149 (MQ=255)
                                                    tGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGa                                                                                                <  1:491259/149‑1 (MQ=255)
                                                     gccgccGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAg                                                                                               <  1:861187/149‑1 (MQ=255)
                                                        gccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCtt                                                                                            >  3:42332/1‑149 (MQ=255)
                                                           gcgGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCt                                                                                         <  3:557226/149‑1 (MQ=255)
                                                             ggATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGc                                                                                       <  2:804204/149‑1 (MQ=255)
                                                              gATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCg                                                                                      >  1:437247/1‑149 (MQ=255)
                                                                  aaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCAc                                                                                  >  3:60012/1‑149 (MQ=255)
                                                                   aCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTg                                                                                                            >  2:421028/1‑122 (MQ=255)
                                                                   aCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTg                                                                                                            <  1:421028/122‑1 (MQ=255)
                                                                       ggCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGgt                                                                                                          >  1:209932/1‑120 (MQ=255)
                                                                       ggCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGgt                                                                                                          <  2:209932/120‑1 (MQ=255)
                                                                            gcagcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAcgg                                                                        >  1:757350/1‑149 (MQ=255)
                                                                            gcagcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAcgg                                                                        >  1:80656/1‑149 (MQ=255)
                                                                                  tGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTc                                                                                           <  2:772669/124‑1 (MQ=255)
                                                                                  tGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTc                                                                                           >  1:772669/1‑124 (MQ=255)
                                                                                     gACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGgcgc                                                               >  3:618899/1‑149 (MQ=255)
                                                                                       ccGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGt                                                             <  3:630514/149‑1 (MQ=255)
                                                                                       ccGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCCGCGGCGCGt                                                             <  1:947045/149‑1 (MQ=255)
                                                                                        cGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGtt                                                            >  4:921755/1‑149 (MQ=255)
                                                                                         gATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTa                                                           <  2:950981/149‑1 (MQ=255)
                                                                                         gATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTa                                                           >  4:827557/1‑149 (MQ=255)
                                                                                          aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTAc                                                                                                                                                <  4:215252/63‑1 (MQ=255)
                                                                                          aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTAc                                                                                                                                                >  3:215252/1‑63 (MQ=255)
                                                                                          aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                          >  3:795991/1‑149 (MQ=255)
                                                                                          aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                          >  1:596310/1‑149 (MQ=255)
                                                                                          aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                          <  2:596310/149‑1 (MQ=255)
                                                                                          aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                          <  3:367773/149‑1 (MQ=255)
                                                                                           tctcTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATg                                                         <  1:477039/149‑1 (MQ=255)
                                                                                             tctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGc                                                       >  4:939444/1‑149 (MQ=255)
                                                                                              ctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCt                                                      >  3:579057/1‑149 (MQ=255)
                                                                                              ctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCt                                                      <  2:204599/149‑1 (MQ=255)
                                                                                                 tGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTaa                                                   <  4:1039219/149‑1 (MQ=255)
                                                                                                 tGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTaa                                                   <  4:135490/149‑1 (MQ=255)
                                                                                                    gTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTg                                                >  4:257347/1‑149 (MQ=255)
                                                                                                    gTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTg                                                >  1:770726/1‑149 (MQ=255)
                                                                                                    gTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTg                                                >  4:565521/1‑149 (MQ=255)
                                                                                                     ttAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGc                                               >  2:72422/1‑149 (MQ=255)
                                                                                                     ttAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACCGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGc                                               >  1:631843/1‑149 (MQ=255)
                                                                                                        aTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGt                                            >  1:41060/1‑149 (MQ=255)
                                                                                                           gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGa                                                                                                >  1:659236/1‑94 (MQ=255)
                                                                                                           gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGa                                                                                                <  2:659236/94‑1 (MQ=255)
                                                                                                           gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGa                                         <  3:214013/149‑1 (MQ=255)
                                                                                                            aTGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGat                                        <  2:958466/149‑1 (MQ=255)
                                                                                                                cgccgACGACCGGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGAGAAACGTTGGCGCTGGGATACGCTTGTGATTACTCAAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGTATATtc                                    <  4:392473/149‑1 (MQ=255)
                                                                                                                    gacgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGt                                <  3:939444/149‑1 (MQ=255)
                                                                                                                    gacgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGt                                >  2:49109/1‑149 (MQ=255)
                                                                                                                     acgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGTATCGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGCATATTCTCGtt                               <  4:100882/149‑1 (MQ=255)
                                                                                                                          cTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCt                          >  3:857816/1‑149 (MQ=255)
                                                                                                                            ggTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTgg                        <  1:17839/149‑1 (MQ=255)
                                                                                                                            ggTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTgg                        <  1:347902/149‑1 (MQ=255)
                                                                                                                             gTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTggg                       <  2:867170/149‑1 (MQ=255)
                                                                                                                             gTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTggg                       <  2:161565/149‑1 (MQ=255)
                                                                                                                              tGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGc                      <  1:702210/149‑1 (MQ=255)
                                                                                                                               gCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                  <  1:5237/136‑1 (MQ=255)
                                                                                                                               gCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                  >  2:5237/1‑136 (MQ=255)
                                                                                                                                  tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAc                                                                          <  3:1015788/93‑1 (MQ=255)
                                                                                                                                  tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAAc                                                                          >  4:1015788/1‑93 (MQ=255)
                                                                                                                                  tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGt                                                             <  4:98843/106‑1 (MQ=255)
                                                                                                                                  tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGt                                                             >  3:98843/1‑106 (MQ=255)
                                                                                                                                        gCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGa            <  4:42332/149‑1 (MQ=255)
                                                                                                                                          aGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATg          >  4:1033918/1‑149 (MQ=255)
                                                                                                                                          aGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATg          >  1:391338/1‑149 (MQ=255)
                                                                                                                                            gATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACACTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGtt                               >  4:279598/1‑126 (MQ=255)
                                                                                                                                            gATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACACTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGtt                               <  3:279598/126‑1 (MQ=255)
                                                                                                                                             aTCGACTCGTACCGGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCa                                                                                   >  4:1006760/1‑73 (MQ=38)
                                                                                                                                             aTCGACTCGTACCGGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCa                                                                                   <  3:1006760/73‑1 (MQ=38)
                                                                                                                                                gACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTggc    >  3:591379/1‑149 (MQ=255)
                                                                                                                                                  ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtgt                                                                               <  4:581990/72‑1 (MQ=255)
                                                                                                                                                  ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtgt                                                                               >  3:581990/1‑72 (MQ=255)
                                                                                                                                                  ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTggcgg  >  2:891215/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGG  >  NC_000913/374535‑374829

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: