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

GCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTC  >  NC_000913/1794039‑1794301                                                                                                                                      |                                                                                                                                  gCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACg                                                                                                                              >  2:526524/1‑139 (MQ=255)  cGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGt                                                                                                                             <  2:1510188/139‑1 (MQ=255)   gAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            <  2:1571630/139‑1 (MQ=255)    aaCTTCTCCAGCGTGGTCACTTCACCGTTGATATATTTCACTACGGTCGGCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTcc                                                                                                                           >  2:1070270/1‑139 (MQ=255)          tcCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                     >  2:174244/1‑139 (MQ=255)                   gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATGATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATaa                                                                                                            >  1:1467378/1‑139 (MQ=255)                   gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATaa                                                                                                            >  2:401099/1‑139 (MQ=255)                    tCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                           <  1:401099/139‑1 (MQ=255)                                               cGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            <  1:87790/95‑1 (MQ=255)                                               cGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            >  2:87790/1‑95 (MQ=255)                                                ggTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTg                                                                               >  1:869086/1‑139 (MQ=255)                                                 gtcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            <  1:453363/93‑1 (MQ=255)                                                 gtcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            >  2:453363/1‑93 (MQ=255)                                                  tcgtcAGAATGGAATCTTGCATTGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCg                                                                             <  2:577849/139‑1 (MQ=255)                                                  tcgtcAGAATGGAATCTTGCATTGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCg                                                                             >  1:577849/1‑139 (MQ=255)                                                  tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTg                                                                                                                         >  2:1401209/1‑95 (MQ=255)                                                  tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTg                                                                                                                         <  1:1401209/95‑1 (MQ=255)                                                    gtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCg                                                                                              <  2:862915/120‑1 (MQ=255)                                                    gtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCg                                                                                              >  1:862915/1‑120 (MQ=255)                                                     tcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTg                                                                          >  1:162499/1‑139 (MQ=255)                                                         aaTGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                      <  2:1248829/139‑1 (MQ=255)                                                                tCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            >  2:1025572/1‑78 (MQ=255)                                                                tCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                            <  1:1025572/78‑1 (MQ=255)                                                                   tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAgg                                                            >  1:1193881/1‑139 (MQ=255)                                                                   tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAgg                                                            >  2:679915/1‑139 (MQ=255)                                                                          tttCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGctct                                                     >  1:1079484/1‑139 (MQ=255)                                                                              tcCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACGATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGTTCTGCgg                                                 <  2:1079484/139‑1 (MQ=255)                                                                                  gggTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAcc                                                                                                          >  1:1119400/1‑78 (MQ=255)                                                                                  gggTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAcc                                                                                                          <  2:1119400/78‑1 (MQ=255)                                                                                   ggTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGTAGCTCTGCGGCAGCt                                            >  1:875086/1‑139 (MQ=255)                                                                                       tgGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATg                                        <  1:526524/139‑1 (MQ=255)                                                                                         gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                           >  1:817904/1‑70 (MQ=255)                                                                                         gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                           <  2:817904/70‑1 (MQ=255)                                                                                         gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGcc                                      >  2:1483387/1‑139 (MQ=255)                                                                                            aGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAgg                                   <  1:1568514/139‑1 (MQ=255)                                                                                             gTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACC‑GATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 >  1:764486/1‑139 (MQ=255)                                                                                              tAAAAATCTTAGCTTTTAATGATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 <  2:1467378/139‑1 (MQ=255)                                                                                              tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                 <  2:522368/139‑1 (MQ=255)                                                                                                        aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAcc                                                                                                          >  1:210597/1‑56 (MQ=255)                                                                                                        aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAcc                                                                                                          <  2:210597/56‑1 (MQ=255)                                                                                                        aGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAata                       >  2:80023/1‑139 (MQ=255)                                                                                                         gCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGTAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAatat                      <  2:875086/139‑1 (MQ=255)                                                                                                                     aGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAgcagca          <  1:685850/139‑1 (MQ=255)                                                                                                                       aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt        >  1:1391289/1‑139 (MQ=255)                                                                                                                           tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                            >  1:1000283/1‑51 (MQ=255)                                                                                                                           tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                            <  2:1000283/51‑1 (MQ=255)                                                                                                                           tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc    >  1:1303266/1‑139 (MQ=255)                                                                                                                           tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc    <  1:174244/139‑1 (MQ=255)                                                                                                                             tttttGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTc  >  1:54831/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                  GCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCTC  >  NC_000913/1794039‑1794301

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