Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,638,026	(C)6→5	coding (627/1179 nt)	bamB ←	outer membrane protein assembly factor BamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,638,021	0	C	.	100.0%	79.5 / NA	20	coding (632/1179 nt)	bamB	outer membrane protein assembly factor BamB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (6/14);  total (6/14)

GACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTG  >  NC_000913/2637897‑2638145                                                                                                                             |                                                                                                                             gacgGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGcccccacga                                                                                                                      <  1:433616/133‑5 (MQ=255)   cgGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACCGCGGTGATGCGACCATAATCG‑TCCCCacgacg                                                                                                                   >  2:202383/1‑133 (MQ=255)             gtCAACATCGCTCAGTCGGTCAATTTCGATAGTACCGGTCGTCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCaa                                                                                                         <  2:240744/133‑1 (MQ=255)                      gCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTg                                                                                                <  1:30474/133‑1 (MQ=255)                       cTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGt                                                                                               <  1:41248/133‑1 (MQ=255)                       cTCAGACGGTCAATTCCGGTAGCACCGGTCACCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGt                                                                                               <  1:202383/133‑1 (MQ=255)                             cGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGgcgc                                                                                         >  1:123678/1‑133 (MQ=255)                                tCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATAAGCACGGCGCTGCCGCGCCAATTATCG‑CCCCAACGCCGCCCGCACCAAAAGCGGTTGCCGGCGCAGa                                                                                      >  2:339182/1‑133 (MQ=255)                                         gTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACg                                                                             >  2:601534/1‑133 (MQ=255)                                                  gTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAagagag                                                                    <  2:455619/133‑1 (MQ=255)                                                    cGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCg                                                                  >  2:588311/1‑133 (MQ=255)                                                            aaaTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCata                                                          <  1:90583/133‑1 (MQ=255)                                                            aaaTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCata                                                          <  1:588311/133‑1 (MQ=255)                                                             aaTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCatat                                                         <  1:601534/133‑1 (MQ=255)                                                                 cGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAg                                                     <  1:537236/133‑1 (MQ=255)                                                                      ctcccAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTaa                                                >  1:563811/4‑133 (MQ=255)                                                                           aaaTCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTc                                           <  2:614159/133‑1 (MQ=255)                                                                                          gTTCCATCAGCACTGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGCGCGAAGGCGTATGGAGTTTAACTGTGCATTGGACAGCGCCg                            >  2:172579/1‑133 (MQ=255)                                                                                                       tGCGCTGACGCGACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCa               <  1:339182/133‑1 (MQ=255)                                                                                                           cTGACGCGAACATTAACG‑CCCCCAGGGCGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAgcgc           <  2:286107/133‑1 (MQ=255)                                                                                                                  gACCATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAAc    <  2:123678/133‑1 (MQ=255)                                                                                                                    ccATTATCG‑CCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTg  <  1:619557/133‑1 (MQ=255)                                                                                                                             |                                                                                                                             GACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTG  >  NC_000913/2637897‑2638145

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