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,7460AG98.3% 170.5 / NA 59intergenic (‑513/‑2161)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (37/21);  minor base T (0/1);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCT  >  NC_000913_3_bme_pgi/4233599‑4233891
                                                                                                                                                   |                                                                                                                                                 
tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCGAAGTCACAATTCTCAAAATCagagga                                                                                                                                                 >  1:667342/1‑150 (MQ=255)
     tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCTCAATTCTCAAAATCAGAGGAGTAtt                                                                                                                                            >  2:235822/1‑150 (MQ=255)
       cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCACAGTCACAATTCTCAAAATCAGAGGAGTATTGc                                                                                                                                          <  1:621774/150‑1 (MQ=255)
            ggTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                                                                     >  1:559033/1‑150 (MQ=255)
              tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGAATGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGc                                                                                                                                   >  1:677539/1‑150 (MQ=255)
               cAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCa                                                                                                                                  >  2:368467/1‑150 (MQ=255)
                 ttagtTCAGCTAATAAATTCTTCACAGAGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                                                                                <  1:428270/149‑1 (MQ=255)
                 gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                                                                                <  1:11980/150‑1 (MQ=255)
                              taaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAg                                                                                                                   >  1:922390/1‑150 (MQ=255)
                               aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAATATCAGAGGAGTATTGCTAATGGCAAGAGACGcgacgacga                                                                                                                  >  1:52702/1‑147 (MQ=255)
                               aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGc                                                                                                                  <  2:24939/150‑1 (MQ=255)
                                aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                                                                 >  1:390149/1‑150 (MQ=255)
                                   gCTTCACTGCGCTAAGGGATTACACTCAACAAAACGCTAACCGCACTAAAACCATCACAATTTTCTGGGACTGGCGCTACAATCATCCAAAGTCACAATTCTCAAAACCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTGGa                                                                                                              <  2:723626/150‑1 (MQ=255)
                                      tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGAGCTAAAATCTTCCAAAGTCACAATTCTAAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGc                                                                                                           >  1:509641/1‑150 (MQ=255)
                                         cTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                                                                        <  2:339309/150‑1 (MQ=255)
                                            cgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGt                                                                                                     >  1:746867/1‑150 (MQ=255)
                                              atatgggATTACACTCAAGAATACGCTAACGGCACTACAACCATCACATTTTTCTGTGACTGTCCCTACAATCTTCCAATGTCACAATTCTCAAAATCAGAGGAGTTTTGCTAATGGCTAGAGACGCGACGAAGCTTGAAGCAACGGTCg                                                                                                   <  2:256947/146‑1 (MQ=255)
                                               tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGc                                                                                                  <  2:559033/150‑1 (MQ=255)
                                                   ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                              >  1:87036/1‑150 (MQ=255)
                                                         cacTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaa                                                                                        >  1:894260/1‑150 (MQ=255)
                                                           cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACTGTCGCAAAGCTGaaaaa                                                                                      >  2:923597/1‑150 (MQ=255)
                                                           cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                      >  2:471409/1‑150 (MQ=255)
                                                             cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCGAAACTCACAATTGTCAAAATCAGAGGAGTATTCCTAATGGCAAGAGACGCGACGAAGCTGGAACCAACGGTCGCAAAGCTGCAAAAGc                                                                                    >  1:318779/1‑150 (MQ=255)
                                                             cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                                                                    >  1:673003/1‑150 (MQ=255)
                                                             cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                                                                    >  1:5937/1‑150 (MQ=255)
                                                             cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATATTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGAAAGAGACGCTACGAAGCTTGAAGCAACGGTCTCAAAGCTGAAAAAGc                                                                                    >  1:479861/1‑150 (MQ=255)
                                                                   tACGCTATAGGCACTAAAACCATCTCATTTTTCTGTGACTGGCGCTATTATCTTCCAAAGTCACAATTCTCTAAATCAGAGGAGTTTTGCTAATTGCATGAGACGCGACGAAGCTTGATGCAACGGTCGCAAAGCTGAAAAAGCTTTggg                                                                              <  2:712883/150‑1 (MQ=255)
                                                                         aaCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                        >  2:865946/1‑150 (MQ=255)
                                                                           cGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTATTGCAAATAGACGCGACGAAGCTTGAAGCAACGCTCGAAAAGCTGATAAAGAATTGTGCgaatacg                                                                      >  1:607015/1‑144 (MQ=255)
                                                                             gCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGCGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                    <  1:87604/150‑1 (MQ=255)
                                                                                cTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                                 >  1:417170/1‑150 (MQ=255)
                                                                                   aaaCCATCACATATTTCTGTGAGTGGCGCTACAATCTTCCAAAGTCAAAATTCTCAAACTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACCGTCGCAATGCTGAAAAAGCTTTGGGCGGATTCGTCACcgcg                                                              >  1:158596/1‑150 (MQ=255)
                                                                                     aCCATCACATTTTTCTGTGACTGGCGCAACAATCTTCCAACGGCACAATTATCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGGTGAAAAAGCATTGGGCGGTATCGGCCCcgcgcg                                                            <  2:52702/150‑1 (MQ=255)
                                                                                      ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                           <  1:142398/150‑1 (MQ=255)
                                                                                        aTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata                                                         >  1:701040/1‑150 (MQ=255)
                                                                                          aacaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGCCGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATg                                                       <  2:722113/149‑1 (MQ=255)
                                                                                             aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgc                                                    >  1:760055/1‑150 (MQ=255)
                                                                                             aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgc                                                    >  1:764182/1‑150 (MQ=255)
                                                                                                     gtgACTGGCGCTACAATCTTCCAAAGTCACACTTCTCAAAATCCGAGGAGTTTTGCTAATGGCAAGAGACGAGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCTGCAtt                                            <  2:509641/150‑1 (MQ=255)
                                                                                                     gtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTGATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCTTTGGGCGGAATCGGCCCCGCGCGCTATGCGCGCGGGAtt                                            >  2:711398/1‑150 (MQ=255)
                                                                                                     gtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTCAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt                                            >  2:136352/1‑150 (MQ=255)
                                                                                                                   aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                              <  2:417170/150‑1 (MQ=255)
                                                                                                                      catCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGCCGCGACGAAGCTTGAAGCAACGGTCGCAAAGCGGAACAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggcc                           <  1:896043/148‑1 (MQ=255)
                                                                                                                      cTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggcc                           >  2:909531/1‑150 (MQ=255)
                                                                                                                              gTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                   >  1:536430/1‑150 (MQ=255)
                                                                                                                              gTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCCGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                   >  1:579970/1‑150 (MQ=255)
                                                                                                                               tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                  <  2:746867/150‑1 (MQ=255)
                                                                                                                                cacaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGACGGAATCTGCCAAGCGCAATAAGCGAGCCGCATACAGCGCCGCCCCCGGACGGTGCGGGcc                 >  2:829366/1‑149 (MQ=255)
                                                                                                                                  caATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGat               >  1:318534/1‑150 (MQ=255)
                                                                                                                                  aaaTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGat               <  1:368467/149‑1 (MQ=255)
                                                                                                                                   aaTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata              >  2:342929/1‑150 (MQ=255)
                                                                                                                                    aTTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAGGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGCCCCCGGCCGGTTCGGGCGatat             >  2:381494/1‑150 (MQ=255)
                                                                                                                                      tctcAAAATCAGAGGAGTATTGCTAAGGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTc           <  2:361818/150‑1 (MQ=255)
                                                                                                                                        tcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc         <  2:390149/150‑1 (MQ=255)
                                                                                                                                          aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                   >  2:364476/1‑138 (MQ=255)
                                                                                                                                          aaaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                   <  1:364476/138‑1 (MQ=255)
                                                                                                                                           aaaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGCCCCCGGCCGGTTCGGGCGATATTCGctct      >  1:585014/1‑150 (MQ=255)
                                                                                                                                              tCAGATGAGTATTGCTATTGGCAAGAGACGCTTCGTAGCTTGAAGCAAGGGTCGCAAAGCTGAAAAATCATTGGGCGGATTCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGcc   <  2:318779/150‑1 (MQ=255)
                                                                                                                                              tCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGcc   >  1:281707/1‑150 (MQ=255)
                                                                                                                                               cAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGAGCCGCCCCCGGACGGTTCGGGCGATATTCGCTCTGccg  >  2:438607/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCT  >  NC_000913_3_bme_pgi/4233599‑4233891

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