Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,680 T→C intergenic (‑447/‑2227) 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,6800TC98.0% 161.2 / ‑6.8 50intergenic (‑447/‑2227)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (32/17);  minor base A (1/0);  total (33/17)
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.80e-01

TACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGG  >  NC_000913_3_bme_pgi/4233537‑4233824
                                                                                                                                               |                                                                                                                                                
taCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGCTTTACACTCAACATTACGCTCACGGCACCAAAAc                                                                                                                                             >  2:733702/1‑149 (MQ=255)
 aCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAAAATTACGCTAACGGCACCAAAAcc                                                                                                                                            >  2:714398/1‑149 (MQ=255)
   aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAAGAAATGCTTCACTGCGCTAAGGGTTTAAACTCAACATTCCGCTGACGGGACCGAAAcact                                                                                                                                          >  1:300268/1‑146 (MQ=255)
   aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATACACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACACAAACCAt                                                                                                                                          >  2:479738/1‑149 (MQ=255)
   aTTTAGAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAAACAt                                                                                                                                          <  2:377985/149‑1 (MQ=255)
     ttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCGAAACCATca                                                                                                                                        >  2:103315/1‑149 (MQ=255)
                gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCtgt                                                                                                                             >  1:367599/1‑149 (MQ=255)
                     cgCGGGGCGGTTTTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTg                                                                                                                        <  1:799148/148‑1 (MQ=255)
                             ggTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTGCa                                                                                                                >  1:255841/1‑149 (MQ=255)
                                  tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCtt                                                                                                           >  2:775416/1‑149 (MQ=255)
                                   cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                          <  2:537910/149‑1 (MQ=255)
                                       gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGAGCTAAGGGTTTATACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                      <  2:721904/149‑1 (MQ=255)
                                         tGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGt                                                                                                    >  2:62383/1‑149 (MQ=255)
                                           gggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTca                                                                                                  <  1:450412/149‑1 (MQ=255)
                                             gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAACGTcaca                                                                                                >  1:128341/1‑149 (MQ=255)
                                             gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTcaca                                                                                                >  1:131229/1‑149 (MQ=255)
                                                aTGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGATCGGGGTGTACACTCAACATTACGCTAACGGCACCAAAACCAGCACATTTTTCTGTGACTGGTGCTACAAGCTTCCAAAGTCACAAtg                                                                                             >  1:491768/1‑148 (MQ=255)
                                                   cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                          >  2:793215/1‑149 (MQ=255)
                                                       tttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTAACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTCCAATCTTGCAAAGTCACAATTCTCaaaa                                                                                      >  1:764409/1‑149 (MQ=255)
                                                        ttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAAt                                                                                     <  1:62383/149‑1 (MQ=255)
                                                          ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                   <  1:501090/149‑1 (MQ=255)
                                                              tCCACTACTGGCGGTCACTAGTTCAGCTAATAAATGCTTCACTTCGCTAAGGGTTTACACTCAACATGATGCTAACGGCACCAAAACCATCACGTTTTTCTGTGACTGGTGCTAGAATCTTCCAAAGTGACAAGTCTCAAAAGCagaag                                                                               >  2:686454/1‑149 (MQ=255)
                                                               ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                              <  1:219055/149‑1 (MQ=255)
                                                               ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                              <  1:392832/149‑1 (MQ=255)
                                                                 aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCAAAATTCTCAAAATCAGAAGAGt                                                                            >  1:472361/1‑149 (MQ=255)
                                                                     cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTg                                                                        >  2:145453/1‑149 (MQ=255)
                                                                      tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGAATTGc                                                                       >  2:527178/1‑149 (MQ=255)
                                                                        gCGGTCAGTAGTTCAGATAATAAATGCTTCACTGCGTTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTa                                                                     <  1:670667/149‑1 (MQ=255)
                                                                                     cAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcg                                                        >  2:563935/1‑149 (MQ=255)
                                                                                      cgCTCATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgc                                                       <  1:11790/148‑1 (MQ=255)
                                                                                        cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                     >  1:437988/1‑149 (MQ=255)
                                                                                                tGCTTCACTGCGATAAGGATTTACAATCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTCTAATCTTCCAAAGTCACAATTCTCAAAACTGGAGTAGTATTACTAATGGCAAGAGACGCGACGAAGCtt                                             <  2:142712/149‑1 (MQ=255)
                                                                                                   ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                          <  1:810203/149‑1 (MQ=255)
                                                                                                           gcTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAaga                                                            <  1:287976/123‑1 (MQ=255)
                                                                                                           gcTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAaga                                                            >  2:287976/1‑123 (MQ=255)
                                                                                                           gcTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGt                                  >  2:364191/1‑149 (MQ=255)
                                                                                                             tAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGtag                                >  1:336518/1‑147 (MQ=255)
                                                                                                                gggTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCACAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACAGTAGCaa                             >  2:232030/1‑149 (MQ=255)
                                                                                                                  gTTTACACTCAACATTACGCTAACGGCCCCCACTCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTACTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAg                           >  2:391846/1‑149 (MQ=255)
                                                                                                                  gTTTACACTCAACATGACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAg                           <  2:524610/149‑1 (MQ=255)
                                                                                                                   tttACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGACGCTTGAAGCAACGGTAGCaagcc                          >  2:205867/1‑146 (MQ=255)
                                                                                                                       cacTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGaa                      <  2:540508/149‑1 (MQ=255)
                                                                                                                          tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                     <  2:526336/115‑1 (MQ=255)
                                                                                                                          tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGcga                                                     >  1:526336/1‑115 (MQ=255)
                                                                                                                                    gCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCGTGAAGCAACGATAGCAAAGCTGAAAAAGCATTGGGCg         >  2:36610/1‑149 (MQ=255)
                                                                                                                                        aCGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAAt     >  1:82942/1‑149 (MQ=255)
                                                                                                                                        aCGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAAt     >  2:785882/1‑149 (MQ=255)
                                                                                                                                        aCGGCACCAAAACCATCACATTTTTCTGAGACTGGTGCCACCATCGTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGTCGAAGCTTGAAGCAACGGTAGCATAGCTGAACCAGCATTGGGCGGagt     >  1:32544/1‑147 (MQ=255)
                                                                                                                                          ggCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCCACGCAGCTTGAAGCAACGATCGCAACGCTGAAAAAGCATGTGCCGGAATCg   >  2:227051/1‑149 (MQ=255)
                                                                                                                                           gCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCgg  >  1:149871/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGG  >  NC_000913_3_bme_pgi/4233537‑4233824

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