Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,236,508	A→G	T83A (ACT→GCT)	yjbF →	lipoprotein YjbF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,236,508	0	A	G	100.0%	83.1 / NA	25	T83A (ACT→GCT)	yjbF	lipoprotein YjbF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/16);  total (9/16)

CTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACG  >  NC_000913/4236376‑4236642                                                                                                                                     |                                                                                                                                       cTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGt                                                                                                                                 <  1:77359/140‑1 (MQ=255)   gACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTa                                                                                                                               <  1:428472/140‑1 (MQ=255)   gACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTa                                                                                                                               <  2:365544/140‑1 (MQ=255)              atatTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGGGCTGGCCTTCCCTGAAAACGGACAAAAAAAATGGGTCACTAAGGATCAGGCGCCCCCCTTTACCCACAATGGc                                                                                                                    >  2:603250/1‑140 (MQ=255)               tatTCAAAATATGCCCTTCGCCCGCCAGTACCTGCAGCTTAATGGCGGCCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGcc                                                                                                                   <  2:193701/140‑1 (MQ=255)                 ttCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGt                                                                                                                 <  2:736234/140‑1 (MQ=255)                  tCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTc                                                                                                                >  1:367943/1‑140 (MQ=255)                  tCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTc                                                                                                                >  2:537942/1‑140 (MQ=255)                   cAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCt                                                                                                               <  2:591573/140‑1 (MQ=255)                      aaTATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATCGGTCACTCCGGATCAGTCTGCTCTCTTTACACAACATGGCCGTCTggg                                                                                                            >  2:231991/1‑139 (MQ=255)                          tGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAg                                                                                                        >  1:696732/1‑140 (MQ=255)                           gCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCCCTCAGGATGAGGGTGCTCTCGTTACACAACATGGCCGTCTGGTGAAGa                                                                                                       >  2:235131/1‑140 (MQ=255)                           gCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGa                                                                                                       >  1:514422/1‑140 (MQ=255)                                          gTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGgcggcg                                                                                        <  2:367943/140‑1 (MQ=255)                                           tACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGa                                                                                       >  1:45229/1‑140 (MQ=255)                                              aTGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACaa                                                                                    <  2:503322/140‑1 (MQ=255)                                                gCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAAcc                                                                                  >  2:37251/1‑140 (MQ=255)                                                    cTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGA‑‑CTGCTCGGCGGCGACAACCTGATTg                                                                            >  2:718553/1‑140 (MQ=255)                                                    cTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGA‑‑CTGCTCGGCGGCGACAACCTGATTg                                                                            <  1:718553/140‑1 (MQ=255)                                                     ttAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACGCAACATGGGCGTCTGGAGAAGACTCCGCTCGGCGGCGACAACCTGAtt                                                                             >  2:235063/1‑140 (MQ=255)                                                       aaTGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGa                                                                           <  1:631210/140‑1 (MQ=255)                                                       aaTGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGa                                                                           <  2:696732/140‑1 (MQ=255)                                                       aaTGGCGGGCCGCAGTTATTTGGGGTGCTGGTCTTCGCTGAAGACGGATAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGa                                                                           <  2:193098/140‑1 (MQ=255)                                                        aTGGCGGGCCGCAGTTACTTGTCGTGCTGGCATTCTCTGTAGACGGACTACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGaa                                                                          <  2:206697/140‑1 (MQ=255)                                                           gCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTg                                                                       <  1:537942/140‑1 (MQ=255)                                                                                                          caaAAATGGGTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGa                        <  2:644122/140‑1 (MQ=255)                                                                                                                  ggTCACTCAGGATCAGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCaa                >  2:451864/1‑140 (MQ=255)                                                                                                                                aGGCTGCTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACg  <  1:603250/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACG  >  NC_000913/4236376‑4236642

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