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.99.0% 888.2 / 1.6 196intergenic (‑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 . (83/111);  minor base G (0/2);  total (83/113)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

GGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGT  >  NC_000913_3_bme_pgi/4235308‑4235502
                                                                                                   |                                                                                               
ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                 <  1:1302815/100‑1 (MQ=255)
ggAACTCTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                 <  1:3786648/100‑1 (MQ=255)
 gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGGCGCGACGCTGCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                <  2:3189089/100‑2 (MQ=255)
  aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatc                                                                                               <  1:1141840/100‑3 (MQ=255)
  aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatc                                                                                               <  1:3516327/100‑3 (MQ=255)
      tttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                          >  1:991862/1‑100 (MQ=255)
       ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                          <  2:1261554/99‑1 (MQ=255)
       ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                          <  1:3419962/99‑1 (MQ=255)
         gCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgtc                                                                                        <  1:339422/99‑1 (MQ=255)
          ccGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgtc                                                                                        <  1:2320202/98‑1 (MQ=255)
          ccGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGa                                                                                      <  1:3118443/100‑1 (MQ=255)
           cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGat                                                                                     <  1:4199017/100‑1 (MQ=255)
           cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGCCGa                                                                                      >  2:1699442/1‑99 (MQ=255)
           cGGTTGTATCGGGCAAGGAAGGGGCGAGCGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgtc                                                                                        >  1:3123769/1‑97 (MQ=255)
            ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    <  2:3320033/100‑1 (MQ=255)
             gTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    <  2:2739801/99‑1 (MQ=255)
              ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATg                                                                                  <  1:4549135/100‑1 (MQ=255)
                 tATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGTGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTa                                                                                <  2:2757577/99‑1 (MQ=255)
                  aTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgt                                                                                         >  2:719053/1‑89 (MQ=255)
                   tCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgtc                                                                                        <  1:720299/89‑1 (MQ=255)
                     gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGTGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                           >  1:4631976/1‑100 (MQ=255)
                     gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                           >  1:228899/1‑100 (MQ=255)
                      ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCgg                                                                          <  1:2778034/100‑1 (MQ=255)
                      ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCgg                                                                          <  2:2827095/100‑1 (MQ=255)
                      ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                           >  2:2603675/1‑99 (MQ=255)
                      ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                           >  1:148723/1‑99 (MQ=255)
                        cAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGa                                                                         >  1:3301462/1‑99 (MQ=255)
                        cAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGa                                                                         >  2:1898387/1‑99 (MQ=255)
                          aggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATa                                                                       <  1:3104119/99‑1 (MQ=255)
                           ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAg                                                                     >  2:2996980/1‑100 (MQ=255)
                            gaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                    <  2:823503/100‑1 (MQ=255)
                               ggggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                 >  1:2796331/1‑100 (MQ=255)
                                gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                <  2:3674761/100‑1 (MQ=255)
                                 ggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                >  2:1447908/1‑99 (MQ=255)
                                 ggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                >  1:3799182/1‑99 (MQ=255)
                                 ggCGAGTGGTCGCGACGCTTCCACGCAAGGACTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                >  1:3710819/1‑99 (MQ=255)
                                  gCGAGTCGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                              <  1:3794369/100‑1 (MQ=255)
                                   cGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                             <  2:2308085/100‑1 (MQ=255)
                                     aGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                            <  1:2393436/99‑1 (MQ=255)
                                       tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                         <  2:2140816/100‑1 (MQ=255)
                                       tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                         <  1:2037119/100‑1 (MQ=255)
                                       tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                         <  1:2900087/100‑1 (MQ=255)
                                        ggTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                         >  1:3398986/1‑99 (MQ=255)
                                         gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCa                                                       >  2:2509796/1‑100 (MQ=255)
                                         gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCa                                                       >  1:4383292/1‑100 (MQ=255)
                                           cgcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                     >  1:1711643/1‑100 (MQ=255)
                                            gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                    >  2:1333057/1‑100 (MQ=255)
                                               aCGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                  <  2:1837470/99‑1 (MQ=255)
                                                cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                <  1:2367665/100‑1 (MQ=255)
                                                   ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             >  2:2188648/1‑100 (MQ=255)
                                                   ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             >  2:1587565/1‑100 (MQ=255)
                                                   ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             >  1:4427234/1‑100 (MQ=255)
                                                     ccACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             <  1:1066509/98‑1 (MQ=255)
                                                       aCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                         <  1:4067102/100‑1 (MQ=255)
                                                       aCGCAAGGGCTAGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                         <  2:3040714/100‑1 (MQ=255)
                                                        cGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                        >  1:3336131/1‑100 (MQ=255)
                                                        cGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                        >  1:3834739/1‑100 (MQ=255)
                                                        cGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                        >  1:4032555/1‑100 (MQ=255)
                                                         gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                       <  2:1114283/100‑1 (MQ=255)
                                                         gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                       <  1:1986380/100‑1 (MQ=255)
                                                         gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                       <  2:662353/100‑1 (MQ=255)
                                                             gggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATATGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                   <  2:419035/100‑1 (MQ=255)
                                                              ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                   <  2:3815720/99‑1 (MQ=255)
                                                              ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  <  1:641264/100‑1 (MQ=255)
                                                              ggCTTGTCGCGCATATGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  <  2:4646391/100‑1 (MQ=255)
                                                                   gTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                             >  1:2226396/1‑100 (MQ=255)
                                                                      gcgcATCTGCATGTGCGAAGGAAAGCTTGATCATCGTAGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                          <  1:3942089/100‑1 (MQ=25)
                                                                      gcgcATCTGCATGCGCGAAGGAAAGCTTGCTCATCGTCGATATGTACGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTg                                     <  2:1294902/89‑1 (MQ=37)
                                                                      gcgcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                           >  2:600367/1‑99 (MQ=255)
                                                                        gcATCTGCGTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGAcgcgacc                        <  2:4485469/100‑4 (MQ=25)
                                                                         cATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcg                       <  2:4356742/100‑1 (MQ=255)
                                                                           tCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                      >  2:2803097/1‑99 (MQ=255)
                                                                           tCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                      >  2:2448165/1‑99 (MQ=255)
                                                                            cTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                    >  2:1960380/1‑100 (MQ=255)
                                                                             tGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                   >  1:4273095/1‑100 (MQ=255)
                                                                              gCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                  >  1:1744908/1‑100 (MQ=255)
                                                                              gCATGCGCGAAGGAAAGCTTGATCATCGTCAATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                  >  1:543138/1‑100 (MQ=255)
                                                                               cATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                  <  1:1248370/99‑1 (MQ=255)
                                                                                aTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgca                <  2:3939245/100‑4 (MQ=21)
                                                                                   cgcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAtn              <  1:252122/99‑2 (MQ=255)
                                                                                   cgcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt               >  1:909255/1‑98 (MQ=255)
                                                                                   cgcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa             <  2:782004/100‑1 (MQ=255)
                                                                                     cgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg           <  1:2575767/100‑1 (MQ=255)
                                                                                     cgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg           <  1:1574825/100‑1 (MQ=255)
                                                                                       aaGGAAAGCTTGATCGTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGTCTGATGCGACGCGGTCGCGTCTTAt               <  2:4266853/94‑1 (MQ=18)
                                                                                       aaGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc         <  1:1019085/100‑1 (MQ=255)
                                                                                       aaGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc         <  1:1346944/100‑1 (MQ=255)
                                                                                       aaGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc         <  1:4469194/100‑1 (MQ=255)
                                                                                        aGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt        <  2:4030378/100‑1 (MQ=255)
                                                                                        aGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcat        <  1:153886/100‑3 (MQ=18)
                                                                                           aaaGCTTGATCATCGTCGATATGTAGGCCGGATCCGGCGTTCACGCCGCATCCGGGAACCGATGACTGATGCGACGCGGTCGCGTCTTCTCAGGCCTACa     <  2:3973762/100‑1 (MQ=25)
                                                                                           aaaGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa     <  2:3593116/100‑1 (MQ=255)
                                                                                           aaaGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa     <  2:327099/100‑1 (MQ=255)
                                                                                             aGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg   <  2:3886975/100‑1 (MQ=255)
                                                                                             aGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg   <  1:4523356/100‑1 (MQ=255)
                                                                                              gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt  >  1:2509934/1‑100 (MQ=255)
                                                                                              gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt  >  2:567126/1‑100 (MQ=255)
                                                                                              gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt  <  1:343021/100‑1 (MQ=255)
                                                                                              gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt  >  1:585452/1‑100 (MQ=255)
                                                                                              gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt  >  1:585611/1‑100 (MQ=255)
                                                                                                   |                                                                                               
GGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGT  >  NC_000913_3_bme_pgi/4235308‑4235502

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: