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%	167.3 / NA	52	*327* (TAG→TAA)	btuC	vitamin B12 ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/28);  total (24/28)

GAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTC  >  NC_000913/1794041‑1794301                                                                                                                                    |                                                                                                                                  gAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            <  2:298030/139‑1 (MQ=255)            gCGTGGTCACTTCACCGTTGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAAc                                                                                                                 <  2:6677/139‑1 (MQ=255)                     cTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAg                                                                                                        <  1:922034/139‑1 (MQ=255)                            gTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                                <  2:684647/124‑1 (MQ=255)                            gTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                                >  1:684647/1‑124 (MQ=255)                                                tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTa                                                                                                                                 <  1:426450/87‑1 (MQ=255)                                                  gtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            >  1:279368/1‑90 (MQ=255)                                                  gtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            <  2:279368/90‑1 (MQ=255)                                                  gtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                            <  2:532913/122‑1 (MQ=255)                                                  gtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                            >  1:532913/1‑122 (MQ=255)                                                   tcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCt                                                                                                                          >  2:430626/1‑91 (MQ=255)                                                   tcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCt                                                                                                                          <  1:430626/91‑1 (MQ=255)                                                          ggAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCgg                                                                                             >  2:39855/1‑113 (MQ=255)                                                          ggAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCgg                                                                                             <  1:39855/113‑1 (MQ=255)                                                                 tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACaa                                                                                                               >  2:313247/1‑88 (MQ=255)                                                                 tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACaa                                                                                                               <  1:313247/88‑1 (MQ=255)                                                                   cATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                      >  1:474819/1‑127 (MQ=255)                                                                   cATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                      <  2:474819/127‑1 (MQ=255)                                                                        tttCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGctct                                                     >  1:898553/1‑139 (MQ=255)                                                                                 ggTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCt                                            >  2:109213/1‑139 (MQ=255)                                                                                   tgtgGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTaa                                          >  1:919181/1‑139 (MQ=255)                                                                                      ggTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGc                                       >  1:722992/1‑139 (MQ=255)                                                                                          aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                   <  1:77620/139‑1 (MQ=255)                                                                                          aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                   <  2:118099/139‑1 (MQ=255)                                                                                           gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                     <  1:339553/56‑1 (MQ=255)                                                                                           gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                     >  2:339553/1‑56 (MQ=255)                                                                                           gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                           >  2:201838/1‑66 (MQ=255)                                                                                           gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                           <  1:201838/66‑1 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCgg                                                 <  1:150770/123‑1 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCgg                                                 >  2:150770/1‑123 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 <  2:898553/139‑1 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 <  1:149261/139‑1 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgc                                  >  1:859372/1‑138 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgc                                  <  2:859372/138‑1 (MQ=255)                                                                                                 aTCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATaa                                                                                                            >  1:792119/1‑59 (MQ=255)                                                                                                 aTCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATaa                                                                                                            <  2:792119/59‑1 (MQ=255)                                                                                                           ttAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGcc                                      >  1:639545/1‑119 (MQ=255)                                                                                                           ttAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGcc                                      <  2:639545/119‑1 (MQ=255)                                                                                                           ttAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 <  1:682881/124‑1 (MQ=255)                                                                                                           ttAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 >  2:682881/1‑124 (MQ=255)                                                                                                            tAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                            <  1:471847/64‑1 (MQ=255)                                                                                                            tAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                            >  2:471847/1‑64 (MQ=255)                                                                                                                     aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTa                            >  2:60773/1‑119 (MQ=255)                                                                                                                     aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTa                            <  1:60773/119‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCAcc                                                                                       >  1:292722/1‑56 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCAcc                                                                                       <  2:292722/56‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 <  1:184115/110‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 >  2:184115/1‑110 (MQ=255)                                                                                                                         tCTGTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc    >  1:631579/1‑139 (MQ=255)                                                                                                                           tttttGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc    >  2:629184/1‑137 (MQ=255)                                                                                                                           tttttGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc    <  1:629184/137‑1 (MQ=255)                                                                                                                           tttttGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTc  <  2:722992/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                  GAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTC  >  NC_000913/1794041‑1794301

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