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

CATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAATACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCACT  >  NC_000913/1663971‑1664230                                                                                                                              |                                                                                                                                       cATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCa                                                                                                                           <  2:192224/139‑1 (MQ=255)    ggCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTCTTGTTGAGCGGTCCAACACTTGCTGCCGCAAtt                                                                                                                        >  2:139071/1‑139 (MQ=255)      cgcTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCt                                                                                                                      >  2:720904/1‑139 (MQ=255)               tCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGc                                                                                                             <  2:637681/139‑1 (MQ=255)                           ggATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACAtt                                                                                                 >  2:134752/1‑139 (MQ=255)                            gatgatGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACAttt                                                                                                >  1:708128/1‑139 (MQ=255)                              tgatgaCCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTaa                                                                                              <  1:32903/139‑1 (MQ=255)                                   aCCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACca                                                                                         <  1:739814/139‑1 (MQ=255)                                      tGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCacac                                                                                      >  1:156496/1‑139 (MQ=255)                                       gCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACca                                                                                         >  2:595542/1‑135 (MQ=255)                                       gCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACca                                                                                         <  1:595542/135‑1 (MQ=255)                                        cGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCa                                                                                    >  1:328488/1‑139 (MQ=255)                                             aTCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGc                                                                               <  1:201106/139‑1 (MQ=255)                                             aTCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGc                                                                               >  2:201106/1‑139 (MQ=255)                                                          cGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACca                                                                                         >  1:58064/1‑116 (MQ=255)                                                          cGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACca                                                                                         <  2:58064/116‑1 (MQ=255)                                                                           aCGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCggg                                                                                                                 >  1:181844/1‑75 (MQ=255)                                                                           aCGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCggg                                                                                                                 <  2:181844/75‑1 (MQ=255)                                                                            cGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACAtt                                                                                                 <  1:375296/90‑1 (MQ=255)                                                                            cGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACAtt                                                                                                 >  2:375296/1‑90 (MQ=255)                                                                                 acacTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTaa                                                                                              <  2:284009/88‑1 (MQ=255)                                                                                 acacTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTaa                                                                                              >  1:284009/1‑88 (MQ=255)                                                                                 acacTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCgg                                           >  1:436170/1‑139 (MQ=255)                                                                                       cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                      >  2:312473/1‑74 (MQ=255)                                                                                       cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                      <  1:312473/74‑1 (MQ=255)                                                                                       cGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGAt                                     >  1:259983/1‑139 (MQ=255)                                                                                         cTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCtt                                                  <  1:654699/124‑1 (MQ=255)                                                                                         cTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCtt                                                  >  2:654699/1‑124 (MQ=255)                                                                                                               tgttgAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGc             <  2:436170/139‑1 (MQ=255)                                                                                                                gttgAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCg            <  2:708128/139‑1 (MQ=255)                                                                                                                   gAGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCt         <  1:134752/139‑1 (MQ=255)                                                                                                                    aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                      >  1:319711/1‑45 (MQ=255)                                                                                                                    aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                      <  2:319711/45‑1 (MQ=255)                                                                                                                    aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAGTACCACGCCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTa        <  1:139071/139‑1 (MQ=255)                                                                                                                    aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCt         <  2:625895/138‑1 (MQ=255)                                                                                                                    aGCGGTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCt         >  1:625895/1‑138 (MQ=255)                                                                                                                        gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTg                                                         >  2:635404/1‑86 (MQ=255)                                                                                                                        gTCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTg                                                         <  1:635404/86‑1 (MQ=255)                                                                                                                         tCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                      >  2:225585/1‑40 (MQ=255)                                                                                                                         tCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGtt                                                                                                      <  1:225585/40‑1 (MQ=255)                                                                                                                         tCCAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCAc   >  2:482555/1‑139 (MQ=255)                                                                                                                          ccAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCACt  <  1:768492/139‑1 (MQ=255)                                                                                                                          ccAACACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCACt  <  1:482555/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                       CATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTGAGCGGTCCAATACTTGCTGCCGCAATTCTGCGGGCAGCTCGCGTTACATTTAATACCACACCATTTGCCATCATTAGCGTTCTGGCGCTGATCGCTTGTGCCGGGGTGATTGTGCTGCAAGGTTTGTCTCTGGCGTCTATTCACT  >  NC_000913/1663971‑1664230

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