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%	116.8 / NA	33	T480A (ACC→GCC)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/21);  total (12/21)

CGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTC  >  NC_000913/4445419‑4445677                                                                                                                                   |                                                                                                                                 cGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGcc                                                                                                                                <  1:77920/133‑1 (MQ=255)            aGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGa                                                                                                                     <  2:55940/133‑1 (MQ=255)                 aTTAGTTGGCGCGGTGAGCTAGCGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAATCCGCGGAAACCTGTACGGCGGCGCCCGGCAGATGGACGTGc                                                                                                                >  2:205708/1‑133 (MQ=255)                        ggCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCTCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGtt                                                                                                         >  1:358798/1‑133 (MQ=255)                        ggCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGtt                                                                                                         >  2:530715/1‑133 (MQ=255)                        ggCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGtt                                                                                                         >  2:47191/1‑133 (MQ=255)                             ggTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGc                                                                                                    >  1:137018/1‑133 (MQ=255)                                 aGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGAc                                                                                                <  1:530715/133‑1 (MQ=255)                                    tAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAGAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCgg                                                                                             <  2:12268/133‑1 (MQ=255)                                    tAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCgg                                                                                             <  2:346628/133‑1 (MQ=255)                                           tAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGtt                                                                                      <  2:137018/133‑1 (MQ=255)                                                           gCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGa                                                                      >  1:445917/1‑133 (MQ=255)                                                                             tGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGa                                                    <  1:47191/133‑1 (MQ=255)                                                                                 cGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAgg                                                <  1:116682/133‑1 (MQ=255)                                                                                  gTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGc                                               <  2:435584/133‑1 (MQ=255)                                                                                    cGAAACTCAATGGCTTGATTAAACTCTGTGGTCGCCTGTACAGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  2:292154/133‑1 (MQ=255)                                                                                    cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  1:454918/133‑1 (MQ=255)                                                                                    cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  1:234964/133‑1 (MQ=255)                                                                                    cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  2:358798/133‑1 (MQ=255)                                                                                    cGAAACTCAATGGCTTGATTAAAAACCGCGGTAGCCTGTACTGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                             <  2:206923/133‑1 (MQ=255)                                                                                        aCTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGt                                         <  2:445917/133‑1 (MQ=255)                                                                                            aaTGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTtata                                     <  2:557252/133‑1 (MQ=255)                                                                                              tGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATg                                   <  1:515111/133‑1 (MQ=255)                                                                                              tGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATg                                   >  2:95431/1‑133 (MQ=255)                                                                                                      ttAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                           <  1:70618/133‑1 (MQ=255)                                                                                                           aCCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTCATATGCAGTCGATGATcc                      >  1:258841/1‑133 (MQ=255)                                                                                                              cgcgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAg                   >  2:427439/1‑133 (MQ=255)                                                                                                               gcgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                           <  1:579806/124‑1 (MQ=255)                                                                                                               gcgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                           >  2:579806/1‑124 (MQ=255)                                                                                                                cgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                           >  1:363826/1‑123 (MQ=255)                                                                                                                cgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                           <  2:363826/123‑1 (MQ=255)                                                                                                                  gTAGCCTGGACGGCGCCGCCTGGCAGATGGTCGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCt               <  1:205708/133‑1 (MQ=255)                                                                                                                               cggcGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTc  >  2:469547/1‑133 (MQ=255)                                                                                                                               cggcGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTc  >  2:492614/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                                 CGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTC  >  NC_000913/4445419‑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: ‑