Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,564,097 (A)7→6 intergenic (‑92/+37) yzgL ← / ← glgP putative uncharacterized protein YzgL/glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,564,0910A.76.2% 33.3 / 1.8 22intergenic (‑86/+43)yzgL/glgPputative uncharacterized protein YzgL/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base . (5/11);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

CAGTATGGGGATGCCGAAAAAAGTCATTACCAGCGAGGTCAAAATCCATTTTCTGTTTTGCATTATTCTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAC  >  NC_000913/3563943‑3564450
                                                                                                                                                         |                                                                                                                                                                                                                                                                                                                                                                            
agtaaTGGGGATGCCGAAAAAAGTCATTACCAGCGAGGTCAAAATCCATTTTCTGTTTTGCATTATTCTTTCCATTCTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGa                                                                                                                                                                                                                                                                                                                 <  1:1465168/201‑1 (MQ=255)
                                 gaggttCAAAATCCATTTTCTGTTTTGCATTATTCTTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑‑‑AAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCttt                                                                                                                                                                                                                                                                                <  1:1465171/200‑1 (MQ=255)
                                                  ttCTGTTTTGCATTATTCTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTATCCTGACATCGACATAGCTg                                                                                                                                                         >  1:1465174/1‑305 (MQ=255)
                                                    cTGTTTTGCATTATTCTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCG‑‑‑‑AAAAACGGGACCCTTTGGCCCCGTTCTATTTATTTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTagat                                                                                                                                                                                           <  1:1465172/269‑2 (MQ=255)
                                                               tattCTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑‑‑‑AAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGa                                                                                                                                                                        <  1:1465173/276‑1 (MQ=255)
                                                                       tCCATTCTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCtt                                                                                                                                                                                                                                                                                                                                                             <  1:566142/94‑1 (MQ=255)
                                                                            ttcttttGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTCCGAATGATGGCGG‑‑AAAAAAACGGGACCCTTTGGCCCCGTTCtattt                                                                                                                                                                                                                                                                                                                                             <  1:566141/101‑1 (MQ=255)
                                                                             cTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCACTCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcgc                                                                                            >  1:566144/1‑341 (MQ=255)
                                                                             cTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGc                                                                                                                                                                                                                                                                                                        >  1:1465175/1‑138 (MQ=255)
                                                                                       ggTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGACG‑‑GAAAAAACGGGACCCTTTGGCCCCGTTCtaatt                                                                                                                                                                                                                                                                                                                                             <  1:566143/95‑1 (MQ=255)
                                                                                                      cccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑‑‑AAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCg                                                                                                                        <  1:1465176/287‑1 (MQ=255)
                                                                                                             ccccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑‑‑AAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCagaagag                                                                                                                                                                                                                                                             <  1:1465177/148‑1 (MQ=255)
                                                                                                              cccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAAACGGGACCCTTTGGCCCCGt                                                                                                                                                                                                                                                                                                                                                    >  1:566146/1‑65 (MQ=255)
                                                                                                                  tGCCCCTTCTGGACACTTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCa                                                                                                      >  1:1465182/1‑293 (MQ=255)
                                                                                                                            tGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCAt                                                                                                                                                                                                                               >  1:1465184/1‑163 (MQ=255)
                                                                                                                            tGGACACTTTTCCGAAATGAT‑‑‑GG‑‑AGAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATc                                                                                                       <  1:1465178/282‑1 (MQ=255)
                                                                                                                             ggACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTg                                                                                                                                                                                                                       <  1:1465179/171‑1 (MQ=255)
                                                                                                                              gACACTTTTCCGAAATGATGGCGG‑‑AAAAAACGGGACCCTTTGGCCCCGTTCTATTTAATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGc                                                                                                                                                                                                                                                      <  1:1465180/141‑1 (MQ=255)
                                                                                                                                    ttttCCGAAATGATGGCGG‑‑‑‑AAAAACGGGACCCTTTGGCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTg                                                                                                                                                                         <  1:1465181/210‑1 (MQ=255)
                                                                                                                                     tttCCGAAATGATGGCGGAAAAAAAAACGGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc  >  1:566148/1‑379 (MQ=255)
                                                                                                                                     tttCCGAAATGATGGCGGAAAAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACttt                                                                                                                                                                               >  1:566147/1‑205 (MQ=255)
                                                                                                                                             atgatgGCGG‑‑‑‑AAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCg                                                                                                                                                                                                <  1:1465183/178‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                                                                                                                                                                                                                                            
CAGTATGGGGATGCCGAAAAAAGTCATTACCAGCGAGGTCAAAATCCATTTTCTGTTTTGCATTATTCTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGG‑‑AAAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAC  >  NC_000913/3563943‑3564450

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