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

CGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCT  >  NC_000913/1794037‑1794300                                                                                                                                        |                                                                                                                                 cGGCGCACTTCTCCAACGTGGTCACTTCACCGTCGATATCTTTTACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTaa                                                                                                                               >  2:1069381/1‑139 (MQ=255)   gCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACg                                                                                                                             >  1:607447/1‑139 (MQ=255)     gAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTc                                                                                                                           <  1:268371/139‑1 (MQ=255)                 cGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACa                                                                                                               <  1:172786/139‑1 (MQ=255)                     gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTaa                                                                                                                 >  1:1431075/1‑133 (MQ=255)                     gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTaa                                                                                                                 <  2:1431075/133‑1 (MQ=255)                     gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATaa                                                                                                           >  1:954016/1‑139 (MQ=255)                     gTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATaa                                                                                                           >  2:227824/1‑139 (MQ=255)                      tCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                    >  2:1416674/1‑129 (MQ=255)                      tCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCttt                                                                                                                    <  1:1416674/129‑1 (MQ=255)                           tCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGAt                                                                                                     >  2:659457/1‑139 (MQ=255)                                gTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                            <  2:567806/127‑1 (MQ=255)                                gTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATa                                                                                                            >  1:567806/1‑127 (MQ=255)                                       tCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgcgc                                                                                         <  1:227824/139‑1 (MQ=255)                                                 cGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGtt                                                                               >  2:440632/1‑139 (MQ=255)                                                   gtcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                          >  1:59635/1‑110 (MQ=255)                                                   gtcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                          <  2:59635/110‑1 (MQ=255)                                                    tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCg                                                                            <  1:440632/139‑1 (MQ=255)                                                    tcgtcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCg                                                                            >  1:1221916/1‑139 (MQ=255)                                                       tcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                          <  2:378876/106‑1 (MQ=255)                                                       tcAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAAc                                                                                                          >  1:378876/1‑106 (MQ=255)                                                         aGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAc                                                                       <  2:602480/139‑1 (MQ=255)                                                                     tGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAgg                                                           >  1:1190441/1‑139 (MQ=255)                                                                                      gtgtGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTa                                          >  1:353890/1‑139 (MQ=255)                                                                                         tgGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                     >  1:90913/1‑109 (MQ=255)                                                                                         tgGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCa                                                                     <  2:90913/109‑1 (MQ=255)                                                                                          ggTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTg                                                                         <  1:405755/104‑1 (MQ=255)                                                                                          ggTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTg                                                                         >  2:405755/1‑104 (MQ=255)                                                                                           gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                           >  1:1077670/1‑85 (MQ=255)                                                                                           gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                           <  2:1077670/85‑1 (MQ=255)                                                                                           gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCgg                                                <  2:728256/128‑1 (MQ=255)                                                                                           gTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCgg                                                >  1:728256/1‑128 (MQ=255)                                                                                                tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGGTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                <  2:954016/139‑1 (MQ=255)                                                                                                tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATATCCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                <  1:861238/139‑1 (MQ=255)                                                                                                tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                <  2:1013128/139‑1 (MQ=255)                                                                                                tAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                <  1:659457/139‑1 (MQ=255)                                                                                                                      tAGACCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCagcagc          >  1:97244/1‑139 (MQ=255)                                                                                                                         aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCg                   >  1:102567/1‑127 (MQ=255)                                                                                                                         aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCg                   <  2:102567/127‑1 (MQ=255)                                                                                                                         aCCGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAAt       >  1:550395/1‑139 (MQ=255)                                                                                                                           cGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCTCCTAACGTTGCGGTGAc                                                                       <  2:819123/73‑1 (MQ=255)                                                                                                                           cGTCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAc                                                                       >  1:819123/1‑73 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                           >  2:542893/1‑51 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGgc                                                                                           <  1:542893/51‑1 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAc                                                                       <  2:171913/71‑1 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGAc                                                                       >  1:171913/1‑71 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                >  1:655261/1‑110 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGgcg                                <  2:655261/110‑1 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATgcgc   >  2:1128075/1‑139 (MQ=255)                                                                                                                             tCTTTTTGCGGTTAACGTCCCGCTTTTAACAATAAc                                                                                                          <  2:800199/36‑1 (MQ=37)                                                                                                                              cTTTTTGCGGTTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCt  >  1:489933/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                 CGGCGAACTTCTCCAGCGTGGTCACTTCACCGTCGATATCTTTCACTACGGTCGTCAGAATGGAATCTTGCATCGTTTCTCCTGGGTGTGGTCAGTAAAAATCTTAGCTTTTAATCATAGACCGTCTTTTTGCGGCTAACGTCCTGCTTTTAACAATAACCAGATAAACACCGGCGCACCTAACGTTGCGGTGACCACGCCAATAGGCAGCTCTGCGGCAGCTAATGCCAGGCGCGCTACAATATCGGCCAGCAGCAATGCGCT  >  NC_000913/1794037‑1794300

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