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%	132.6 / NA	37	I163T (ATA→ACA)	ynfH	oxidoreductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/17);  total (20/17)

TTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAATACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCAC  >  NC_000913/1663965‑1664229                                                                                                                                    |                                                                                                                                      tttCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTtgct                                                                                                                                >  1:367052/1‑139 (MQ=255)          ggCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCa                                                                                                                          >  1:476301/1‑136 (MQ=255)          ggCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCa                                                                                                                          <  2:476301/136‑1 (MQ=255)            cgcTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCt                                                                                                                     >  1:286460/1‑139 (MQ=255)                 ggCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCggg                                                                                                                >  2:389926/1‑139 (MQ=255)                            cGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                     <  1:388693/139‑1 (MQ=255)                            cGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                     <  1:389926/139‑1 (MQ=255)                                         aCCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACca                                                                                        <  2:367052/139‑1 (MQ=255)                                               gtgtATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCAt                                                                                  <  2:738517/139‑1 (MQ=255)                                                   aTCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                                                       >  1:412796/1‑114 (MQ=255)                                                   aTCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTcgcg                                                                                                       <  2:412796/114‑1 (MQ=255)                                                          cGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCAtt                                                                       >  1:429319/1‑139 (MQ=255)                                                            acacCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACa                                                                                                  <  1:77930/110‑1 (MQ=255)                                                            acacCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACa                                                                                                  >  2:77930/1‑110 (MQ=255)                                                                  tGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCa                                                                                                                          >  1:105133/1‑80 (MQ=255)                                                                  tGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCa                                                                                                                          <  2:105133/80‑1 (MQ=255)                                                                           ggCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGAt                                                      >  2:600473/1‑139 (MQ=255)                                                                                aaCGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCtt                                                 <  2:429319/139‑1 (MQ=255)                                                                                          cTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGt                                       >  1:123509/1‑139 (MQ=255)                                                                                            aCGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGa                                     <  2:123509/139‑1 (MQ=255)                                                                                             cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGAt                                    >  1:694782/1‑139 (MQ=255)                                                                                             cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGAt                                    >  1:454641/1‑139 (MQ=255)                                                                                                  gCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGa                                     <  2:886742/133‑1 (MQ=255)                                                                                                  gCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGa                                     >  1:886742/1‑133 (MQ=255)                                                                                                       cttcCTGACGGTATTGTTGAGCGGTCCAACACTtgc                                                                                                                                 <  2:389094/36‑1 (MQ=255)                                                                                                       cttcCTGACGGTATTGTTGAGCGGTCCAACACTtgc                                                                                                                                 >  1:389094/1‑36 (MQ=255)                                                                                                       cttcCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCaa                          >  1:190173/1‑139 (MQ=255)                                                                                                                  tATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGt                                       <  2:672619/115‑1 (MQ=255)                                                                                                                  tATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGt                                       >  1:672619/1‑115 (MQ=255)                                                                                                                     tgttgAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGc            <  2:694782/139‑1 (MQ=255)                                                                                                                          aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTc                                                                                                          >  1:446577/1‑40 (MQ=255)                                                                                                                          aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTc                                                                                                          <  2:446577/40‑1 (MQ=255)                                                                                                                          aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGt                                       >  1:649331/1‑107 (MQ=255)                                                                                                                          aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGt                                       <  2:649331/107‑1 (MQ=255)                                                                                                                              gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCa   <  2:484666/139‑1 (MQ=255)                                                                                                                               tCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCAc  >  1:193913/1‑139 (MQ=255)                                                                                                                                 cAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGAGTGATTGTGCTGCAAGGTTTGTCTCTGGCGt          >  1:139925/1‑129 (MQ=255)                                                                                                                                    |                                                                                                                                      TTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAATACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCAC  >  NC_000913/1663965‑1664229

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