Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,746 A→G intergenic (‑513/‑2161) 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,233,7460AG99.0% 741.5 / ‑4.2 205intergenic (‑513/‑2161)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (105/98);  total (106/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 212 total aligned reads displayed.

GGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCG  >  NC_000913_3_bme_pgi/4233650‑4233842
                                                                                                |                                                                                                
 gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAgaggag                                                                                               >  1:3371073/1‑99 (MQ=255)
 gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAgaggag                                                                                               >  2:1485109/1‑99 (MQ=255)
 gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTCTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACCATTCTCAAAATCCgaggag                                                                                               >  2:4924499/1‑99 (MQ=255)
      cacTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGc                                                                                         <  1:756932/100‑1 (MQ=255)
       acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCGAAATCAGAGGAGTATTGc                                                                                         >  2:3429898/1‑99 (MQ=255)
       acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGc                                                                                         >  2:2637263/1‑99 (MQ=255)
          cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAAt                                                                                     >  1:4124988/1‑100 (MQ=255)
           aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAAt                                                                                     <  2:426317/99‑1 (MQ=255)
            aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                   <  1:1264366/100‑1 (MQ=255)
            aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                   >  1:3453066/1‑100 (MQ=255)
            aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                    <  2:3447434/99‑1 (MQ=255)
             cATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGCAGTATTGCTAATGGc                                                                                  <  2:5391222/100‑1 (MQ=255)
             cATTACGCTAACGGCACTAAAACCATCACAGTTTGCTGGGACTGGCGCTACGATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGc                                                                                  <  2:4112009/100‑1 (MQ=255)
              aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGc                                                                                  <  2:4723642/99‑1 (MQ=255)
              aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGc                                                                                  >  2:2219104/1‑99 (MQ=255)
                  cGCTAACGGCACTAAAACCGTCACATCTTTCTGTGACTGCCACTACAATCTTCCAAAGTCACAATTCGCAAAATCAGAGGAGTATTACTAGTGGCaa                                                                                >  1:5076709/1‑97 (MQ=255)
                  cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCaa                                                                                >  1:4543300/1‑97 (MQ=255)
                  cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCaa                                                                                >  1:290002/1‑97 (MQ=255)
                  cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAaga                                                                              <  2:1442859/99‑1 (MQ=255)
                  cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGccgctc                                                                             <  1:4366459/100‑6 (MQ=255)
                  cGCTAACGGCACTAAAACCATCACATTTTGCTGTGACTGGCGCTAGAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTg                                                                                          <  2:1732547/87‑1 (MQ=255)
                    cTAACGGCACTAAAACCATCACATTTTTTTTTTTTTTTCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAgaga                                                                            >  2:673650/1‑99 (MQ=255)
                    cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAgaga                                                                            >  2:1076451/1‑99 (MQ=255)
                      aaCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcg                                                                          >  2:1217647/1‑99 (MQ=255)
                       aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgcg                                                                        <  1:232893/100‑1 (MQ=255)
                       aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgcg                                                                        <  2:2295378/100‑1 (MQ=255)
                        cGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcga                                                                       >  1:3753986/1‑100 (MQ=255)
                        cGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcga                                                                       >  2:4864280/1‑100 (MQ=255)
                         ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                      <  2:1788786/100‑1 (MQ=255)
                         ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                      <  1:1252768/100‑1 (MQ=255)
                         ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                      <  1:1322815/100‑1 (MQ=255)
                         ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                      <  1:14721/100‑1 (MQ=255)
                           cACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgacga                                                                    <  1:1916601/100‑1 (MQ=255)
                             cTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAg                                                                  >  2:4173381/1‑100 (MQ=255)
                                aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                <  1:2196411/99‑1 (MQ=255)
                                aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                <  2:5292863/99‑1 (MQ=255)
                                aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                <  1:4784444/99‑1 (MQ=255)
                                  aCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGa                                                             <  1:359126/100‑1 (MQ=255)
                                   ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTg                                                              >  2:4677084/1‑98 (MQ=255)
                                   ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                            <  2:1652335/100‑1 (MQ=255)
                                   ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                            <  1:3259084/100‑1 (MQ=255)
                                     aTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGc                                                          <  1:2848017/100‑1 (MQ=255)
                                     aTCACATTTTTATGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGc                                                          <  2:3533542/100‑1 (MQ=255)
                                       cacaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCaa                                                        >  1:3811038/1‑100 (MQ=255)
                                        gcaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATGCTCATAAACAGAGGAGTATTGCTAATGGCAAGGGAGGCGACGAAGCTTGAAGCAAc                                                       <  2:4468405/99‑1 (MQ=255)
                                        acaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTTTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                       <  2:5319811/100‑1 (MQ=255)
                                          aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACgg                                                     <  2:635517/100‑1 (MQ=255)
                                          aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACgg                                                     <  2:2128919/100‑1 (MQ=255)
                                          aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACgg                                                     <  1:4225731/100‑1 (MQ=255)
                                           tttttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGt                                                    <  2:4100400/100‑1 (MQ=255)
                                              ttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCg                                                  <  2:7811/99‑1 (MQ=255)
                                               tCTGTGACGGGCGCTACAATCTTCCACAGTCACAATGCGCAAAATCAGAGGAGTATTGCTAGTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGc                                                 <  2:2155148/99‑1 (MQ=255)
                                                 tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCCACGGTCg                                                  >  2:2748420/1‑96 (MQ=255)
                                                 tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCa                                                >  1:583802/1‑98 (MQ=255)
                                                 tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGAATTGCTAATGGCAAGAGACGCAACGAAGCTTGAAGCAACGGTCGCa                                                >  2:861178/1‑98 (MQ=255)
                                                   tgACTGGCGCTACAATCTTCCAAAGTCCCAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAg                                             >  1:3244046/1‑99 (MQ=255)
                                                   tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                            >  1:4742473/1‑100 (MQ=255)
                                                     aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                           >  2:3590523/1‑99 (MQ=255)
                                                     aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTg                                          <  1:448033/100‑1 (MQ=255)
                                                         gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                      <  1:3988353/100‑1 (MQ=255)
                                                            cTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                    >  2:1982923/1‑99 (MQ=255)
                                                            cTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                    >  2:2256160/1‑99 (MQ=255)
                                                            cTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                    >  1:891534/1‑99 (MQ=255)
                                                               cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                <  1:4121278/100‑1 (MQ=255)
                                                               cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                >  1:1161745/1‑100 (MQ=255)
                                                                aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAATCAt                                 >  1:4636076/1‑98 (MQ=255)
                                                                aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                               <  1:2469191/100‑1 (MQ=255)
                                                                 aTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAATCATTgg                              >  1:4394825/1‑100 (MQ=255)
                                                                   cTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                            >  1:2542038/1‑100 (MQ=255)
                                                                    ttCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                            >  2:4958609/1‑99 (MQ=255)
                                                                    ttCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAAGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                           >  2:3733592/1‑100 (MQ=255)
                                                                     tCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                           >  1:2659776/1‑99 (MQ=255)
                                                                     tCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                           >  1:5288672/1‑99 (MQ=255)
                                                                     tCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                           >  2:3891282/1‑99 (MQ=255)
                                                                      ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGa                                                             <  2:2731584/64‑1 (MQ=255)
                                                                      ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                            <  2:1263015/97‑1 (MQ=255)
                                                                      ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                          >  1:3333686/1‑99 (MQ=255)
                                                                      ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                          <  2:478198/99‑1 (MQ=255)
                                                                      ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                           <  1:866370/98‑1 (MQ=255)
                                                                      ccAAAGTCACAATTCTCAAAAGCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                          <  2:477878/99‑1 (MQ=255)
                                                                       cAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCGTTGGGCgg                          >  2:5494159/1‑98 (MQ=255)
                                                                       cAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                          >  2:1498687/1‑98 (MQ=255)
                                                                       cAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGaa                        <  1:2096850/100‑1 (MQ=255)
                                                                        aaaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGACAAAGCATTGGGCGGAAt                       >  2:2875862/1‑100 (MQ=255)
                                                                          aGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCg                     <  2:2339672/100‑1 (MQ=255)
                                                                            tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                   >  2:3334566/1‑100 (MQ=255)
                                                                             cacaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAATGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                  >  1:1486910/1‑100 (MQ=255)
                                                                               caATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                >  1:5174583/1‑100 (MQ=255)
                                                                                aaTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg               <  1:5126131/100‑1 (MQ=255)
                                                                                    ctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGCAAAAGCATTGGGCGGAATCGGCCCcgctc            >  1:5157627/1‑97 (MQ=255)
                                                                                    ctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGCATCGGCCCcgcgc            >  2:3039160/1‑99 (MQ=255)
                                                                                     tcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc            >  2:418594/1‑98 (MQ=255)
                                                                                       aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata        <  1:1726498/100‑1 (MQ=255)
                                                                                       aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata        >  1:4481055/1‑100 (MQ=255)
                                                                                       aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGat         <  1:1596836/99‑1 (MQ=255)
                                                                                        aaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGatat       <  1:2287817/100‑1 (MQ=255)
                                                                                        aaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTGGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata        <  2:3290343/99‑1 (MQ=255)
                                                                                          aTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc     >  2:1346206/1‑100 (MQ=255)
                                                                                           tCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcg    >  1:1515060/1‑100 (MQ=255)
                                                                                           tCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcg    <  1:353502/100‑1 (MQ=255)
                                                                                                |                                                                                                
GGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCG  >  NC_000913_3_bme_pgi/4233650‑4233842

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