Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,608,076	A→G	T92A (ACT→GCT)	holD →	DNA polymerase III subunit psi

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,608,076	0	A	G	100.0%	56.1 / NA	17	T92A (ACT→GCT)	holD	DNA polymerase III subunit psi
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/7);  total (10/7)

TGATCCTTTAGTGAGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCC  >  NC_000913/4607952‑4608203                                                                                                                             |                                                                                                                                tgaTCCTTTACTGAGCGATGTTCTTCGCGCATTAACCGTCACCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGaa                                                                                                                         <  2:184973/133‑1 (MQ=255)    tCCTTTAGTGAGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCg                                                                                                                      >  2:1939/1‑133 (MQ=255)     ccTTTAGTGAGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGc                                                                                                                     >  1:231182/1‑133 (MQ=255)         tAGTGAGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATc                                                                                                                 >  2:524111/1‑133 (MQ=255)             gAGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTg                                                                                                             >  1:57948/1‑133 (MQ=255)              aGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTgg                                                                                                            >  1:308651/1‑133 (MQ=255)                cGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGaa                                                                                                          >  1:114573/1‑133 (MQ=255)                 gATGTTCTGCTCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAg                                                                                                         <  2:231182/133‑1 (MQ=255)                                         aGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGgcgc                                                                                 <  1:1939/133‑1 (MQ=255)                                                                    cgaGAAAAACTCGCGGTGCTGCCTTACGCCAGGCACTGCAAGAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACgg                                                      <  2:300871/131‑1 (MQ=255)                                                                      aGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCa                                                    >  2:533462/1‑133 (MQ=255)                                                                        aaaaaaTCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCAcg                                                  <  2:57948/133‑1 (MQ=255)                                                                          aaaaTCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCAcgcg                                                <  1:524111/133‑1 (MQ=255)                                                                                                        tGCAACAGTTGGCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATg                  >  1:534020/1‑133 (MQ=255)                                                                                                            aCAGTTGGCGGTTGGGTGCTGACGCACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAaca              <  2:234382/133‑1 (MQ=255)                                                                                                                  ggCGGTTGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGAttt        >  1:347705/1‑133 (MQ=255)                                                                                                                        tGGGTGCTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTccc  >  2:62001/1‑133 (MQ=255)                                                                                                                             |                                                                                                                                TGATCCTTTAGTGAGCGATGTTCTGCGCGCATTAACCGTCAGCCCCGACCAGGTGCTGCAACTGACGCCAGAAAAAATCGCGATGCTGCCGCAAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCC  >  NC_000913/4607952‑4608203

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