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

GAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGC  >  NC_000913/1794041‑1794299                                                                                                                                    |                                                                                                                                gAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          <  2:807588/139‑1 (MQ=255)  aaCTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTcc                                                                                                                         <  2:553164/139‑1 (MQ=255)                 gacACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAAGCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          <  1:1042605/121‑1 (MQ=255)                 gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATAGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                          >  2:1042605/1‑123 (MQ=255)                                 tatCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          <  1:752918/139‑1 (MQ=255)                                           tACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACg                                                                                >  2:1557980/1‑139 (MQ=255)                                                tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                              <  2:1305719/104‑1 (MQ=255)                                                tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                              >  1:1305719/1‑104 (MQ=255)                                                            aaTCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTa                                                                                                                               >  1:769454/1‑75 (MQ=255)                                                            aaTCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTa                                                                                                                               <  2:769454/75‑1 (MQ=255)                                                              tCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAAc                                                                                                               >  2:313929/1‑89 (MQ=255)                                                              tCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAAc                                                                                                               <  1:313929/89‑1 (MQ=255)                                                                 tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAgg                                                          >  2:136261/1‑139 (MQ=255)                                                                          tctcCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGc                                                 >  2:1542036/1‑139 (MQ=255)                                                                            tcCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCgg                                               <  2:1331830/139‑1 (MQ=255)                                                                                 ggTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCt                                                                                                                        >  2:56254/1‑61 (MQ=255)                                                                                 ggTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCt                                                                                                                        <  1:56254/61‑1 (MQ=255)                                                                                  gtgtGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTa                                         >  1:689323/1‑139 (MQ=255)                                                                                  gtgaGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGGAGCTCTGCGGCAGCTa                                         >  1:1150113/1‑139 (MQ=255)                                                                                     tgGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATg                                      >  2:243349/1‑139 (MQ=255)                                                                                          aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                 <  2:689323/139‑1 (MQ=255)                                                                                          aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                 <  1:136261/139‑1 (MQ=255)                                                                                           gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgc                                >  1:1486933/1‑139 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                               <  1:1557980/139‑1 (MQ=255)                                                                                            tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                               <  1:1542036/139‑1 (MQ=255)                                                                                                      aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                           <  1:268628/53‑1 (MQ=255)                                                                                                      aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                           >  2:268628/1‑53 (MQ=255)                                                                                                      aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAg                                            <  2:1219197/116‑1 (MQ=255)                                                                                                      aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAg                                            >  1:1219197/1‑116 (MQ=255)                                                                                                      aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAata                     >  2:1464462/1‑139 (MQ=255)                                                                                                      aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAata                     >  2:351309/1‑139 (MQ=255)                                                                                                           ttAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                               <  2:544965/124‑1 (MQ=255)                                                                                                           ttAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                               >  1:544965/1‑124 (MQ=255)                                                                                                             aaTCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATATCCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcgc                              <  2:1212799/123‑1 (MQ=255)                                                                                                             aaTCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          <  1:813885/63‑1 (MQ=255)                                                                                                             aaTCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                          >  2:813885/1‑63 (MQ=255)                                                                                                             aaTCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcgc                              >  1:1212799/1‑123 (MQ=255)                                                                                                                  tAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                    >  2:1450329/1‑80 (MQ=255)                                                                                                                  tAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                    <  1:1450329/80‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTGACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACg                                                                                <  2:1080981/61‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTGACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACg                                                                                >  1:1080981/1‑61 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                               <  2:424836/110‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                               >  1:424836/1‑110 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  >  2:1247629/1‑139 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  <  1:243349/139‑1 (MQ=255)                                                                                                                         tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc  <  1:1464462/139‑1 (MQ=255)                                                                                                                          cTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                    >  2:508262/1‑72 (MQ=255)                                                                                                                          cTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                    <  1:508262/72‑1 (MQ=255)                                                                                                                                    |                                                                                                                                GAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGC  >  NC_000913/1794041‑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: ‑