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.5% 739.7 / ‑4.8 203intergenic (‑513/‑2161)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (100/102);  total (101/102)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 215 total aligned reads displayed.

GGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCG  >  NC_000913_3_bme_pgi/4233649‑4233842
                                                                                                 |                                                                                                
 ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAgaggag                                                                                               >  1:4840436/1‑100 (MQ=255)
 ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAgaggag                                                                                               >  2:4802450/1‑100 (MQ=255)
   tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTa                                                                                             <  2:1629965/100‑1 (MQ=255)
       cacTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAGGTCACAATGCTCAAGATCAgaggag                                                                                               <  2:390728/94‑1 (MQ=255)
        acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCTAAGTCACAATTCTCAAAATCAGAGGAGTATTGCt                                                                                        <  2:381207/100‑1 (MQ=255)
        acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGc                                                                                         <  2:2325423/99‑1 (MQ=255)
         cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCCCAATTCTCAAAATCAGAGGAGTATTGCt                                                                                        >  1:4120053/1‑99 (MQ=255)
            aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                    <  1:1024384/100‑1 (MQ=255)
            aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                    >  2:1717689/1‑100 (MQ=255)
             aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                   <  2:286822/100‑1 (MQ=255)
             aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACAGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATgg                                                                                   >  1:5284975/1‑100 (MQ=255)
               aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGCGGAGTATTGCTAATGGc                                                                                  >  2:3044458/1‑99 (MQ=255)
                ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGc                                                                                  >  1:2955899/1‑98 (MQ=255)
                  aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                               <  1:5202777/99‑1 (MQ=255)
                  aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                               <  1:4219264/99‑1 (MQ=255)
                  aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                               <  1:312902/99‑1 (MQ=255)
                  aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGccgct                                                                              >  1:3830792/1‑96 (MQ=255)
                   cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCaa                                                                                >  2:4179119/1‑97 (MQ=255)
                   cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAaga                                                                              <  2:5281776/99‑1 (MQ=255)
                     cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAc                                                                           >  1:1752985/1‑100 (MQ=255)
                     cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAc                                                                           >  1:3387932/1‑100 (MQ=255)
                       aaCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACACTCTTCCAAAGTCACAATTCTCAAACTCAGAGGAGTATTGCTAATGGCAAGAGAcg                                                                          >  2:2100054/1‑99 (MQ=255)
                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgcg                                                                        <  1:4792068/100‑1 (MQ=255)
                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgcg                                                                        <  1:637204/100‑1 (MQ=255)
                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgc                                                                         <  2:378696/99‑1 (MQ=255)
                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgc                                                                         >  1:1323928/1‑99 (MQ=255)
                        aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgc                                                                         >  1:752515/1‑99 (MQ=255)
                        aCGGCACTAAAACCATCACATTTTTCTGTGACTAGCGCTACAATCTTCCAAAGTCGCAATTCTCTAAATCAGAGGAGTAGTGCTAATGTCAAGAGTcgcg                                                                        <  2:4418030/100‑1 (MQ=255)
                          ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                      <  2:2757390/100‑1 (MQ=255)
                           gCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGcgac                                                                      >  2:4826252/1‑99 (MQ=255)
                              cTAAAACCATCACATTTTTCTTTTACTTTCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAg                                                                  >  2:2079507/1‑100 (MQ=255)
                                aaaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTCCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                >  2:5058414/1‑100 (MQ=255)
                                 aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                <  1:1932782/99‑1 (MQ=255)
                                 aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                <  1:2672882/99‑1 (MQ=255)
                                   aCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTg                                                              >  1:5076949/1‑99 (MQ=255)
                                   aCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGa                                                             <  1:658175/100‑1 (MQ=255)
                                    ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTg                                                              >  2:1166895/1‑98 (MQ=255)
                                     cATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           <  2:3723222/100‑1 (MQ=255)
                                     cATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           <  2:2603295/100‑1 (MQ=255)
                                     cATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           <  2:346344/100‑1 (MQ=255)
                                     cATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           <  2:2485005/100‑1 (MQ=255)
                                     cATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           <  2:5212752/100‑1 (MQ=255)
                                      aTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           >  2:23946/1‑99 (MQ=255)
                                      aTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           >  1:4566223/1‑99 (MQ=255)
                                      aTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                           >  1:355595/1‑99 (MQ=255)
                                      aTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGc                                                          <  1:4660091/100‑1 (MQ=255)
                                       tCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCa                                                         <  1:5308803/100‑1 (MQ=255)
                                         acaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                       >  2:728475/1‑100 (MQ=255)
                                         acaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                       >  1:1418561/1‑100 (MQ=255)
                                          caTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACg                                                      >  1:697615/1‑100 (MQ=255)
                                               ttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCg                                                  <  2:1465198/99‑1 (MQ=255)
                                                 cTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCa                                                >  2:293300/1‑99 (MQ=255)
                                                    tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                            >  1:3623813/1‑100 (MQ=255)
                                                    tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                            <  2:1235362/100‑1 (MQ=255)
                                                    tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                            <  2:4119993/100‑1 (MQ=255)
                                                    tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                            <  2:4062594/100‑1 (MQ=255)
                                                    tgACTGGCGCTACAATCTTCCAAAGNCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCan                                               <  1:4713864/97‑2 (MQ=255)
                                                      aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                           <  2:3297140/99‑1 (MQ=255)
                                                      aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                           <  1:1844973/99‑1 (MQ=255)
                                                      aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                           <  2:1889267/99‑1 (MQ=255)
                                                      aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTg                                          <  1:2791447/100‑1 (MQ=255)
                                                      aCTGGCGCTACAAGCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                           <  2:1535010/99‑1 (MQ=255)
                                                         ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaa                                        >  1:2566052/1‑99 (MQ=255)
                                                          gcgcTACAATCTTCCAAAGTCACAATTGTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaa                                       >  1:5157689/1‑99 (MQ=255)
                                                          gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                      <  1:3161403/100‑1 (MQ=255)
                                                             cTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                   <  1:1625777/100‑1 (MQ=255)
                                                                cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                >  1:4390886/1‑100 (MQ=255)
                                                                cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                <  2:648785/100‑1 (MQ=255)
                                                                cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                 <  1:5121826/99‑1 (MQ=255)
                                                                cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                 <  2:4953633/99‑1 (MQ=255)
                                                                cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                 <  2:3064062/99‑1 (MQ=255)
                                                                 aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                               >  2:217225/1‑100 (MQ=255)
                                                                 aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                               <  1:2873697/100‑1 (MQ=255)
                                                                 aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                               >  1:2480550/1‑100 (MQ=255)
                                                                    cTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                            >  2:2384010/1‑100 (MQ=255)
                                                                       ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                            <  2:1411972/97‑1 (MQ=255)
                                                                       ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                          >  2:3635822/1‑99 (MQ=255)
                                                                       ccAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                           <  1:2223322/98‑1 (MQ=255)
                                                                           aGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCg                     <  2:89419/100‑1 (MQ=255)
                                                                           aGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCg                     <  1:156910/100‑1 (MQ=255)
                                                                            gTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                    <  1:1219522/100‑1 (MQ=255)
                                                                             tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                   <  1:3104043/100‑1 (MQ=255)
                                                                               acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGccc                 <  2:2539820/100‑1 (MQ=255)
                                                                                caATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                >  1:2188898/1‑100 (MQ=255)
                                                                                caATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                >  1:3733257/1‑100 (MQ=255)
                                                                                 aaTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg               >  2:3746514/1‑100 (MQ=255)
                                                                                 aaTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg               >  1:393296/1‑100 (MQ=255)
                                                                                    tctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc            >  1:2517749/1‑100 (MQ=255)
                                                                                    tctcAAAATCAGAGGAGTATTGCTAATCGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc            >  1:4046584/1‑100 (MQ=255)
                                                                                    tcccAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGAAGAAGCCTGAAGCTTCGGTCGCAAAGGTGAAAAAACATAGGGCGGAACCGGCCCggcgc            >  2:2324317/4‑100 (MQ=255)
                                                                                     ctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc            >  2:1018065/1‑99 (MQ=255)
                                                                                       cAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa          >  1:1998860/1‑99 (MQ=255)
                                                                                        aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata        <  1:1154479/100‑1 (MQ=255)
                                                                                        aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata        <  1:5269347/100‑1 (MQ=255)
                                                                                        aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata        <  2:2709954/100‑1 (MQ=255)
                                                                                         aaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGatat       >  2:4703749/1‑100 (MQ=255)
                                                                                          aaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATg      <  2:4261019/100‑1 (MQ=255)
                                                                                             cAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcg    >  2:3339799/1‑99 (MQ=255)
                                                                                              agaggaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcg  >  1:2033224/1‑100 (MQ=255)
                                                                                              agaggaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcg  >  1:4192867/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCG  >  NC_000913_3_bme_pgi/4233649‑4233842

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: