Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,794,172	C→T	*327* (TAG→TAA)	btuC ←	vitamin B12 ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,794,172	0	C	T	100.0%	134.5 / NA	42	*327* (TAG→TAA)	btuC	vitamin B12 ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/20);  total (22/20)

AACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGC  >  NC_000913/1794042‑1794299                                                                                                                                   |                                                                                                                                aaCTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTcc                                                                                                                         >  2:380996/1‑139 (MQ=255)   cTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTg                                                                                                                       >  1:284118/1‑139 (MQ=255)           gCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGA‑CGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                              <  1:1170796/139‑1 (MQ=255)            cGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                              >  2:56813/1‑139 (MQ=255)               ggTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                           <  1:1100537/139‑1 (MQ=255)                               atatCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTaa                                                                                                                <  1:516637/118‑1 (MQ=255)                               atatCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTaa                                                                                                                >  2:516637/1‑118 (MQ=255)                               atatCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCgg                                                                                           >  1:1255726/1‑139 (MQ=255)                                tatCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          <  2:1186640/139‑1 (MQ=255)                                tatCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          <  2:796492/139‑1 (MQ=255)                                       cACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                           <  2:596935/115‑1 (MQ=255)                                       cACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                           >  1:596935/1‑115 (MQ=255)                                          tACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAcc                                                                                                        <  1:172522/115‑1 (MQ=255)                                          tACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAcc                                                                                                        >  2:172522/1‑115 (MQ=255)                                            cGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                   >  1:92481/1‑102 (MQ=255)                                            cGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                   <  2:92481/102‑1 (MQ=255)                                                    aGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          <  2:681636/87‑1 (MQ=255)                                                    aGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          >  1:681636/1‑87 (MQ=255)                                                    aGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAc                                                                      <  1:801535/139‑1 (MQ=255)                                                         ggAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGc                                                                 >  2:866450/1‑139 (MQ=255)                                                                tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          >  1:400394/1‑75 (MQ=255)                                                                tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          <  2:400394/75‑1 (MQ=255)                                                                tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAgg                                                          >  2:1025481/1‑139 (MQ=255)                                                                           tcCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCgg                                               <  2:156531/139‑1 (MQ=255)                                                                              tGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAg                                            <  1:380996/139‑1 (MQ=255)                                                                                         aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                 <  1:1025481/139‑1 (MQ=255)                                                                                              aaaTCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcgcgc                            <  2:1255726/139‑1 (MQ=255)                                                                                                     aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          <  2:468191/38‑1 (MQ=38)                                                                                                     aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          >  1:468191/1‑38 (MQ=255)                                                                                                                  aGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATAACGGCcagcagga        <  2:284118/139‑3 (MQ=255)                                                                                                                    aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGTGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt      >  1:10214/1‑139 (MQ=255)                                                                                                                    aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt      >  1:917369/1‑139 (MQ=255)                                                                                                                       gTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          <  2:352187/52‑1 (MQ=255)                                                                                                                       gTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          >  1:352187/1‑52 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGTGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  <  2:10214/139‑1 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                 >  1:1093746/1‑108 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                 <  2:1093746/108‑1 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  >  2:369069/1‑139 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  >  2:1291937/1‑139 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  >  1:121939/1‑139 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  >  2:790943/1‑139 (MQ=255)                                                                                                                        tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  >  1:1048089/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                AACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGC  >  NC_000913/1794042‑1794299

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