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%	123.4 / NA	39	*109* (TAG→TAA)	ynfA	UPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

CGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCA  >  NC_000913/1655217‑1655469                                                                                                                                   |                                                                                                                           cGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTgcg                                                                                                                    >  1:384742/1‑139 (MQ=255)            ggAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGc                                                                                                        <  2:207897/139‑1 (MQ=255)              aGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                      <  1:560705/139‑1 (MQ=255)                   tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                 >  2:484755/1‑124 (MQ=255)                   tCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTATCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                 <  1:484755/124‑1 (MQ=255)                          aaaaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                      <  2:516100/127‑1 (MQ=255)                          aaaaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                      >  1:516100/1‑127 (MQ=255)                              ggCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGcc                                                                                      >  2:335508/1‑139 (MQ=255)                                         gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTGTCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCaa                                                                           >  2:477227/1‑139 (MQ=255)                                               tgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACg                                                                     >  2:9276/1‑139 (MQ=255)                                               tgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACg                                                                     <  2:738464/139‑1 (MQ=255)                                                    ggTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAccc                                                                >  1:107538/1‑139 (MQ=255)                                                       tGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                                             <  1:9276/139‑1 (MQ=255)                                                         aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                      <  2:655629/96‑1 (MQ=255)                                                         aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                      >  1:655629/1‑96 (MQ=255)                                                         aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                           >  1:130122/1‑139 (MQ=255)                                                         aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                           >  1:363443/1‑139 (MQ=255)                                                                aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                  >  1:487834/1‑125 (MQ=255)                                                                aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAc                                                                  <  2:487834/125‑1 (MQ=255)                                                                   aaaTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                            >  2:445688/1‑77 (MQ=255)                                                                   aaaTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                            <  1:445688/81‑5 (MQ=255)                                                                              cTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGTCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCAc                                      <  2:384742/139‑1 (MQ=255)                                                                                tttGTCAGCTGCCACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTTCGCCCCC‑GCCCGCAACAATGATCAACATGCCGCAATGCGCAATCAACGCACCCGTCCAGTCATAAAGAGGCAGTTTCACGc                                    >  1:308801/1‑139 (MQ=255)                                                                                tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCGAAGCGCAATCAACGCACCCg                                                               >  1:363189/1‑112 (MQ=255)                                                                                tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCGAAGCGCAATCAACGCACCCg                                                               <  2:363189/112‑1 (MQ=255)                                                                                                           tatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg         >  2:260222/1‑139 (MQ=255)                                                                                                           tatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACg         <  1:477227/139‑1 (MQ=255)                                                                                                             ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaaca                                                                                                    <  1:652638/46‑1 (MQ=255)                                                                                                             ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaaca                                                                                                    >  2:652638/1‑46 (MQ=255)                                                                                                             ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                   >  1:254398/1‑127 (MQ=255)                                                                                                             ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                   <  2:254398/127‑1 (MQ=255)                                                                                                             ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGcc       >  2:599732/1‑139 (MQ=255)                                                                                                                aTAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTg    <  2:130122/139‑1 (MQ=255)                                                                                                                aTAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTg    <  1:335508/139‑1 (MQ=255)                                                                                                                  aaaaTGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCa  >  2:420909/1‑139 (MQ=255)                                                                                                                    aaTGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACCATGCGCAGCCACATCACCACCGTg    <  2:490916/135‑1 (MQ=255)                                                                                                                    aaTGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAATGCGCAGCCACATCAACACCGTg    >  1:490916/1‑135 (MQ=255)                                                                                                                          tttATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGAt                                                                                               >  2:642540/1‑38 (MQ=255)                                                                                                                          tttATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGAt                                                                                               <  1:642540/38‑1 (MQ=255)                                                                                                                                   |                                                                                                                           CGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCA  >  NC_000913/1655217‑1655469

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: ‑