Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,920,598	T→C	E97E (GAA→GAG)	atpF ←	ATP synthase Fo complex subunit b

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,920,598	0	T	C	100.0%	83.7 / NA	25	E97E (GAA→GAG)	atpF	ATP synthase Fo complex subunit b
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/11);  total (14/11)

GAACGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCA  >  NC_000913/3920473‑3920727                                                                                                                              |                                                                                                                                  gAACGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCACCTGGGCAACGCTTGTAGTACGCTCCTGtt                                                                                                                            >  1:233704/1‑133 (MQ=255)                  ctcGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCt                                                                                                           <  2:509110/133‑1 (MQ=255)                       cgcCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACTCTCCTGTACTGCCTACGCTTTAGCTTCGTc                                                                                                      >  1:201639/1‑133 (MQ=255)                         ccAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCa                                                                                                    >  1:368378/1‑133 (MQ=255)                            gCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCtt                                                                                                          >  1:154685/1‑124 (MQ=255)                            gCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCtt                                                                                                          <  2:154685/124‑1 (MQ=255)                             cAACAGCCAGGATAGCAACTTGCTTACGC‑GCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATc                                                                                               <  2:550137/133‑1 (MQ=255)                                    cAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAg                                                                                         >  1:332474/1‑133 (MQ=255)                                       gATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAgcgg                                                                                      >  1:349068/1‑133 (MQ=255)                                        aTAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAgcggc                                                                                     <  2:408607/133‑1 (MQ=255)                                                       cGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCt                                                                      <  1:268296/133‑1 (MQ=255)                                                         cAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTc                                                                     <  2:502833/132‑1 (MQ=255)                                                         cAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTc                                                                     >  1:502833/1‑132 (MQ=255)                                                            ctctTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCgatg                                                                 >  1:458124/1‑133 (MQ=255)                                                                  aCGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTAcc                                                           <  2:565209/133‑1 (MQ=255)                                                                         cGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCtt                                                    >  1:468773/1‑133 (MQ=255)                                                                                    ggCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgctttcgct                                         >  1:621198/1‑133 (MQ=255)                                                                                      cTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCcgctttcgcttt                                       >  1:517533/1‑133 (MQ=255)                                                                                              tCCGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCTCTTTCGCTTTTTTCAGCt                               >  1:232561/1‑133 (MQ=255)                                                                                               ccGCCTGCGCCTGGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTg                              <  2:621198/133‑1 (MQ=255)                                                                                                 gcctgcgcctgGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggt                            <  2:517533/133‑1 (MQ=255)                                                                                                      cgcctgGGCCACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTAGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTggtcggtc                       >  1:558740/1‑133 (MQ=255)                                                                                                          tgGGCCACGATTTTAGTACGCTCCCGTTCTGCATCCGCTTTGGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCgcgc                   <  2:232561/133‑1 (MQ=255)                                                                                                              ccACGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGc               <  2:332474/133‑1 (MQ=255)                                                                                                                aCGATTTTAGTACGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCt             >  1:11240/1‑133 (MQ=255)                                                                                                                           aCGCTCCTGTTCTGCCTCCGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGcaaggtca  <  2:458124/133‑1 (MQ=255)                                                                                                                              |                                                                                                                                  GAACGTTCGATGATCTTCTCGGCGCCAGCAACAGCCAGGATAGCAACTTGCTTACGCAGCTCTTCACGGGCACGTTTACGCTCGGCTTCAATTTCCGCCTGCGCCTGGGCCACGATTTTAGTACGTTCCTGTTCTGCCTCAGCTTTCGCTTCGTCCAGAATCTGCGAGCGGCGTTTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCA  >  NC_000913/3920473‑3920727

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