Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,235,408 Δ1 bp intergenic (‑2175/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,235,4070A.98.5% 586.5 / 3.0 135intergenic (‑2174/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (58/73);  minor base G (1/1);  total (60/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_bme_pgi/4235259‑4235550
                                                                                                                                                     |                                                                                                                                               
tCAATGCCTGCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                 <  2:247040/149‑1 (MQ=255)
 cAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                                                                <  1:562305/149‑2 (MQ=255)
 cAATGCCTTCGATCAATGCGTCGTGGAACTGGGTAAGGAACTGTCGACGGAACTTTTGCCAGTTGTATCCGGCAAGGAAGGGGCGACAGGACGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                                                                <  1:1180/149‑2 (MQ=255)
 cAATGCCTACGAGCAATGGGGCGTCGAACTGGGTAAGGCACGGGCGATGCAACTTTTGCATGGTCGATCTGGCAAGGAAGGGGCTAGTGTTCGAGACGCTTCCACGCAAGGGCTTGTCGCTCATCTGCATGCTCGAAGGAAAGCTTGat                                                                                                                                                <  2:35551/149‑2 (MQ=255)
   aTGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatca                                                                                                                                              >  2:122075/1‑146 (MQ=255)
     gCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAAGAAAGCTTGatcatc                                                                                                                                           >  1:22232/1‑149 (MQ=255)
      ccTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                                                                          <  1:258525/149‑1 (MQ=255)
       cTTCGAGCAATGGTGCGTGGAACTGGGTACTGACGTGGTAACTTTGCTTTTGCCGGTTTGATCGGGCAAGGAAGGGGCGAGTGGTCGCGCCGATTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgt                                                                                                                                         <  2:572855/149‑1 (MQ=255)
          cGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGa                                                                                                                                      <  1:462849/149‑1 (MQ=255)
           tattaaTGGGGCGAGGAACGGGGAAAGGAACTGGAGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGat                                                                                                                                     <  1:82285/145‑1 (MQ=255)
           gATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGGATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGat                                                                                                                                     <  1:318103/149‑1 (MQ=255)
            aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                                                                    >  1:53983/1‑149 (MQ=255)
               aaTGGGGGCTGGAAATGGGGAAGGAACCGGCGAAGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGt                                                                                                                                 <  1:284905/149‑1 (MQ=255)
               aaTGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAGGGGGCGAGTGGTCGCGACGGTGCCACGCAAGGGCTTGTCGCGCAAATGCATGCGCGAAGGAAAGCTTGATCAACGTCGATATGt                                                                                                                                 >  2:136557/1‑149 (MQ=255)
                aTGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTa                                                                                                                                >  2:566133/1‑149 (MQ=255)
                     gCGGGGAACTTGGAAAGTATATGGCGACGGAACTTTTGCTGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                                                                           <  1:72324/149‑1 (MQ=255)
                       gTGGCACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGACGAGGCGGGTGGTAGCGACGCTCTCACGCAAGGGCTTGTCGCGCATCTGAGAGCGCGAAGGGACGCTTGATCATCGTCGATATATAgacacga                                                                                                                         >  2:366698/1‑143 (MQ=255)
                             cTGGGTAATGAGCTGTCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                   <  1:135311/149‑1 (MQ=255)
                               gggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                 <  2:53983/149‑1 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                >  1:551019/1‑149 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  1:449376/149‑1 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACACTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGACCGGATAAGGCGtt                                                                                                                >  1:224408/1‑149 (MQ=255)
                                   aaGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                             <  2:297402/149‑1 (MQ=255)
                                    aGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      >  2:334280/1‑139 (MQ=255)
                                    aGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  1:334280/139‑1 (MQ=255)
                                     ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           <  2:548240/149‑1 (MQ=255)
                                        acttgtgACGGAATGTTTGCCGTTTGGTACGGGAAAGGAAGGGGCGAGTGGTCGCTACGCTTCCACGCAAGGGCTTGTCGCGCATTTGTATGCGCGGAGGATAGCTTGATCATCGTCGTTATGTTGGCCGGATAAGGCGTTCAcgccgc                                                                                                        <  2:219763/143‑1 (MQ=255)
                                          tGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  2:335947/139‑1 (MQ=255)
                                          tGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                >  1:335947/1‑139 (MQ=255)
                                           ggCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                     >  2:321811/1‑149 (MQ=255)
                                             cGACTGAACTTCTGCCGGTTTTGTCGGGGAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGTTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                   <  2:43771/149‑1 (MQ=38)
                                               aCGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                 >  2:495710/1‑149 (MQ=255)
                                                     cTTTTGCCGGTTGTATGGGGCAAGGAAGGGGCGTGTGGTCGCGACGCTTCCACGTAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                           <  2:142682/149‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                          <  1:320360/149‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                          >  1:378648/1‑149 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                          >  1:501516/1‑149 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAACGAt                                                                                          >  2:344214/1‑149 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGAAGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                           <  2:474864/148‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGAAAGGAAGGGGCGAGTGGTCGCGACGCTTCCAAGCAAGGGCATGACGCGAGTCTGCATGCGCGAAGGACAGCTTGCACATCGTCGATATGTGAGCCGGATCAGGCGTTGAGGCCGTATCCGGCAACCGAt                                                                                          >  2:321876/1‑149 (MQ=255)
                                                        ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                        <  2:49932/149‑1 (MQ=255)
                                                          gCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGGAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                      <  1:397350/149‑1 (MQ=255)
                                                               ttGTATCGGGCAAGCAAGGGGCGAGTGGTCGCTACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                 <  1:368278/149‑1 (MQ=255)
                                                                  tATCGGTCAAGGAAGGGGCGAGTGGTCGCGACTCTTCCACGCCAGGGCTTGTCGCGCATCTGAATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              <  2:152002/149‑1 (MQ=38)
                                                                   aTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                             <  1:293925/149‑1 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCCTCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAAGCgg                                                                          >  1:267519/1‑149 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                          >  1:307472/1‑149 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                          >  1:191529/1‑149 (MQ=255)
                                                                       ggaaaagaaggGGCTAGTGGTCGATACGCTTCCACGTAAGGGTTTGTCGCTTATCTGCATGCGAGAAGGAAAGCTTGACCAACCTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                         <  1:201875/143‑1 (MQ=21)
                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                         >  1:196097/1‑149 (MQ=255)
                                                                        caaacgaaggCGCGATCGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGAGAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        <  2:75639/144‑1 (MQ=38)
                                                                            ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGGCGCGt                                                                    <  2:551019/149‑1 (MQ=255)
                                                                             gaaggGGCGAGTGGTCGCGACGCTTCCACGCANGGGCTTNTCGCGCATNTGCATGCGCGAAGGAAAGCTTNATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   <  1:329032/149‑1 (MQ=255)
                                                                                ggggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:501490/1‑149 (MQ=255)
                                                                                 gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                               <  1:495710/149‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACTTTCCCACGTAAGGGCTTGTCGCGTTTCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATCTAGGCCGGATAAGGCGTTCACGCCGTATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             <  1:456060/149‑1 (MQ=37)
                                                                                    ggAGTGGTCGCGACGCTTCCACGCAAGTGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGTTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            <  1:99959/148‑1 (MQ=255)
                                                                                      ggTGGTCGCGACGCCTCCACGCAAGGGCTTGGCGCGGCTCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  2:477649/148‑1 (MQ=37)
                                                                                      aGTGGTCGCTACGTGTCCACGCAAGTGCTTTTCGCGCATGTGCATGCGGAAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  1:114449/149‑1 (MQ=25)
                                                                                       gTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCCGATGCGACGCGGTCGCGTCTTATCAGGcc                                                         <  1:497433/149‑1 (MQ=255)
                                                                                          ggtgcgAAGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                      <  1:366575/146‑1 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGGCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                      >  1:6549/1‑149 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              <  1:526416/124‑1 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              >  2:526416/1‑124 (MQ=255)
                                                                                            cgcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                    <  1:61709/149‑1 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                    >  2:161920/1‑84 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                    <  1:161920/84‑1 (MQ=255)
                                                                                              ttgcttttCCACGCAAGGGCTTGTTGCCGAGCGGGATGCGCGAAGGAGAGCCGCATCATCGTCGATATGTAAGCCAGAGAATGCTTTCACGCCGCATCCGGCAACCGATGCTTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  <  2:436854/143‑1 (MQ=18)
                                                                                                 cGCTTCCACGCAAGGGGTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               <  2:20300/149‑1 (MQ=255)
                                                                                                   cTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                             <  2:199792/149‑1 (MQ=255)
                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCTCCGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            <  2:337826/149‑1 (MQ=255)
                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAAGTAAGCTGGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGATCGATGc                                            >  1:143015/1‑149 (MQ=25)
                                                                                                     tCCACGCAAGGGTTTGTCGCGCATCTGCTTGTTGCTCAGAAACCTGGATCTTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                           <  1:315591/149‑1 (MQ=18)
                                                                                                      ccACGCAAGGGCTTTTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                          <  1:288393/149‑1 (MQ=255)
                                                                                                       caagCAAGGGTTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                         <  1:421095/146‑1 (MQ=255)
                                                                                                        aCGCAAGGGTTTTTCGGGTATCTGCATGCGCGAATGAAAGCTTGTTCATCGGCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                        <  1:274317/149‑1 (MQ=21)
                                                                                                        aCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                        <  1:163264/149‑1 (MQ=255)
                                                                                                          gCAAGGGTTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCACCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  1:569892/149‑1 (MQ=255)
                                                                                                          gCAAGGGGTTGTCCCGCATTTGCATGGTCGGACGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTTACGTCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  1:12749/149‑1 (MQ=21)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGTATGCGGGAAGGAAAGCGTGATCACCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  1:233517/149‑1 (MQ=25)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  1:363514/149‑1 (MQ=255)
                                                                                                                gCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGCTCAGCTACGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGCTAGGTACAt                                >  1:294379/1‑149 (MQ=25)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                              >  2:447318/1‑149 (MQ=255)
                                                                                                                     tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAACGTAt                           >  1:521094/1‑149 (MQ=255)
                                                                                                                      cgcgCATCTGCATGCGCGATGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAATGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                          <  1:566638/149‑1 (MQ=255)
                                                                                                                      cgcgCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGat                                                                                                                                     >  1:256795/1‑42 (MQ=255)
                                                                                                                       gcgcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCAATGCCGATATGTACAACGTATTc                         >  1:62546/1‑149 (MQ=255)
                                                                                                                       gcgcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTAGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTTCATCGTATTc                         >  1:18431/1‑149 (MQ=255)
                                                                                                                        cgcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                        <  2:294379/149‑1 (MQ=255)
                                                                                                                                cATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                               <  2:11978/102‑1 (MQ=255)
                                                                                                                                cATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                               >  1:11978/1‑102 (MQ=255)
                                                                                                                                 aTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               >  1:349272/1‑149 (MQ=255)
                                                                                                                                 aTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  1:365928/149‑1 (MQ=255)
                                                                                                                                   gcgcGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc             >  1:1754/1‑149 (MQ=255)
                                                                                                                                      gcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt           >  2:242994/1‑149 (MQ=255)
                                                                                                                                        gAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg         <  1:122075/149‑1 (MQ=255)
                                                                                                                                           ggAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCAc      >  1:291574/1‑149 (MQ=255)
                                                                                                                                            gAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATTTACATCGTATTCGGCAATTAATACATAGCACg     >  1:302161/1‑149 (MQ=255)
                                                                                                                                            gAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACg     >  2:72309/1‑149 (MQ=255)
                                                                                                                                              aaGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat   <  2:366900/149‑1 (MQ=255)
                                                                                                                                               aGCTTGATCATCTTCTATATGTCGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAACCGATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAGGTCGATGCCTATATGTACATCGTATTCGGCAATTAATACATAGCACgatt  <  2:535401/149‑1 (MQ=21)
                                                                                                                                                     |                                                                                                                                               
TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_bme_pgi/4235259‑4235550

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