Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,347	C→T	*109* (TAG→TAA)	ynfA ←	UPF0060 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,347	0	C	T	100.0%	137.6 / NA	43	*109* (TAG→TAA)	ynfA	UPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/23);  total (20/23)

AAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGC  >  NC_000913/1655220‑1655480                                                                                                                                |                                                                                                                                      aaaaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                            >  1:539313/1‑139 (MQ=255)               aTCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTg                                                                                                                                 >  1:990112/1‑120 (MQ=255)               aTCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTg                                                                                                                                 <  2:990112/120‑1 (MQ=255)                              gTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCa                                                                                              >  1:117524/1‑139 (MQ=255)                                 tatTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                       <  2:1076453/112‑5 (MQ=255)                                 tatTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                       >  1:1076453/1‑108 (MQ=255)                                     tGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                       <  2:539313/139‑1 (MQ=255)                                           gtgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                          <  1:1333467/99‑2 (MQ=255)                                           gtgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                          >  2:1333467/1‑98 (MQ=255)                                               aTGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAACTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                             <  1:651592/139‑1 (MQ=255)                                                   ttGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAgcgc                                                                                        <  2:1060886/124‑1 (MQ=255)                                                   ttGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAgcgc                                                                                        >  1:1060886/1‑124 (MQ=255)                                                    tGACAAAAGATGAAATAGAATACCTCTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATCTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                        <  2:618299/139‑1 (MQ=255)                                                      aCAAAAGATGAAATAGAATGCCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                       >  1:1022418/1‑122 (MQ=255)                                                      aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                                                       <  2:1022418/122‑1 (MQ=255)                                                      aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                      >  1:54089/1‑139 (MQ=255)                                                      aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                      <  2:54089/139‑1 (MQ=255)                                                             aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATaa                                                               <  2:117524/139‑1 (MQ=255)                                                                  aTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGt                                                          <  2:1241204/139‑1 (MQ=255)                                                                             tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTg                                                                                                                                 <  1:1172501/58‑1 (MQ=255)                                                                             tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTg                                                                                                                                 >  2:1172501/1‑58 (MQ=255)                                                                              ttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACg                                                >  2:1422236/1‑137 (MQ=255)                                                                              ttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACg                                                <  1:1422236/137‑1 (MQ=255)                                                                                gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGAt                                                                                                          <  1:900648/77‑1 (MQ=255)                                                                                gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGAt                                                                                                          >  2:900648/1‑77 (MQ=255)                                                                                     cTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAGAGAGTCAGTTTCACGCCATCcaca                                       >  1:857310/1‑139 (MQ=255)                                                                                         cacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                          <  1:702107/53‑2 (MQ=255)                                                                                         cacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                          >  2:702107/1‑52 (MQ=255)                                                                                                    atcttatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGc                                                                                               <  2:1206031/68‑1 (MQ=255)                                                                                                    atcttatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGc                                                                                               >  1:1206031/1‑68 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                                 >  1:1371747/1‑45 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                                 <  2:1371747/45‑1 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                        >  1:1000007/1‑86 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                        >  1:274420/1‑86 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                        <  2:274420/86‑1 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATcc                                          >  2:269776/1‑116 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATcc                                          <  1:269776/116‑1 (MQ=255)                                                                                                         atTGATAAAATGGATTTATGTTATTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                                        <  2:1000007/86‑1 (MQ=255)                                                                                                            gATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGtt                                                     <  1:1312004/102‑1 (MQ=255)                                                                                                            gATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGtt                                                     >  2:1312004/1‑102 (MQ=255)                                                                                                               aaaaTGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATg                                                                                                   <  2:528710/53‑1 (MQ=255)                                                                                                               aaaaTGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATg                                                                                                   >  1:528710/1‑53 (MQ=255)                                                                                                                          aTGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGc  <  1:1384227/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                      AAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGC  >  NC_000913/1655220‑1655480

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