Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,556,610	T→C	intergenic (‑290/‑383)	yjiC ← / → iraD	uncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,556,610	0	T	C	100.0%	107.8 / NA	29	intergenic (‑290/‑383)	yjiC/iraD	uncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/13);  total (16/13)

CAAAAAACGCTCATTTAAAATTATTTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGTTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTC  >  NC_000913/4556477‑4556745                                                                                                                                      |                                                                                                                                        cAAAAAACGCTCATTTAAAATTATTTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCa                                                                                                                                    >  1:1137181/1‑139 (MQ=255) cAAAAAACGCTCATTTAAAATTATTTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCa                                                                                                                                    >  1:1039753/1‑139 (MQ=255)               ttAAAATTATTTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTa                                                                                                                      <  1:1205013/139‑1 (MQ=255)                       attTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTaa                                                                                                              <  2:1137181/139‑1 (MQ=255)                       attTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTAAATTATTTACAAACTaa                                                                                                              <  2:1039753/139‑1 (MQ=255)                                            atatCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGc                                                                                         >  2:736696/1‑139 (MQ=255)                                                                                   aGGTCTATATTAATAATGCATGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAg                                                  >  2:878267/1‑139 (MQ=255)                                                                                        tataTTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTagca                                             >  1:1288473/1‑139 (MQ=255)                                                                                          tatTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTagcagc                                           >  2:152975/1‑139 (MQ=255)                                                                                             taataaTGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGt                                        <  1:318872/139‑1 (MQ=255)                                                                                               ataatGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGa                                      >  1:212914/1‑139 (MQ=255)                                                                                               ataatGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGa                                      >  2:603473/1‑139 (MQ=255)                                                                                                taatGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTag                                               <  1:248521/129‑1 (MQ=255)                                                                                                taatGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTag                                               >  2:248521/1‑129 (MQ=255)                                                                                                taatGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGAc                                     >  2:32230/1‑139 (MQ=255)                                                                                                taatGCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGAc                                     <  2:228697/139‑1 (MQ=255)                                                                                                    gCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGt                                 >  2:18515/1‑139 (MQ=255)                                                                                                    gCCTGTGAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGt                                 >  1:696081/1‑139 (MQ=255)                                                                                                         tgAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGAc                                     <  2:690842/130‑1 (MQ=255)                                                                                                         tgAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGAc                                     >  1:690842/1‑130 (MQ=255)                                                                                                          gAATGGTATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTc                           <  1:152975/139‑1 (MQ=255)                                                                                                                tATTTTTGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTAc                     >  1:1310615/1‑139 (MQ=255)                                                                                                                   ttttGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTa                  <  2:696081/139‑1 (MQ=255)                                                                                                                    tttGATGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTaa                 <  2:413028/139‑1 (MQ=255)                                                                                                                        aTGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTAt             <  1:760515/139‑1 (MQ=255)                                                                                                                         tGTATTTGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTAtt            >  2:796033/1‑139 (MQ=255)                                                                                                                             tttGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGa        <  1:736696/139‑1 (MQ=255)                                                                                                                             tttGATATGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGa        <  1:18515/139‑1 (MQ=255)                                                                                                                                   atGCTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGGAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTc  >  2:786584/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        CAAAAAACGCTCATTTAAAATTATTTGCATGCAATTTAAAAGCATATCTTATTACTAATTGGAATTTGATGTTGCTATATTGAGGTCTATATTAATAATGCCTGTGAATGGTATTTTTGATGTATTTGATATGTTATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTC  >  NC_000913/4556477‑4556745

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