Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,793 C→A intergenic (‑560/‑2114) 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,7930CA98.1% 150.0 / NA 53intergenic (‑560/‑2114)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (26/26);  total (27/26)
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

GGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGA  >  NC_000913_3_bme_pgi/4233649‑4233936
                                                                                                                                                |                                                                                                                                               
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)
                    gCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCGTGAAGCAACGATAGCAAAGCTGAAAAAGCATTGGGCg                                                                                                                         >  2:36610/1‑149 (MQ=255)
                        aCGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                                     >  2:785882/1‑149 (MQ=255)
                        aCGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                                     >  1:82942/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)
                             accaaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCCCAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                                                                                                                <  1:224308/146‑1 (MQ=255)
                             accaaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGACGAGTATTGCTAATGGCAACAGACGCGACGAAGCTTGAACCAACGGTAGCACACCCGAAAACGCATTGGCCGGAATCGGcc                                                                                                                >  2:369109/4‑149 (MQ=255)
                               caaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCAAAATTGTCTAAATCAGAAGAGTCTTGCTAATGGCACGAGACGCGACGAAGCTTGAAGCAACGGTAGCAGAGCTGAAAAAGCCTTGGTCGGAATCGGgccc                                                                                                              >  2:151183/2‑149 (MQ=255)
                                            tttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGc                                                                                                                          <  1:16907/124‑1 (MQ=255)
                                            tttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGc                                                                                                                          >  2:16907/1‑124 (MQ=255)
                                              tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGACCATTAATGCGACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc                                                                                               >  1:405241/1‑149 (MQ=255)
                                                       cTGGTGCTACAATCTTCCAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAg                                                                                      <  1:733702/149‑1 (MQ=255)
                                                       cTGGTGCCACAATCTTCCAAAGTCACAATTTTCAAAATAAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAg                                                                                      <  2:469033/149‑1 (MQ=255)
                                                        tGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc                                                                                     >  2:691302/1‑149 (MQ=255)
                                                                gaaTCTTCCAAAGTCACAATTCTCAACATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTCGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAccc                                                                             <  2:169434/148‑1 (MQ=255)
                                                                cAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAccc                                                                             >  1:559620/1‑149 (MQ=255)
                                                                 aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGACAAAGCATTGGGCGCAATCGGCCCCGCGCGATATGCGCGCGGCCTTCAGCGCCGAcccc                                                                            >  2:434281/1‑149 (MQ=255)
                                                                 aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAcccc                                                                            >  2:781184/1‑149 (MQ=255)
                                                                    gttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggc                                                                         <  2:547276/148‑1 (MQ=255)
                                                                    cTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGTCATTCAGCGCCGCACccggc                                                                         >  2:25344/1‑149 (MQ=255)
                                                                         aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAAAGCCGACCCCGGCCGGtt                                                                    >  2:758278/1‑149 (MQ=255)
                                                                          aaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                                                   <  2:559620/149‑1 (MQ=255)
                                                                             tCAAAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                                                                <  1:453637/149‑1 (MQ=255)
                                                                              cacaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                                                               <  1:103315/149‑1 (MQ=255)
                                                                               acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCTGCGCCGACCCCCGCCGGTTCGGGCg                                                              >  2:462885/1‑149 (MQ=255)
                                                                                     ctcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTc                                                        <  2:437988/149‑1 (MQ=255)
                                                                                                aagaGTATTGCTCATGGGGATAGGCACGACGAAGCTTGCAGAGCCGGCAGCAAAGCCAAAAATTCATCGTATGGTAGCTGCCCCGCGCGATATGCGGGCGGCATTCAGCGCCGACCCCGGCGGGTTCGGGCGATCTTAGCTCTCCTTgg                                             >  2:177522/1‑149 (MQ=255)
                                                                                                 agagTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGga                                            <  1:714398/149‑1 (MQ=255)
                                                                                                     tATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgatgat                                        >  1:512447/1‑149 (MQ=255)
                                                                                                      aTTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtc                                       <  2:82942/149‑1 (MQ=255)
                                                                                                          cTAATGGCAAGAGACGGGACGAAGGTTGAAGCAACGGTAGCAAATCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGTGATATTCGCTCTGCCTGGATGAtcttct                                   <  2:255841/149‑1 (MQ=255)
                                                                                                          cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTATCACTGGATGAtcttct                                   >  2:550048/1‑149 (MQ=255)
                                                                                                          cTAAGGGCAAGAGCCGCGACGAGGCTTGAAGGAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtcttct                                   <  2:304915/149‑1 (MQ=255)
                                                                                                            aaTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATAGGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt                                 <  2:322925/149‑1 (MQ=255)
                                                                                                             aTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGtt                                <  1:391846/149‑1 (MQ=255)
                                                                                                                   agagACGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCAGTGGGCCAAAACGGCCCCGCGCGAGATGCGTGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTg                          <  2:719076/149‑1 (MQ=255)
                                                                                                                    gagaCGCGACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTgg                         <  1:691302/149‑1 (MQ=255)
                                                                                                                     aggcgcgAGGAAGCGTGAAGCAACGGTAGGAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGt                        <  2:367599/146‑1 (MQ=255)
                                                                                                                         gcgACGAAGCTTGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTAAACGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa                    <  1:785882/149‑1 (MQ=255)
                                                                                                                                   ttGAAGCAACGGTAGCAACGCTGAAACGCATTGGGCGGACTCCGCCCCGAGCGATATCCGCGCGCCCTTCAGCGCCGACCCCGGCCGGTGCGGGCGATATTCGCCCTGCCCGGGTTCTCTTCTGTTCGCCTGGTCGAAGTGcccgctca         >  2:377198/1‑143 (MQ=255)
                                                                                                                                   ttGAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGGCGCAGAGCCCGGta          >  2:274620/1‑148 (MQ=255)
                                                                                                                                     gAAGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACGGGTAGAAGTGCCggggtga        >  1:585501/1‑145 (MQ=255)
                                                                                                                                      aaGCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAac       <  1:781184/149‑1 (MQ=255)
                                                                                                                                        gCAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAacga     <  1:364191/149‑1 (MQ=255)
                                                                                                                                         cAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAacgac    >  2:305205/1‑149 (MQ=255)
                                                                                                                                           aCGGTAGGAAGGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga  <  1:145453/149‑1 (MQ=255)
                                                                                                                                           aCGGTAGAAAAGCGTAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga  <  1:527178/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGA  >  NC_000913_3_bme_pgi/4233649‑4233936

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