Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,772 G→T intergenic (‑539/‑2135) 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,7720GT100.0% 81.8 / NA 23intergenic (‑539/‑2135)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/8);  total (15/8)

CTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAA  >  NC_000913_3_bme_pgi/4233625‑4233918
                                                                                                                                                   |                                                                                                                                                  
cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACGTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctac                                                                                                                                                  >  2:71022/1‑150 (MQ=255)
          cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        >  1:137980/1‑150 (MQ=255)
            tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGc                                                                                                                                      <  2:56403/150‑1 (MQ=255)
                        gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCaaa                                                                                                                          >  1:64408/1‑150 (MQ=255)
                                       aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                                                                                                           >  2:85090/1‑150 (MQ=255)
                                                  ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                                                <  1:138178/150‑1 (MQ=255)
                                                   gCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                                               <  2:137980/150‑1 (MQ=255)
                                                    cACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                                                                                                        >  4:15566/1‑140 (MQ=255)
                                                    cACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                                                                                                        <  3:15566/140‑1 (MQ=255)
                                                                    tttttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcg                                                                              >  4:79684/1‑150 (MQ=255)
                                                                          tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAt                                                                        >  2:150133/1‑150 (MQ=255)
                                                                           gtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt                                                                       >  1:132739/1‑150 (MQ=255)
                                                                            tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTc                                                                      >  1:141759/1‑150 (MQ=255)
                                                                              aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAg                                                                    >  2:97604/1‑150 (MQ=255)
                                                                                                 aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                                 >  2:1005/1‑150 (MQ=255)
                                                                                                   aGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCgg                                               <  3:63698/150‑1 (MQ=255)
                                                                                                    gTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                                              <  1:85090/150‑1 (MQ=255)
                                                                                                    gTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                                              <  1:1005/150‑1 (MQ=255)
                                                                                                       acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                           >  3:102108/1‑150 (MQ=255)
                                                                                                              tcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCGCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc                                    <  1:97604/150‑1 (MQ=255)
                                                                                                                          gagTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgatg                        >  4:67723/1‑150 (MQ=255)
                                                                                                                                  cTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTg                >  1:108844/1‑150 (MQ=255)
                                                                                                                                         cAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACt         >  1:125745/1‑150 (MQ=255)
                                                                                                                                               aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGa   >  3:69782/1‑150 (MQ=255)
                                                                                                                                                cgctACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa  <  1:71022/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAA  >  NC_000913_3_bme_pgi/4233625‑4233918

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