Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,664,096	T→C	I163T (ATA→ACA)	ynfH →	oxidoreductase, membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,664,096	0	T	C	100.0%	135.1 / NA	38	I163T (ATA→ACA)	ynfH	oxidoreductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/18);  total (20/18)

AGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAATACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTG  >  NC_000913/1663969‑1664197                                                                                                                                |                                                                                                      aGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCg                                                                                            >  1:820083/1‑139 (MQ=255)                      ttCGTCTGGATGATGACCTGCGTGTATCAGATCGACACCGTGCCATCCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTtgcg                                                                       <  2:731250/139‑1 (MQ=255)                        cGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                     <  1:875737/139‑1 (MQ=255)                           cTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACa                                                                  <  2:926235/139‑1 (MQ=255)                              gatgatGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACAttt                                                               >  1:680839/1‑139 (MQ=255)                                        tGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCacac                                                     >  2:616146/1‑139 (MQ=255)                                                        acacCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAg                                     <  2:680839/139‑1 (MQ=255)                                                              tGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                     >  1:1094225/1‑101 (MQ=255)                                                              tGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                     <  2:1094225/101‑1 (MQ=255)                                                                        gCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                       >  1:14679/1‑89 (MQ=255)                                                                        gCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                       >  2:848651/1‑89 (MQ=255)                                                                        gCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                       <  2:14679/89‑1 (MQ=255)                                                                        gCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                       <  1:848651/89‑1 (MQ=255)                                                                        gCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATc                     >  2:635340/1‑139 (MQ=255)                                                                                  tACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCg                                   >  1:732614/1‑115 (MQ=255)                                                                                  tACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCg                                   <  2:732614/115‑1 (MQ=255)                                                                                        aCGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGa                       >  2:582405/1‑121 (MQ=255)                                                                                        aCGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGa                       <  1:582405/121‑1 (MQ=255)                                                                                         cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGAt    >  2:1108300/1‑139 (MQ=255)                                                                                         cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGAt    >  2:149488/1‑139 (MQ=255)                                                                                          gCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGa                       >  1:390444/1‑119 (MQ=255)                                                                                          gCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGa                       <  2:390444/119‑1 (MQ=255)                                                                                              gCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTa                                                                    <  2:116852/70‑1 (MQ=255)                                                                                              gCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTa                                                                    >  1:116852/1‑70 (MQ=255)                                                                                                   cttcCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGc                   <  1:726876/114‑1 (MQ=255)                                                                                                   cttcCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGc                   >  2:726876/1‑114 (MQ=255)                                                                                                    ttcCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                       <  2:347343/61‑1 (MQ=255)                                                                                                    ttcCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                       >  1:347343/1‑61 (MQ=255)                                                                                                                        cGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATtg  >  2:230156/1‑110 (MQ=255)                                                                                                                        cGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATtg  <  1:230156/110‑1 (MQ=255)                                                                                                                          gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                     <  1:684035/41‑1 (MQ=255)                                                                                                                          gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                     >  2:684035/1‑41 (MQ=255)                                                                                                                          gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTa                                                                    <  1:76067/42‑1 (MQ=255)                                                                                                                          gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTa                                                                    >  2:76067/1‑42 (MQ=255)                                                                                                                          gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATtg  <  2:526198/108‑1 (MQ=255)                                                                                                                          gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATtg  >  1:526198/1‑108 (MQ=255)                                                                                                                             cAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACAcc                                                    >  1:881351/1‑55 (MQ=255)                                                                                                                             cAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACAcc                                                    <  2:881351/55‑1 (MQ=255)                                                                                                                                |                                                                                                      AGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAATACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTG  >  NC_000913/1663969‑1664197

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