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% 962.2 / 1.6 207intergenic (‑2174/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (106/98);  minor base G (0/2);  total (106/101)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

GAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGT  >  NC_000913_3_bme_pgi/4235309‑4235502
                                                                                                  |                                                                                               
gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                <  1:3324722/100‑2 (MQ=255)
gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                <  1:5438914/100‑2 (MQ=255)
gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                <  1:4182297/100‑2 (MQ=255)
gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                <  2:2420887/100‑2 (MQ=255)
 aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatc                                                                                               <  1:3753333/100‑3 (MQ=255)
      ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                          <  1:5324458/99‑1 (MQ=255)
      ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                          <  2:4857459/99‑1 (MQ=255)
      ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                          <  2:5169759/99‑1 (MQ=255)
       tGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgt                                                                                         <  2:1062984/99‑1 (MQ=255)
           ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    >  2:282698/1‑100 (MQ=255)
           ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    >  2:1190545/1‑100 (MQ=255)
           ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    >  2:1251498/1‑100 (MQ=255)
           ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    <  2:1574492/100‑1 (MQ=255)
           ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    >  2:630303/1‑100 (MQ=255)
           ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                    >  1:2644455/1‑100 (MQ=255)
             ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATg                                                                                  <  1:3918805/100‑1 (MQ=255)
                   cGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGcc                                                                            >  2:5845991/1‑100 (MQ=255)
                   cGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGcc                                                                            >  1:2936662/1‑100 (MQ=255)
                    gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGcc                                                                            <  2:4453754/99‑1 (MQ=255)
                    gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGcc                                                                            <  1:2060664/99‑1 (MQ=255)
                    gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                           <  2:2557740/100‑1 (MQ=255)
                    gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCACCTGCATGCGCGACGGAAAGCTTGATCATCGTCGCTATGTAGGCCg                                                                           >  2:4038485/1‑100 (MQ=255)
                     ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCg                                                                           >  1:5762913/1‑99 (MQ=255)
                      gCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGa                                                                         <  1:4200012/100‑1 (MQ=255)
                         aggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGAt                                                                        >  1:3596016/1‑98 (MQ=255)
                         aggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGAt                                                                        >  1:3899622/1‑98 (MQ=255)
                          ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAg                                                                     <  1:5013562/100‑1 (MQ=255)
                          ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAg                                                                     <  1:1685374/100‑1 (MQ=255)
                             aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                  >  2:1684319/1‑100 (MQ=255)
                               gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                <  2:1895586/100‑1 (MQ=255)
                               gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                 <  2:2221249/99‑1 (MQ=255)
                               gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                 <  2:4307145/99‑1 (MQ=255)
                                ggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                >  2:530035/1‑99 (MQ=255)
                                ggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                               <  2:3521764/100‑1 (MQ=255)
                                  cGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                             >  1:795274/1‑100 (MQ=255)
                                  cGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                             >  1:945500/1‑100 (MQ=255)
                                    aGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCGTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                           >  1:6043995/1‑100 (MQ=255)
                                      tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                          <  1:3044983/99‑1 (MQ=255)
                                      tGGGCGCGACGCTTCCACGCAAGGGCTTGTCGCGGATCTGCATGCGCGAAGGAAAGCTTGATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                          <  2:3692/99‑1 (MQ=38)
                                       ggTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                         >  1:3413647/1‑99 (MQ=255)
                                        gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCa                                                       >  1:2316321/1‑100 (MQ=255)
                                        gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACACCGCa                                                       >  1:5074653/1‑100 (MQ=255)
                                          cgcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                     >  2:5357845/1‑100 (MQ=255)
                                          cgcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                     >  1:3272967/1‑100 (MQ=255)
                                          cgcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                     >  1:5330705/1‑100 (MQ=255)
                                            cgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                   >  1:371052/1‑100 (MQ=255)
                                              aCGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                  <  2:1239279/99‑1 (MQ=255)
                                               cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                 >  1:3504073/1‑99 (MQ=255)
                                               cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                <  1:1993224/100‑1 (MQ=255)
                                                 cTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                              <  1:1079477/100‑1 (MQ=255)
                                                 cTTCCACGCAAGGGCTTGTCGCGCATCTGCATCCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                              >  2:4613470/1‑100 (MQ=255)
                                                  ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             >  1:5591174/1‑100 (MQ=255)
                                                   tCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             >  1:2137894/1‑99 (MQ=255)
                                                    ccACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             <  1:1362291/98‑1 (MQ=255)
                                                    ccACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                             <  1:877813/98‑1 (MQ=255)
                                                     cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                          >  2:2678044/1‑100 (MQ=255)
                                                     cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                          >  2:2353974/1‑100 (MQ=255)
                                                     cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                          >  2:879924/1‑100 (MQ=255)
                                                     cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                          >  1:1887031/1‑100 (MQ=255)
                                                     cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                          >  1:4633351/1‑100 (MQ=255)
                                                     cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                          >  1:4818987/1‑100 (MQ=255)
                                                      aCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                         >  1:149503/1‑100 (MQ=255)
                                                      aCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                         >  2:598049/1‑100 (MQ=255)
                                                        gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcn                                       <  2:3784027/100‑2 (MQ=255)
                                                        gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                        <  1:4534267/99‑1 (MQ=255)
                                                         cAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                      >  2:3187078/1‑100 (MQ=255)
                                                             ngCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  >  2:5506825/2‑100 (MQ=255)
                                                             ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  >  2:5795145/1‑100 (MQ=255)
                                                             ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  <  1:1847815/100‑1 (MQ=255)
                                                             ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  <  1:4289271/100‑1 (MQ=255)
                                                             ggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                  <  1:208851/100‑1 (MQ=255)
                                                             ggCGTGTCGGGCAGCTGCATGCGCGAAGGAAAGCTGGGTCATCGTCGGTAGGTAGGCCGGAGAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGagg                                  <  2:2696565/100‑3 (MQ=21)
                                                                   tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                            <  1:1483662/100‑1 (MQ=255)
                                                                   tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                             <  2:2110070/99‑1 (MQ=255)
                                                                   tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                             <  2:1480530/99‑1 (MQ=255)
                                                                    cgcgCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                           >  1:1642336/1‑100 (MQ=255)
                                                                      cgcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                         <  1:1255925/100‑1 (MQ=255)
                                                                      cgcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                         >  1:924830/1‑100 (MQ=255)
                                                                        cATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcg                       <  1:5383435/100‑1 (MQ=255)
                                                                           cTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                      >  2:3053056/1‑98 (MQ=255)
                                                                            tGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                    >  2:399287/1‑99 (MQ=255)
                                                                            tGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                   >  1:4221564/1‑100 (MQ=255)
                                                                               aTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                <  1:308425/100‑1 (MQ=255)
                                                                                tGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt               >  1:5709911/1‑100 (MQ=255)
                                                                                 gcgcGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt               >  2:2503595/1‑99 (MQ=255)
                                                                                  cgcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt               >  2:349081/1‑98 (MQ=255)
                                                                                  cgcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt               >  1:1909112/1‑98 (MQ=255)
                                                                                  cgcgAAGGAAAGCTTGATCAGCGTCGATGTGTAGGCCGGATAAGGCGTTCACGCCGCGTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                   <  2:1731645/94‑1 (MQ=21)
                                                                                    cgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg           <  1:676535/100‑1 (MQ=255)
                                                                                    cgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg           <  1:284944/100‑1 (MQ=255)
                                                                                    cgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg           >  1:409173/1‑100 (MQ=255)
                                                                                    cgAAGGAAAGCCTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg           <  1:3969326/100‑1 (MQ=255)
                                                                                     gAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc          >  1:5496826/1‑100 (MQ=255)
                                                                                     gAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc          <  2:1305314/100‑1 (MQ=255)
                                                                                     gAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc          >  1:1164258/1‑100 (MQ=255)
                                                                                      aaGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc         <  2:1396787/100‑1 (MQ=255)
                                                                                        ggAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatc       <  1:649088/100‑4 (MQ=18)
                                                                                         gAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa       >  1:1874893/1‑99 (MQ=255)
                                                                                            aGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg   <  1:4738269/100‑1 (MQ=255)
                                                                                             gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt  >  2:3360535/1‑100 (MQ=255)
                                                                                                  |                                                                                               
GAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGT  >  NC_000913_3_bme_pgi/4235309‑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: