Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,275,204	C→T	intergenic (‑163/‑267)	yjcB ← / → pdeC	uncharacterized protein YjcB/c‑di‑GMP phosphodiesterase PdeC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,275,204	0	C	T	96.7%	81.8 / NA	30	intergenic (‑163/‑267)	yjcB/pdeC	uncharacterized protein YjcB/c‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/14);  minor base A (1/0);  total (16/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 = 9.80e-01

TGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCC  >  NC_000913/4274957‑4275441                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              tGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCat                                                                                                                                                                                                                                            >  1:436900/1‑251 (MQ=255)     gAGTGGGCTGGCCATAAACATCATTTCGGCACTAAGAAGTTGACAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACTCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATAtt                                                                                                                                                                                                                                         <  2:358618/250‑1 (MQ=255)         gtgCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGcc                                                                                                                                                                                                                                    <  1:624312/251‑1 (MQ=255)              ggCCAGAAACATCATGACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGAc                                                                                                                                                                                                                               <  2:33687/251‑1 (MQ=255)                ccAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACCCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAAAACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGcc                                                                                                                                                                                                                                    >  2:2626/1‑244 (MQ=255)                ccAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGcc                                                                                                                                                                                                                                    <  1:2626/244‑1 (MQ=255)                                              gCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCGGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                                                                                                                                                                                <  2:593132/250‑1 (MQ=255)                                                       gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                                                                                                                                                                                       >  1:317739/1‑250 (MQ=255)                                                                              gTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGAc                                                                                                                                                                                                                               >  1:363263/1‑187 (MQ=255)                                                                              gTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTGGTTCACGGGAGGCGTGGGACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGTCCGTGAc                                                                                                                                                                                                                               <  2:363263/187‑1 (MQ=255)                                                                                agcctacAAACAATTCCACTGAAAAGATCCTATGCCCTGGACAGAGATGCGTTGTTCACGGGTGGCGTGGTTCCCTCTGTCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTTTTTGTTTCATTAATTTTGTGAATTTTATCACAATTGTTTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTAAGTGGCATTAAGGGCATACCTTCCCTCCGCctct                                                                                                                                                              <  2:573512/244‑1 (MQ=255)                                                                                   cATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAt                                                                                                                                                          >  1:231440/1‑251 (MQ=255)                                                                                       aaaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCt                                                                                                                                                      >  2:635442/1‑251 (MQ=255)                                                                                        aaccaTTCCTCTGAAAAGGGCCGATTCACGGGACTGCGCGGCTTTGTTCACGGGAGGCGTGTTACACTCTGGCTATCGCGGGGATTGCAGAACACAAAACTGAAACCCTCTGTTTGTTTCATTAATTTTGTGAATTATATCTCAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTTGATTTTCAGTGGCATTAAGGGCAGACCTTCCCTCCGACCCTGCATTCCTg                                                                                                                                                     <  2:449978/247‑1 (MQ=255)                                                                                         aCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGt                                                                                                                                                    >  1:430071/1‑251 (MQ=255)                                                                                            aTTCCTCTGAAAAGCGCCGATGCCCTGGCCAGCGCGGCGTTGTTCACGGGCGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAgg                                                                                                                                                 >  2:105383/1‑251 (MQ=255)                                                                                                                       gACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGTCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAatt                                                                                                                       >  1:473058/1‑250 (MQ=255)                                                                                                                          aGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTAttt                                                                                                                   >  1:74788/1‑251 (MQ=255)                                                                                                                                               ggCGTGGTACACTCTGGCTATCTCGGGGCTTTCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTTTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCATGTGACGATTAGaa                                                                                              <  1:662955/251‑1 (MQ=255)                                                                                                                                                                  agcgcgGGGCTTGAAGAACACAAAAATGAAACTCTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCa                                                                           <  2:637624/249‑1 (MQ=255)                                                                                                                                                                      cgGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGc                                                                        >  1:177113/1‑250 (MQ=255)                                                                                                                                                                                    cacaAAAATGAAACACCCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTACCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTtat                                                         <  2:156366/251‑1 (MQ=255)                                                                                                                                                                                                acacTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCatat                                                                                                                                                                             >  1:522477/1‑123 (MQ=255)                                                                                                                                                                                                acacTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCatat                                                                                                                                                                             <  2:522477/123‑1 (MQ=255)                                                                                                                                                                                                      tgtatgtttCATTAATTTTGTGAACTATATCATAATTGATTGTTTTTGTGTCATATTAGGCCATGACTTTTATTTCTGTACATATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTa                                        <  2:556341/246‑1 (MQ=255)                                                                                                                                                                                                       gtttgttgCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGATTTTTAGTGCTGTACAGATTATGTGGTTTGTCATTGGCATTAAGGGCATTTCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTGCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTaaa                                      <  1:388260/251‑1 (MQ=255)                                                                                                                                                                                                               cATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCgg                              >  2:157435/1‑251 (MQ=255)                                                                                                                                                                                                                                 tatCACAATTAATTATTTGTTAGACATGTTAACCCGTGCCTTTTATTGCTGTACAGATTATGAGGTTTTTCAGTGGCATAAAGGGCATATCTTCCAGCCGCCTCTGCATTCCTGTAGGCTATTAATTTTGCATATCAATGAATTATTTTCATCCAGGAGACCATTAGAAAGGTATAAATTTCAAATAAGGCAACAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAAGCTTTTCACCt            >  2:616816/1‑251 (MQ=255)                                                                                                                                                                                                                                           gattgTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGcc  >  1:279721/1‑251 (MQ=255)                                                                                                                                                                                                                                           gattgTTTGTTAGTCATATTAGGCCGTGACTTTTATTGCTGTACAGATGATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGCATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCTAGGGTTTTCACCTTGCAATGGcc  >  1:369422/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              TGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCC  >  NC_000913/4274957‑4275441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑