Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,480,876	T→G	F383V (TTC→GTC)	yjgN →	conserved inner membrane protein YjgN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,480,876	0	T	G	100.0%	111.6 / NA	32	F383V (TTC→GTC)	yjgN	conserved inner membrane protein YjgN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/21);  total (11/21)

GATTACGGGCGGTCTGGCTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTG  >  NC_000913/4480746‑4481005                                                                                                                                   |                                                                                                                                  gATTACGGGCGGTCTGGCTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGcgtc                                                                                                                                 >  2:390948/1‑133 (MQ=255)     aCGGGCGGTCTGGCTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAg                                                                                                                             >  1:500650/1‑133 (MQ=255)      cGGGCGGTCTGGCTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGc                                                                                                                            >  2:95467/1‑133 (MQ=255)        ggCGGTCTGGCTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCa                                                                                                                          <  1:593778/133‑1 (MQ=255)                  cTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACg                                                                                                                <  1:581621/133‑1 (MQ=255)                    tATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACggg                                                                                                              <  1:95467/133‑1 (MQ=255)                    tATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACggg                                                                                                              <  1:16411/133‑1 (MQ=255)                      tCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGa                                                                                                            >  2:127230/1‑133 (MQ=255)                           tgctgAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTAt                                                                                                       >  2:385529/1‑133 (MQ=255)                            gctgAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATg                                                                                                      <  1:390948/133‑1 (MQ=255)                                 aaaTATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCcttc                                                                                                 <  1:127230/133‑1 (MQ=255)                                 aaaTATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCcttc                                                                                                 <  2:80670/133‑1 (MQ=255)                                     atGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTACCCACTATTTCACGGGGGGTTATGCCTTCTTTa                                                                                             >  2:201179/1‑133 (MQ=255)                                      tGGATGATTGCCTGGCAGGCAAAATATACGTACTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTAc                                                                                            <  2:158372/133‑1 (MQ=255)                                       ggATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTAcc                                                                                           <  1:385529/133‑1 (MQ=255)                                       ggATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTAcc                                                                                           <  1:116906/133‑1 (MQ=255)                                       ggATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTAcc                                                                                           >  2:474673/1‑133 (MQ=255)                                         gtgatTGACTTGCAGACAAAAAATACGTCTTTACTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCAt                                                                                         <  2:230853/132‑1 (MQ=255)                                         atgatTGAATGGCAGGCACAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCAt                                                                                         <  1:201179/133‑1 (MQ=255)                                                aCTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTa                                                                                  <  2:597375/133‑1 (MQ=255)                                                aCTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTa                                                                                  <  1:536974/133‑1 (MQ=255)                                                aCTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTa                                                                                  <  2:53076/133‑1 (MQ=255)                                                     caggcaAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAgaga                                                                             >  1:516648/1‑133 (MQ=255)                                                      aggcaAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAt                                                                            >  1:20939/1‑133 (MQ=255)                                                               aTACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGc                                                                   <  2:20939/133‑1 (MQ=255)                                                                tACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGcc                                                                  <  1:474673/133‑1 (MQ=255)                                                                  cGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTAATGCCTTCCTTATCATTTCTGTAAGAGATACAAAAGGCCgg                                                                >  2:253305/1‑133 (MQ=255)                                                                              gCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCgg                                                    >  1:578764/1‑133 (MQ=255)                                                                                            cttccttTAATCAATAAAGAAGAACAACCAGATAAAGGCGGCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                      >  2:337170/1‑133 (MQ=255)                                                                                                   tAATCAATGAAGAAGCACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTgtgt                               <  1:253305/133‑1 (MQ=255)                                                                                                   tAATCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTgtgt                               <  1:566116/133‑1 (MQ=255)                                                                                                     aTCAATAAAGAAGAACAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGttt                             >  2:348391/1‑133 (MQ=255)                                                                                                               aagaaCAACCAGATAAAGGCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                   <  1:412664/133‑1 (MQ=255)                                                                                                                                ggCGTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctg  <  2:516648/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  GATTACGGGCGGTCTGGCTTATCCACTGCTGAAAATATGGATGATTGACTGGCAGGCAAAAAATACGTATTTGCTGGGCGATTTGGATGACCTTCCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTG  >  NC_000913/4480746‑4481005

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