Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,141,192 Δ1 bp coding (1030/2502 nt) frwA ← putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,141,1890C.96.8% 722.9 / 17.9 156coding (1033/2502 nt)frwAputative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base . (76/74);  total (78/78)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 171 total aligned reads displayed.

CTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGC  >  NC_000913_3_bme_pgi/4141040‑4141338
                                                                                                                                                     |                                                                                                                                                     
cTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTc                                                                                                                                                       >  3:715634/1‑150 (MQ=255)
cTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTc                                                                                                                                                       <  1:968559/150‑1 (MQ=255)
  gcgcTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccc                                                                                                                                                     >  2:187921/1‑150 (MQ=255)
   cgcTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccct                                                                                                                                                    <  4:633838/150‑2 (MQ=255)
   cgcTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccct                                                                                                                                                    >  1:1597372/1‑149 (MQ=255)
      ttCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccctggc                                                                                                                                                 <  3:752315/150‑5 (MQ=255)
            ggAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAAcc                                                                                                                                          >  3:560561/1‑150 (MQ=255)
                tGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccctggc                                                                                                                                                 >  1:1551625/1‑136 (MQ=255)
                 gAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAAt                                                                                                                                     <  1:1080357/150‑1 (MQ=255)
                         ttAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCgg                                                                                                                             <  2:988648/150‑1 (MQ=255)
                             cGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGcccc                                                                                                                               <  2:1356208/144‑1 (MQ=255)
                             cGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGtt                                                                                                                         <  3:626565/150‑1 (MQ=255)
                              ggCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTg                                                                                                                        <  3:700149/150‑1 (MQ=255)
                                    ttCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCg                                                                                                                  >  1:316199/1‑150 (MQ=255)
                                     tCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAAAGAGCTAAGGCTTACCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGa                                                                                                                 >  1:571218/1‑150 (MQ=255)
                                       aaTGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGata                                                                                                               <  1:767457/150‑1 (MQ=255)
                                        aTGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGatat                                                                                                              >  1:1236314/1‑150 (MQ=255)
                                             gATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATc                                                                                                         <  3:262789/150‑1 (MQ=255)
                                              aTACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCg                                                                                                        <  1:1375219/150‑1 (MQ=255)
                                                aCCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGtt                                                                                                      >  1:466639/1‑150 (MQ=255)
                                                    tCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  <  1:1193828/150‑1 (MQ=255)
                                                      aGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgc                                                                                                >  4:370228/1‑150 (MQ=255)
                                                      aGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgc                                                                                                >  2:119311/1‑150 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  <  2:1024811/143‑1 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  >  1:1024811/1‑143 (MQ=255)
                                                             cGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACg                                                                                         <  1:187921/150‑1 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGa                                                                                                                 >  1:1610240/1‑125 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                        >  1:1488341/1‑150 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                        >  1:1356838/1‑150 (MQ=255)
                                                               ggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                       <  3:370228/150‑1 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                                                                      >  1:531931/1‑150 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTACCGTCGATCTAAATCGTTTGTTGCTGCCACGGAg                                                                                      >  1:1202895/1‑150 (MQ=255)
                                                                 cTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGt                                                                                     <  4:527899/150‑1 (MQ=255)
                                                                  ttgttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCTCAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTa                                                                                    <  2:631005/150‑1 (MQ=255)
                                                                  ttgttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTa                                                                                    >  2:359628/1‑150 (MQ=255)
                                                                     ttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAgg                                                                                 >  2:372577/1‑150 (MQ=255)
                                                                        ggTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTTCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCa                                                                              >  4:714827/1‑150 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATc                                                                            >  1:1504164/1‑150 (MQ=255)
                                                                           cAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGGCTCGGAGTAAGGGCATCa                                                                           <  2:842671/150‑1 (MQ=255)
                                                                            aGCCAGACACGCTGATACTCACGCAGGGCGTCCTGTACGAGGGCTACTTGCTAATACTCCCCAGCTCCGGCCTCCCTGGCACAAACACAAACACCCCGGCGCTGCCGTCGATATAAATCGTTTGTTGCCGCCACGGAGTAACAGCATCac                                                                          >  2:492885/1‑149 (MQ=255)
                                                                             gCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAat                                                                         >  1:1330039/1‑150 (MQ=255)
                                                                               cAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAatat                                                                       <  2:1504164/150‑1 (MQ=255)
                                                                                 gACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATcc                                                                     <  4:649894/150‑1 (MQ=255)
                                                                                   cacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCcac                                                                   <  2:1615150/150‑1 (MQ=255)
                                                                                    acGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCcaca                                                                  >  1:964987/1‑150 (MQ=255)
                                                                                      gCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCacac                                                                 >  4:957069/1‑150 (MQ=255)
                                                                                      gCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacc                                                                >  3:440847/1‑150 (MQ=255)
                                                                                       cTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGATGATAATAACTAGCTACGGGTTCCCTGGTTCAACCACAATCTCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacca                                                               <  3:563303/150‑1 (MQ=255)
                                                                                        tGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                              >  2:1557238/1‑150 (MQ=255)
                                                                                             cTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCCCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCCCACCAACCAGc                                                         >  1:1041478/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACTATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  1:1007267/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                                    >  1:182969/1‑138 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                                    <  2:182969/138‑1 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  2:1539266/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  2:1198730/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAAGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  1:941956/1‑150 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  <  2:854202/106‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:466639/150‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:1597372/150‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:1597503/150‑1 (MQ=255)
                                                                                                  gCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTgg                                                    <  3:546869/150‑1 (MQ=255)
                                                                                                   cAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                                   <  4:560561/150‑1 (MQ=255)
                                                                                                    aGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                       >  4:578195/1‑113 (MQ=255)
                                                                                                     gggCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAAt                                                 <  3:190447/150‑1 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                              >  3:988885/1‑150 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                              <  2:316199/150‑1 (MQ=255)
                                                                                                          tCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga                                            >  3:380065/1‑150 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaacg                                         >  2:717436/1‑150 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccct                                                                                                                                                    <  2:592092/42‑2 (MQ=255)
                                                                                                                 cgcgGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAc                                                                                          >  4:499779/1‑97 (MQ=255)
                                                                                                                    gggCTTCTTGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                                  <  1:359628/150‑1 (MQ=255)
                                                                                                                       cttcttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAg                               >  2:649048/1‑150 (MQ=255)
                                                                                                                       cttcttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAg                               <  2:420182/150‑1 (MQ=255)
                                                                                                                          cttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                              <  2:770577/116‑1 (MQ=255)
                                                                                                                          cttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGAAGCAAGGGCTTCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGAt                            >  1:326626/1‑150 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacca                                                               <  1:234957/114‑1 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                                   <  2:1447719/126‑1 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                                   >  1:1447829/1‑126 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGCTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacca                                                               >  2:234957/1‑114 (MQ=255)
                                                                                                                            tGCTGATAATAACGAGCTACGGCTGCTCTCGCTCAACCACAGGCGGCCCGTCGGTGCCGTCGAGATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCa                          <  3:354841/150‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCAc                         <  3:957069/150‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCa                                                                           <  1:826929/99‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                                  <  1:1099462/140‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGACGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                                  >  2:1099462/1‑140 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGtt                                                                                                      >  3:601942/1‑67 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCg                       <  2:1583163/146‑1 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCg                       >  1:1583163/1‑146 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGt                      <  2:1240134/147‑1 (MQ=255)
                                                                                                                                     taaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTAtgtg                 >  1:618825/1‑150 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                                    >  2:948297/1‑99 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGa                >  3:852111/1‑150 (MQ=255)
                                                                                                                                       aCGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaa                                           >  4:361837/1‑122 (MQ=255)
                                                                                                                                        cGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATc                                                                      >  4:615020/1‑94 (MQ=255)
                                                                                                                                        cGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      >  1:78861/1‑110 (MQ=255)
                                                                                                                                        cGAGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:78861/110‑1 (MQ=255)
                                                                                                                                          aGCTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAggtg            <  2:531931/150‑1 (MQ=255)
                                                                                                                                            cTACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                       >  2:1108928/1‑73 (MQ=255)
                                                                                                                                             tACGGCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGgtggtg         <  1:1163932/150‑1 (MQ=255)
                                                                                                                                                ggCTTCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTgccg      >  3:260036/1‑150 (MQ=255)
                                                                                                                                                     ccccTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCcgccgc  <  1:1422547/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGC  >  NC_000913_3_bme_pgi/4141040‑4141338

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: