Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,445,549	A→G	T480A (ACC→GCC)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,445,549	0	A	G	100.0%	99.2 / NA	29	T480A (ACC→GCC)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/16);  total (13/16)

GCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTC  >  NC_000913/4445432‑4445677                                                                                                                      |                                                                                                                                 gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACg                                                                                                                   >  2:93164/1‑133 (MQ=255) gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACg                                                                                                                   >  2:89868/1‑133 (MQ=255) gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACg                                                                                                                   >  2:36537/1‑133 (MQ=255)    aTTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGc                                                                                                                <  1:89868/133‑1 (MQ=255)          tGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGt                                                                                                          <  2:516067/133‑1 (MQ=255)          tGGCGCGGCGAGCTATCGATTACCGGCATTGACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGt                                                                                                          <  1:241676/133‑1 (MQ=255)           ggCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCGGGCGGAAGGACGTACCGGAGtt                                                                                                         >  1:150603/1‑133 (MQ=255)                            cTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACg                                                                                        <  1:93164/133‑1 (MQ=255)                            cTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACg                                                                                        >  1:487036/1‑133 (MQ=255)                                      ttAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGaa                                                                              <  2:331824/133‑1 (MQ=255)                                           acCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAg                                                                         >  2:296322/1‑133 (MQ=255)                                             cGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGGCGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTg                                                                       >  2:487978/1‑133 (MQ=255)                                              gCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGa                                                                      <  2:101008/133‑1 (MQ=255)                                              gCCAAAGAGTACCCGCACTGGCCGTCGACACTCAATGGCTTCATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGa                                                                      <  2:216166/133‑1 (MQ=255)                                                       ttCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACg                                                             >  1:484594/1‑133 (MQ=255)                                                         cccGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGc                                                           >  1:351789/1‑133 (MQ=255)                                                                     gTCGAAACTAAATGGCTTGATTAAAAGCGGCGTTAGCCTGTTCGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGc                                               <  2:246646/133‑1 (MQ=255)                                                                       cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  1:36537/133‑1 (MQ=255)                                                                       cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  1:310024/133‑1 (MQ=255)                                                                       cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  1:303910/133‑1 (MQ=255)                                                                            cTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGtt                                        <  1:394189/133‑1 (MQ=255)                                                                               aaTGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTtata                                     <  1:296322/133‑1 (MQ=255)                                                                                    cTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCCGAAAGGCAACAGTGATATGCAg                                >  1:166287/1‑133 (MQ=255)                                                                                         ttAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                           <  2:277917/133‑1 (MQ=255)                                                                                              aCCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATcc                      <  2:487036/133‑1 (MQ=255)                                                                                                ccGCGGTAGCCTGTACGGCGGCGCTTGGTAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCa                    <  2:172041/133‑1 (MQ=255)                                                                                                  gcgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAgg                  >  2:501480/1‑133 (MQ=255)                                                                                                   cgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAggg                 >  1:506547/1‑133 (MQ=255)                                                                                                                  cggcGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCACTGGATGATCCCAGGACTTCATGTGGAACTc  >  1:242847/1‑133 (MQ=255)                                                                                                                  cggcGCCTGGCAGATGGACGGGCCGGAGTTGAAGCGGACCGGTAACGTTAAACAGAACAAAGTGAACGGTGACGGCACGCTGAAAGGCAACAGGTATATGCAGTGGATGATCCCACCGCGTCATCTGGAACTc  >  2:336495/1‑133 (MQ=255)                                                                                                                      |                                                                                                                                 GCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTC  >  NC_000913/4445432‑4445677

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