Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,041,671	(T)6→7	intergenic (+46/‑108)	dsbA → / → b3861	protein disulfide isomerase I/putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,041,665	1	.	T	100.0%	13.5 / NA	5	intergenic (+40/‑114)	dsbA/b3861	protein disulfide isomerase I/putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (1/4);  total (1/4)

TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCT  >  NC_000913/4041573‑4041711                                                                                              |                                               tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                 <  2:2790914/93‑1 (MQ=255) tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                 <  2:2344911/93‑1 (MQ=255) tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                 <  2:3477049/93‑1 (MQ=255)  ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑t                                               <  1:1139234/93‑1 (MQ=255)  ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑t                                               <  2:5921922/93‑1 (MQ=255)  ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑t                                               <  2:1228696/93‑1 (MQ=255)   tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tt                                              <  2:3164806/93‑1 (MQ=255)   tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tt                                              <  2:1332199/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                             <  1:5900843/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                             <  1:6306427/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                             <  1:5655722/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                             <  1:190703/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                             <  1:1369558/93‑1 (MQ=255)     tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                            <  1:925652/93‑1 (MQ=255)     tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                            <  1:2654235/93‑1 (MQ=255)      gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                           <  2:3700138/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                          <  2:960259/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                          <  2:4862191/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                          <  2:4961724/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                          <  1:5302254/93‑1 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                        <  2:5607268/93‑3 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                       <  2:2030746/93‑4 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                       <  2:5858096/93‑4 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                       <  1:2592198/93‑4 (MQ=255)           gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                      <  1:5876229/93‑5 (MQ=255)             agTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                     <  2:405653/93‑1 (MQ=255)              gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                    <  2:3961382/93‑1 (MQ=255)                aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTta                                  <  2:2534922/93‑1 (MQ=255)                 tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                 <  1:4865164/93‑1 (MQ=255)                                                aaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCt  >  2:1520196/1‑93 (MQ=255)                                                                                              |                                               TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCT  >  NC_000913/4041573‑4041711

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