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,7460AG97.4% 121.7 / NA 38intergenic (‑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 (23/14);  minor base C (1/0);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGC  >  NC_000913_3_bme_pgi/4233604‑4233889
                                                                                                                                              |                                                                                                                                                
tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTAt                                                                                                                                            >  1:390599/1‑149 (MQ=255)
           agtagtTCAGCTAATACATCCTTCCCTGCGCTAAGGCTTCACCCTCACCCTTACGCTAACGGCACTAAAACCAACCCATTTTTCTGTCACTGGCGCTCCAATCTTCCAAAGTCACAATCCTCACAATCACACGAGCCTTCCTaacgccc                                                                                                                                 >  2:445807/1‑144 (MQ=255)
                   acgtaataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgc                                                                                                                         <  1:154132/146‑1 (MQ=255)
                   aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgc                                                                                                                         <  1:574458/149‑1 (MQ=255)
                   aGCTAATAAAGGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAcgc                                                                                                                         >  2:997500/1‑149 (MQ=255)
                         tcaaTGCTTCACTGCGGTAAGGGTTTACCCTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGGCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                                                                   <  1:973956/147‑1 (MQ=255)
                          aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAg                                                                                                                  <  2:742606/149‑1 (MQ=255)
                           aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGc                                                                                                                 >  2:118426/1‑149 (MQ=255)
                            aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACACGACGAAGCt                                                                                                                >  1:292192/1‑149 (MQ=255)
                                       cgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                                                                               >  2:985256/1‑123 (MQ=255)
                                       cgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAag                                                                                                                               <  1:985256/123‑1 (MQ=255)
                                         ttAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGCCGAAGCTTGAAGCAACGGTc                                                                                                   >  2:176996/2‑149 (MQ=255)
                                                tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCACAGc                                                                                            >  1:512679/1‑149 (MQ=255)
                                                      cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAAACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaca                                                                                      >  2:941340/1‑147 (MQ=255)
                                                       tCAACATTACGCTAACGGCACTAAAACCATCCCATTTTTCTGTGACTGGCGCTCCAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCCACGGTCGCACAGCTGaaaaa                                                                                     >  2:933269/1‑149 (MQ=255)
                                                       tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAACGCTGaaaaa                                                                                     >  2:936699/1‑149 (MQ=255)
                                                            aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGGGGTTCTTAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                                                                <  1:997500/149‑1 (MQ=255)
                                                              tacgtaaaCGGCACTAAAACGCTCACATCTTTTTGTGACTGGCGCTACAATCGTCCAAAGTCAAAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTgg                                                                              <  2:199207/143‑1 (MQ=255)
                                                                          aCTAAAACCATCACATTTTTCTGTTACTGGCGCTACAATTTTCCAAAGTCACAATTATCAAAATAGGAGTAGGGTTGCTAATGCCCCGAATGTTGTTAAATCGTGTAGCTACGGTCGTAATTCTGAAAAAGCATGGTGCTGAATCGGcc                                                                  >  2:717954/1‑149 (MQ=255)
                                                                          aCTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                                                                  >  1:415734/1‑149 (MQ=255)
                                                                              aaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgc                                                              >  1:533660/1‑149 (MQ=255)
                                                                                aCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                                                            >  2:517114/1‑149 (MQ=255)
                                                                                aCCATCACATTTTTCTGTGACTGGCGCTACAAGCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                                                            <  1:365393/149‑1 (MQ=255)
                                                                                       caTTTTTCTGGGACTGGCGCTACAAACTTCCAAAGCCACAATTCTCAAAATCAGAGGAGTATTGCGAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc                                                    <  1:189620/149‑1 (MQ=255)
                                                                                             tCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCgg                                              <  2:226323/149‑1 (MQ=255)
                                                                                              cTGTCACTGGCGCTACCATCTTCCAACGTCACAATTCTCAAACTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCGCGCGCATATCGCGCggggc                                           <  2:592908/149‑4 (MQ=255)
                                                                                                              aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGCAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCTATGCGCGCGGCATTCAGCGCCGAcccc                             >  1:78198/1‑149 (MQ=255)
                                                                                                                tCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccgg                           >  2:630711/1‑149 (MQ=255)
                                                                                                                tCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccgg                           >  1:561560/1‑149 (MQ=255)
                                                                                                                tCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccgg                           >  1:336290/1‑149 (MQ=255)
                                                                                                                  ttCCAAAGTCACAATTCTCAAAAACAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggcc                         <  1:118426/149‑1 (MQ=255)
                                                                                                                   tCCAAATTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACTGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccg                        <  1:81513/149‑1 (MQ=255)
                                                                                                                      aaaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGATAACGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCCTTCAGCGCCGACCACGCCCGGtc                     >  2:940090/1‑148 (MQ=255)
                                                                                                                            acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCg               >  1:405721/1‑149 (MQ=255)
                                                                                                                                   tcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGCCGAAGCTTGCAGCAACGGTCGAAAAGATGCAAAACCATTGGGACGACTCCACCAGGCCCGATATCGCCTCGCCATGGAGGGCCGACCCCCGCCGGTTCTGGCGATAttgc        >  2:593932/1‑147 (MQ=255)
                                                                                                                                   tcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCg        >  2:347624/1‑149 (MQ=255)
                                                                                                                                         tCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCGCTGc  >  1:433643/1‑149 (MQ=255)
                                                                                                                                         ccAGAGGAGTATTGCTAATGGAAAGAGACGTGACGAATCTATAAGAACCGGTCGCAAATCTGAAAAAGCCTTGGGCGGAATCGGCCCTGCGCGTTATCCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc  <  1:332068/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGC  >  NC_000913_3_bme_pgi/4233604‑4233889

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