Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,887,571	Δ1 bp	intergenic (+163/+7)	sokX → / ← cysH	small RNA SokX/phosphoadenosine phosphosulfate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,887,570	0	G	.	100.0%	109.3 / NA	26	intergenic (+162/+8)	sokX/cysH	small RNA SokX/phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/11);  total (15/11)

TTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCC  >  NC_000913/2887498‑2887643                                                                         |                                                                          tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATAtt                                                                   >  2:469050/1‑80 (MQ=255)    aaaaTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTAc                                                                <  2:469020/80‑1 (MQ=255)       aTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCtt                                                            <  2:469021/81‑1 (MQ=255)        tCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTc                                                           <  2:469022/81‑1 (MQ=255)        tCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTc                                                           >  2:469047/1‑81 (MQ=255)           ggAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCgtg                                                        <  1:469023/81‑1 (MQ=255)              aTGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTaa                                                     >  2:469052/1‑81 (MQ=255)                gctgctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAAcc                                                   >  1:469045/1‑81 (MQ=255)                    ctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCca                                                 <  2:469024/79‑1 (MQ=255)                        cAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTcc                                           >  1:469044/1‑81 (MQ=255)                         aaaaTAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTccc                                          >  2:469054/1‑81 (MQ=255)                              aaaaaCGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCtttt                                      >  2:469059/1‑80 (MQ=255)                                  aCGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc                                 <  1:469025/81‑1 (MQ=255)                                   cGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc                                 >  2:469053/1‑80 (MQ=255)                                     cccccATCATTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaa                              <  2:469026/81‑1 (MQ=255)                                         catcatTTCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaa                          <  1:469027/81‑1 (MQ=255)                                               ttCTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGt                    >  2:469060/1‑81 (MQ=255)                                                 cTGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGttt                  >  2:469055/1‑81 (MQ=255)                                                  tGACAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTttc                 >  1:469062/1‑81 (MQ=255)                                                    aCAGAGGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGttt                  <  1:469028/78‑1 (MQ=255)                                                       gagGCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTttcttctt            <  2:469029/81‑1 (MQ=255)                                                          gCGTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCg         >  2:469058/1‑81 (MQ=255)                                                            gTTTAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGcc       <  2:469030/81‑1 (MQ=255)                                                              ttAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCAt     >  1:469057/1‑81 (MQ=255)                                                              ttAATTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCAt     >  2:469069/1‑81 (MQ=255)                                                                 aTTTGTCC‑GCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGcc  >  1:469068/1‑81 (MQ=255)                                                                         |                                                                          TTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCC  >  NC_000913/2887498‑2887643

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