Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,057,133	(A)7→8	intergenic (+142/+45)	sibC → / ← serA	small RNA SibC/phosphoglycerate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,057,126	1	.	A	100.0%	62.9 / NA	18	intergenic (+135/+52)	sibC/serA	small RNA SibC/phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (9/9);  total (9/9)

CTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑AAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAT  >  NC_000913/3057053‑3057201                                                                           |                                                                            cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑aaaaaa                                                                       >  1:503121/1‑80 (MQ=255)   gACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑aaaaaaaac                                                                    >  2:503116/1‑79 (MQ=255)    aCTTTGATTTATTTTCTGGAACGGACTCGAAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑aaaaa                                                                        <  1:503087/76‑1 (MQ=255)     cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑aaaaaaaacgg                                                                  >  2:503111/1‑77 (MQ=255)        tGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGc                                                                 <  2:590395/80‑1 (MQ=255)               ttttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGt                                                         <  1:503089/81‑1 (MQ=255)                tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGt                                                         >  1:503119/1‑80 (MQ=255)                      gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGAcc                                                     <  2:503089/78‑1 (MQ=255)                      gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCt                                                    <  1:503091/79‑1 (MQ=255)                         cAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                               >  1:503126/1‑81 (MQ=255)                          aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGt                                              <  1:503092/81‑1 (MQ=255)                           gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGtt                                             <  1:503093/81‑1 (MQ=255)                              tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt                                          >  1:503123/1‑81 (MQ=255)                                  aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                       >  1:503122/1‑80 (MQ=255)                                       aGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag                                  >  2:503127/1‑80 (MQ=255)                                                aaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                        <  1:503094/81‑1 (MQ=255)                                                aaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                         >  1:503134/1‑80 (MQ=255)                                                           ggTGTGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGa             >  2:503136/1‑81 (MQ=255)                                                            gtgtGGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc            >  2:503124/1‑81 (MQ=255)                                                               tgGAGAAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACggg         >  2:503126/1‑81 (MQ=255)                                                                 gagaAGGGATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg       <  1:503095/81‑1 (MQ=255)                                                                      gggATAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt  <  1:503096/81‑1 (MQ=255)                                                                           |                                                                            CTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑AAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAT  >  NC_000913/3057053‑3057201

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