Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,312,0810GC20.4% 545.0 / 130.2 225intergenic (‑38/+21)alsB/alsRD‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR
Reads supporting (aligned to +/- strand):  ref base G (88/91);  new base C (18/28);  total (106/119)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 235 total aligned reads displayed.

GCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTG  >  NC_000913_3_bme_pgi/4311982‑4312180
                                                                                                   |                                                                                                   
         gACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                             <  1:4532540/99‑1 (MQ=255)
          acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                           >  2:1534521/1‑100 (MQ=255)
           caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                          <  2:1260227/100‑1 (MQ=255)
           caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                           >  2:3491377/1‑99 (MQ=255)
            aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                          <  1:3161478/99‑1 (MQ=255)
            aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                         <  2:3004203/100‑1 (MQ=255)
            aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                         <  1:1781054/100‑1 (MQ=255)
              caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGTTATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                       <  2:3066166/100‑1 (MQ=255)
              caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                       <  1:1434545/100‑1 (MQ=38)
              caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                       <  1:772053/100‑1 (MQ=38)
              caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                       <  1:5055658/100‑1 (MQ=38)
              caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                       <  2:720584/100‑1 (MQ=255)
                ttAAGCCCACGAGTGTGCCGCTGACATATTTCAGATATTTATTCATACTAATGTTCTCACGGTAATGGGACAGGTAATTCCGCCATCGCTGAGTTACCtg                                                                                     <  2:3211689/100‑2 (MQ=37)
                ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                     <  1:2219808/100‑2 (MQ=38)
                ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                     <  1:2404616/100‑2 (MQ=38)
                ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtt                                                                                     >  2:4521058/1‑100 (MQ=255)
                  aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGCATTACCtttt                                                                                   >  2:1520146/1‑100 (MQ=255)
                   aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcc                                                                                  <  1:2702212/100‑5 (MQ=37)
                   aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                                   <  1:1109126/99‑4 (MQ=37)
                   aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                                   <  1:2828588/99‑4 (MQ=37)
                   aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                                   <  1:1516238/99‑4 (MQ=37)
                    gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                 >  2:4609532/1‑100 (MQ=255)
                     cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                >  2:1485282/1‑100 (MQ=255)
                     cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                <  1:1960205/100‑1 (MQ=255)
                      ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                >  1:3796302/1‑99 (MQ=255)
                       cACGAGTGTGCCTCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                              <  2:3442635/100‑1 (MQ=255)
                           aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                          <  1:142684/100‑1 (MQ=255)
                             tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                        <  1:835460/100‑1 (MQ=255)
                               tgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                        <  2:1768765/98‑1 (MQ=255)
                               tgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAATCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcg                                                                      >  1:4410799/1‑100 (MQ=255)
                                  cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                   >  2:403231/1‑100 (MQ=255)
                                   gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                                  >  1:1658596/1‑100 (MQ=255)
                                     tGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTg                                                                >  2:3304975/1‑100 (MQ=255)
                                       aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTg                                                                <  2:2797969/98‑1 (MQ=255)
                                        aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                             <  2:3717820/100‑1 (MQ=255)
                                         atatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                              >  1:601551/1‑98 (MQ=255)
                                          tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                           <  1:2760212/100‑1 (MQ=255)
                                           atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                           >  1:479494/1‑99 (MQ=255)
                                            tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                           >  2:561909/1‑98 (MQ=255)
                                            tttCAGATAGTTATTCATAATAATGTTCTCACGGTAAGGGGGCAGGTAACTCCGGGATGGCGGAATTACCTTTTTCTTATGTCAGCGCGCCTGGt                                                              <  2:2298999/95‑1 (MQ=255)
                                             ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                           >  1:3286357/1‑97 (MQ=255)
                                              tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                        >  1:2277250/1‑99 (MQ=255)
                                              tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                        <  2:3525860/99‑1 (MQ=255)
                                              tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                        >  2:5139211/1‑99 (MQ=255)
                                               cAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                       >  2:3266797/1‑99 (MQ=255)
                                               cAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                       <  2:1149770/99‑1 (MQ=255)
                                                cgATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGGGAAGaaa                                                     <  2:4066577/99‑1 (MQ=255)
                                                 gATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                    <  1:2095887/100‑1 (MQ=255)
                                                  atatTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                                   <  1:3989677/100‑1 (MQ=255)
                                                   tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                                  <  2:2212539/100‑1 (MQ=255)
                                                   tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                                   >  1:3582078/1‑99 (MQ=255)
                                                   tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                                   >  2:2880639/1‑99 (MQ=255)
                                                       tattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                              >  1:37590/1‑100 (MQ=255)
                                                        attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                             >  2:4632333/1‑100 (MQ=255)
                                                         ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                            >  1:211894/1‑100 (MQ=255)
                                                           cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                            >  2:1444904/1‑66 (MQ=255)
                                                           cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                          >  2:1586391/1‑100 (MQ=255)
                                                           cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                          >  2:2461095/1‑100 (MQ=255)
                                                           cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                          >  1:4753391/1‑100 (MQ=255)
                                                             taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                 >  2:4809943/1‑91 (MQ=255)
                                                             taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                        >  2:57481/1‑100 (MQ=255)
                                                             taataaTGTTCTCACGGTAAGGGGACAGCTAACTCCGCGATGGCGGAATCAGCTTCGTTTTATTTGGGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                          >  2:4775401/1‑98 (MQ=255)
                                                               ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                      >  1:1379525/1‑100 (MQ=255)
                                                                  atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtt                                   >  1:3019392/1‑100 (MQ=255)
                                                                   tGTTCTCACGGTAAGGGGACAGGTAACTCCACGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttt                                  <  2:2040997/100‑1 (MQ=255)
                                                                       ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                              >  1:3373256/1‑100 (MQ=255)
                                                                       ctcACGGTAAGGGGACAGGGATCTCCGCGATGTCGGAGTTACCTTTTTCTTATTTGAGCGCGCCTGGTGAGAAGTGATCAATAATTGCGCCGGtatttt                               <  2:3029850/99‑1 (MQ=255)
                                                                        tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                              >  1:2656583/1‑99 (MQ=255)
                                                                          aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                              <  2:821540/97‑1 (MQ=255)
                                                                          aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                             <  1:4232160/98‑1 (MQ=255)
                                                                            ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTTTTCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTCTTGCatat                         >  2:5015034/1‑100 (MQ=255)
                                                                              tAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTCCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatta                       <  2:2403646/100‑1 (MQ=255)
                                                                                 ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                    <  2:1903831/100‑1 (MQ=255)
                                                                                   ggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                  <  2:2777553/100‑1 (MQ=255)
                                                                                   ggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                  <  1:4803629/100‑1 (MQ=255)
                                                                                   aaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat                   >  1:350701/5‑99 (MQ=37)
                                                                                      cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc               <  2:4230689/100‑1 (MQ=255)
                                                                                      aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc               >  2:1209093/2‑100 (MQ=38)
                                                                                      aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc               >  1:354808/2‑100 (MQ=38)
                                                                                      aaGGTAATTCCGCCATCGCGGAGATACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTCTTTTGCATATTAATATTAGc               >  2:3822497/2‑100 (MQ=38)
                                                                                       aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCTTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              >  1:1621647/1‑100 (MQ=38)
                                                                                       aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCGTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              >  1:4086235/1‑100 (MQ=38)
                                                                                       aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              >  1:3707357/1‑100 (MQ=38)
                                                                                       aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              >  2:2580377/1‑100 (MQ=38)
                                                                                        ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt             >  1:2584381/1‑100 (MQ=38)
                                                                                        ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt             >  2:2809764/1‑100 (MQ=38)
                                                                                        ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt             >  1:4118709/1‑100 (MQ=255)
                                                                                        ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt             >  1:924780/1‑100 (MQ=255)
                                                                                         gTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg            >  1:1251492/1‑100 (MQ=38)
                                                                                         gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg            <  2:2254674/100‑1 (MQ=255)
                                                                                         gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg            <  2:2654017/100‑1 (MQ=255)
                                                                                          tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt           <  1:4292819/100‑1 (MQ=255)
                                                                                            aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc         <  2:1306456/100‑1 (MQ=255)
                                                                                             cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc         >  1:2786128/1‑99 (MQ=255)
                                                                                             cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc         >  2:249346/1‑99 (MQ=255)
                                                                                              tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCaa       >  1:2615463/1‑100 (MQ=255)
                                                                                                cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg     <  2:4945570/100‑1 (MQ=255)
                                                                                                 gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt    >  1:5268751/1‑100 (MQ=255)
                                                                                                 gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt    >  2:767279/1‑100 (MQ=255)
                                                                                                  cgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt    >  2:4752074/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
GCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTG  >  NC_000913_3_bme_pgi/4311982‑4312180

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