Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,477,281	Δ1 bp	intergenic (‑391/‑29)	b4256 ← / → yjgN	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,477,279	0	T	.	100.0%	151.8 / NA	34	intergenic (‑389/‑31)	b4256/yjgN	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (21/13);  total (21/13)

TACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGT  >  NC_000913/4477204‑4477351                                                                            |                                                                         tACTGATAAGAGTATTGTCACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAt                                                                          >  2:2203571/1‑76 (MQ=255) tACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAt                                                                          >  1:3228304/1‑76 (MQ=255)            gTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAttt                                                              <  1:1284316/76‑1 (MQ=255)             tATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAtttt                                                             <  1:3508686/76‑1 (MQ=255)              aTTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAttttt                                                            >  2:1823507/1‑76 (MQ=255)                tGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTAt                                                          >  1:130956/1‑76 (MQ=255)                 gTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATg                                                         <  1:2119300/76‑1 (MQ=255)                 gTAACGTGGCAGAACATGTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TGCAGGGAATTGTTATg                                                         <  1:2127442/76‑1 (MQ=255)                  tAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATgg                                                        >  1:1969849/1‑76 (MQ=255)                   aaCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGc                                                       <  2:1747762/76‑1 (MQ=255)                      gTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TCAAGGGAATTTTTATGCCTCa                                                    >  2:3855308/1‑76 (MQ=255)                       tGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑CTCAGGGAATTTTTATGGCTCaa                                                   <  1:800703/76‑1 (MQ=255)                         gCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGt                                                 >  1:1895267/1‑76 (MQ=255)                           aGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGtta                                               >  1:4043079/1‑76 (MQ=255)                            gAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGttat                                              >  2:3279002/1‑76 (MQ=255)                                     aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAAt                                     >  2:4286276/1‑76 (MQ=255)                                     aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAAt                                     >  1:220730/1‑76 (MQ=255)                                       agtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATgg                                   >  2:2320240/2‑76 (MQ=255)                                          ctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtg                                <  1:4521844/76‑2 (MQ=255)                                           tgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtgt                               >  2:2587992/1‑74 (MQ=255)                                               tatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCg                          <  2:602583/76‑1 (MQ=255)                                                    aaGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCa                     >  2:2431225/1‑76 (MQ=255)                                                     aGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCAt                    <  1:709517/76‑1 (MQ=255)                                                     aGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCAt                    <  2:3989531/76‑1 (MQ=255)                                                       ttATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTc                  >  2:1134517/1‑76 (MQ=255)                                                       ttATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTc                  >  1:4357762/1‑76 (MQ=255)                                                        tataTCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                 >  1:1888757/1‑76 (MQ=255)                                                        tataTCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                 >  1:1867464/1‑76 (MQ=255)                                                                tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt         <  1:2719885/76‑1 (MQ=255)                                                                    ttttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCAt     >  2:984744/1‑76 (MQ=255)                                                                     tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg    >  1:3662328/1‑76 (MQ=255)                                                                     tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg    <  1:2171652/76‑1 (MQ=255)                                                                      ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATgg   >  1:1240482/1‑76 (MQ=255)                                                                      ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATgg   >  1:1139031/1‑76 (MQ=255)                                                                       tttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt  <  2:2480564/76‑1 (MQ=255)                                                                       tttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt  >  2:2154756/1‑76 (MQ=255)                                                                            |                                                                         TACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGT  >  NC_000913/4477204‑4477351

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