Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,798,659	(T)8→7	coding (654/1074 nt)	waaU ←	lipopolysaccharide core biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,798,652	0	T	.	100.0%	152.0 / NA	34	coding (661/1074 nt)	waaU	lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/18);  total (16/18)

ATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATT  >  NC_000913/3798509‑3798791                                                                                                                                                |                                                                                                                                            atcaAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttt                                                                                                                                        >  1:16590/1‑149 (MQ=255)     aaTGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttgca                                                                                                                                    >  4:108796/1‑146 (MQ=255)                   gATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATtaat                                                                                                                     <  2:16590/149‑1 (MQ=255)                       tCAAGTATTGGTATTGTCAATAATTCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGTCATA‑TTTTTTTGCACCTAATGGATTaataat                                                                                                                  <  8:123914/148‑1 (MQ=255)                                 gtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAAt                                                                                                                             >  6:11964/1‑127 (MQ=255)                                 gtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAAt                                                                                                                             <  5:11964/127‑1 (MQ=255)                                  tattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGt                                                                                                      <  1:42128/149‑1 (MQ=255)                                      gTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAtct                                                                                                  >  3:8061/1‑149 (MQ=255)                                       tCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAtctc                                                                                                 >  7:21071/1‑149 (MQ=255)                                         aaTAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTcc                                                                                                >  7:74609/1‑148 (MQ=255)                                                    cGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCtt                                                                                     <  3:26587/148‑1 (MQ=255)                                                     ggTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATg                                                                                                                            >  3:90496/1‑108 (MQ=255)                                                     ggTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATg                                                                                                                            <  4:90496/108‑1 (MQ=255)                                                         aCCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttattt                                                                               <  4:8061/149‑1 (MQ=255)                                                            cGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCtt                                                                                                        >  7:43725/1‑121 (MQ=255)                                                            cGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCtt                                                                                                        <  8:43725/121‑1 (MQ=255)                                                                 aaTATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTc                                                                        <  8:74609/148‑1 (MQ=255)                                                                        tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTa                                                                >  4:119592/1‑149 (MQ=255)                                                                                      aaGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAg                                                  >  8:17098/1‑149 (MQ=255)                                                                                      aaGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAg                                                  <  7:17098/149‑1 (MQ=255)                                                                                             ttttCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCTTA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCCATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGt                                           >  8:115524/1‑149 (MQ=255)                                                                                                 cACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTc                                                                                                          >  5:35296/1‑82 (MQ=255)                                                                                                 cACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTc                                                                                                          <  6:35296/82‑1 (MQ=255)                                                                                                  aCTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCtttta                                                                                  >  6:67378/1‑105 (MQ=255)                                                                                                  aCTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCtttta                                                                                  <  5:67378/105‑1 (MQ=255)                                                                                                                   ttttATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttat                                                                                 <  4:29189/89‑1 (MQ=255)                                                                                                                   ttttATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttat                                                                                 >  3:29189/1‑89 (MQ=255)                                                                                                                   ttttATTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCAc                      >  8:99370/1‑148 (MQ=255)                                                                                                                       aTTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTtttattttat                                                                            <  2:74759/90‑1 (MQ=255)                                                                                                                       aTTTGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTtttattttat                                                                            >  1:74759/1‑90 (MQ=255)                                                                                                                        tttGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAg                                            <  7:33391/121‑1 (MQ=255)                                                                                                                        tttGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAg                                            >  8:33391/1‑121 (MQ=255)                                                                                                                        tttGCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATaa                 <  3:119592/148‑1 (MQ=255)                                                                                                                           gCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTcc                                                                                                <  4:91681/66‑1 (MQ=255)                                                                                                                           gCTCAAACGTAAGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTcc                                                                                                >  3:91681/1‑66 (MQ=255)                                                                                                                                      aGACGGCATA‑TTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTAtt  <  8:21071/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                            ATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATT  >  NC_000913/3798509‑3798791

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