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)
            ggAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAAcc                                                                                                                                          >  3:560561/1‑150 (MQ=255)
                tGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccctggc                                                                                                                                                 >  1:1551625/1‑136 (MQ=255)
                 gAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAAt                                                                                                                                     <  1:1080357/150‑1 (MQ=255)
                         ttAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCgg                                                                                                                             <  2:988648/150‑1 (MQ=255)
                             cGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGcccc                                                                                                                               <  2:1356208/144‑1 (MQ=255)
                             cGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGtt                                                                                                                         <  3:626565/150‑1 (MQ=255)
                              ggCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTg                                                                                                                        <  3:700149/150‑1 (MQ=255)
                                    ttCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCg                                                                                                                  >  1:316199/1‑150 (MQ=255)
                                     tCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAAAGAGCTAAGGCTT‑ACCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGa                                                                                                                 >  1:571218/1‑150 (MQ=255)
                                       aaTGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGata                                                                                                               <  1:767457/150‑1 (MQ=255)
                                        aTGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGatat                                                                                                              >  1:1236314/1‑150 (MQ=255)
                                             gATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATc                                                                                                         <  3:262789/150‑1 (MQ=255)
                                              aTACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCg                                                                                                        <  1:1375219/150‑1 (MQ=255)
                                                aCCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGtt                                                                                                      >  1:466639/1‑150 (MQ=255)
                                                    tCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  <  1:1193828/150‑1 (MQ=255)
                                                      aGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgc                                                                                                >  4:370228/1‑150 (MQ=255)
                                                      aGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTtgc                                                                                                >  2:119311/1‑150 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  <  2:1024811/143‑1 (MQ=255)
                                                           tACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  >  1:1024811/1‑143 (MQ=255)
                                                             cGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACg                                                                                         <  1:187921/150‑1 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGa                                                                                                                 >  1:1610240/1‑125 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                        >  1:1488341/1‑150 (MQ=255)
                                                              gggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACgg                                                                                        >  1:1356838/1‑150 (MQ=255)
                                                               ggCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                       <  3:370228/150‑1 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                                                                      >  1:531931/1‑150 (MQ=255)
                                                                gCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTACCGTCGATCTAAATCGTTTGTTGCTGCCACGGAg                                                                                      >  1:1202895/1‑150 (MQ=255)
                                                                 cTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGt                                                                                     <  4:527899/150‑1 (MQ=255)
                                                                  ttgttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCTCAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTa                                                                                    <  2:631005/150‑1 (MQ=255)
                                                                  ttgttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTa                                                                                    >  2:359628/1‑150 (MQ=255)
                                                                     ttgGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAgg                                                                                 >  2:372577/1‑150 (MQ=255)
                                                                        ggTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTTCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCa                                                                              >  4:714827/1‑150 (MQ=255)
                                                                          tCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATc                                                                            >  1:1504164/1‑150 (MQ=255)
                                                                           cAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGGCTCGGAGTAAGGGCATCa                                                                           <  2:842671/150‑1 (MQ=255)
                                                                            aGCCAGACACGCTGATACTCACGCAGGGCGTCCTGTACGAGGGCTACTTGCTAATACTCCCCAGCTCCGGCCT‑CCCTGGCACAAACACAAACACCCCGGCGCTGCCGTCGATATAAATCGTTTGTTGCCGCCACGGAGTAACAGCATCac                                                                          >  2:492885/1‑149 (MQ=255)
                                                                             gCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAat                                                                         >  1:1330039/1‑150 (MQ=255)
                                                                               cAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAatat                                                                       <  2:1504164/150‑1 (MQ=255)
                                                                                 gACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATcc                                                                     <  4:649894/150‑1 (MQ=255)
                                                                                   cacGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCcac                                                                   <  2:1615150/150‑1 (MQ=255)
                                                                                    acGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCcaca                                                                  >  1:964987/1‑150 (MQ=255)
                                                                                      gCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCacac                                                                 >  4:957069/1‑150 (MQ=255)
                                                                                      gCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacc                                                                >  3:440847/1‑150 (MQ=255)
                                                                                       cTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGATGATAATAACTAGCTACGGGTT‑CCCTGGTTCAACCACAATCTCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacca                                                               <  3:563303/150‑1 (MQ=255)
                                                                                        tGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                              >  2:1557238/1‑150 (MQ=255)
                                                                                             cTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCCCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCCCACCAACCAGc                                                         >  1:1041478/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACTATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  1:1007267/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                                    >  1:182969/1‑138 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                                    <  2:182969/138‑1 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  2:1539266/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  2:1198730/1‑150 (MQ=255)
                                                                                              tCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAAGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                        >  1:941956/1‑150 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTttgtt                                                                                                  <  2:854202/106‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:466639/150‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:1597372/150‑1 (MQ=255)
                                                                                                aCGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:1597503/150‑1 (MQ=255)
                                                                                                  gCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTgg                                                    <  3:546869/150‑1 (MQ=255)
                                                                                                   cAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                                   <  4:560561/150‑1 (MQ=255)
                                                                                                    aGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                       >  4:578195/1‑113 (MQ=255)
                                                                                                     gggCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAAt                                                 <  3:190447/150‑1 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                              >  3:988885/1‑150 (MQ=255)
                                                                                                        cGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGtt                                              <  2:316199/150‑1 (MQ=255)
                                                                                                          tCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga                                            >  3:380065/1‑150 (MQ=255)
                                                                                                             tGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaacg                                         >  2:717436/1‑150 (MQ=255)
                                                                                                               tACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTccct                                                                                                                                                    <  2:592092/42‑2 (MQ=255)
                                                                                                                 cgcgGGCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCAc                                                                                          >  4:499779/1‑97 (MQ=255)
                                                                                                                    gggCTTCTTGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                                  <  1:359628/150‑1 (MQ=255)
                                                                                                                       cttcttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAg                               >  2:649048/1‑150 (MQ=255)
                                                                                                                       cttcttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAg                               <  2:420182/150‑1 (MQ=255)
                                                                                                                          cttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaa                                                              <  2:770577/116‑1 (MQ=255)
                                                                                                                          cttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGAAGCAAGGGCTTCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGAt                            >  1:326626/1‑150 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacca                                                               <  1:234957/114‑1 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                                   <  2:1447719/126‑1 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGa                                                   >  1:1447829/1‑126 (MQ=255)
                                                                                                                           ttGCTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGCTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACacca                                                               >  2:234957/1‑114 (MQ=255)
                                                                                                                            tGCTGATAATAACGAGCTACGGC‑TGCTCTCGCTCAACCACAGGCGGCCCGTCGGTGCCGTCGAGATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCa                          <  3:354841/150‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCAc                         <  3:957069/150‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCa                                                                           <  1:826929/99‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                                  <  1:1099462/140‑1 (MQ=255)
                                                                                                                              cTGATAATAACGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGACGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGc                                  >  2:1099462/1‑140 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGtt                                                                                                      >  3:601942/1‑67 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCg                       <  2:1583163/146‑1 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCg                       >  1:1583163/1‑146 (MQ=255)
                                                                                                                                   aataaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGt                      <  2:1240134/147‑1 (MQ=255)
                                                                                                                                     taaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTAtgtg                 >  1:618825/1‑150 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                                                                    >  2:948297/1‑99 (MQ=255)
                                                                                                                                      aaCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGa                >  3:852111/1‑150 (MQ=255)
                                                                                                                                       aCGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaa                                           >  4:361837/1‑122 (MQ=255)
                                                                                                                                        cGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATc                                                                      >  4:615020/1‑94 (MQ=255)
                                                                                                                                        cGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      >  1:78861/1‑110 (MQ=255)
                                                                                                                                        cGAGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGtt                                                      <  2:78861/110‑1 (MQ=255)
                                                                                                                                          aGCTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAggtg            <  2:531931/150‑1 (MQ=255)
                                                                                                                                            cTACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGa                                                                                       >  2:1108928/1‑73 (MQ=255)
                                                                                                                                             tACGGCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGgtggtg         <  1:1163932/150‑1 (MQ=255)
                                                                                                                                                ggCTT‑CCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTgccg      >  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: ‑