Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,722,719	C→T	G13D (GGT→GAT)	cobW3 ←	cobalamin (vitamin B12) biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,722,719	0	C	T	100.0%	117.1 / NA	34	G13D (GGT→GAT)	cobW3	cobalamin (vitamin B12) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

AGTTCACTTAAAGGTATGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTACCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGATCTCCTATTGACTGACATGACTCTATATTCAGTTATAGTT  >  CP000730/2722584‑2722845                                                                                                                                        |                                                                                                                               aGTTCACTTAAAGGTATGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATCATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGACCAAATGATTTAAGACAGTCGTTTTATCGcta                                                                                                                           >  2:138526/1‑139 (MQ=255)       cTTAAAGGTATGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTaaa                                                                                                                     >  1:249420/1‑141 (MQ=255)                 tGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCt                                                                                                                        <  2:12027/128‑1 (MQ=255)                 tGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCt                                                                                                                        <  2:11982/128‑1 (MQ=255)                 tGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCt                                                                                                                        >  1:11982/1‑128 (MQ=255)                 tGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCt                                                                                                                        >  1:12027/1‑128 (MQ=255)                 tGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTaaaaa                                                                                                                   <  1:27736/133‑1 (MQ=255)                 tGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTaaaaa                                                                                                                   >  2:27736/1‑133 (MQ=255)                   tCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCtata                                                                                                         <  1:183963/141‑1 (MQ=255)                    cTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTataa                                                                                                        <  2:100647/141‑1 (MQ=255)                        aGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCtata                                                                                                         <  1:165231/136‑1 (MQ=255)                        aGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCtata                                                                                                         >  2:165231/1‑136 (MQ=255)                        aGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAAc                                                                                                    >  2:119618/1‑141 (MQ=255)                            aCTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTAtt                                                                                                >  1:7100/1‑141 (MQ=255)                              tAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTAtttt                                                                                              >  2:56015/1‑141 (MQ=255)                                  gCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATc                                                                                          >  2:81388/1‑141 (MQ=255)                                        tCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGAc                                                                                    >  2:91339/1‑141 (MQ=255)                                            cacTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACtctc                                                                                >  2:249076/1‑141 (MQ=255)                                                  ttttCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTta                                                                          <  1:56015/141‑1 (MQ=255)                                                   tttCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTtat                                                                         <  1:81388/141‑1 (MQ=255)                                                        cAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTc                                                                    <  2:161493/141‑1 (MQ=255)                                                                 tcatGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATa                                                           >  1:126505/1‑141 (MQ=255)                                                                   atgatgACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGc                                                         <  1:119618/141‑1 (MQ=255)                                                                         aCTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTc                                                   <  2:186843/141‑1 (MQ=255)                                                                             ctgGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAag                                               <  1:91409/141‑1 (MQ=255)                                                                             ctgGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAag                                               <  1:136773/141‑1 (MQ=255)                                                                             ctgGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAag                                               <  1:179531/141‑1 (MQ=255)                                                                                     aGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGAtct                                       <  1:105213/141‑1 (MQ=255)                                                                                      gCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGAtctc                                      >  2:120754/1‑141 (MQ=255)                                                                                          ttCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGcc                                                                                                                            <  2:109808/51‑1 (MQ=255)                                                                                          ttCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTATCGcc                                                                                                                            >  1:109808/1‑51 (MQ=255)                                                                                                          gCGAGCAAATGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGATCTCCTATTGACTGTCATGActct                  <  2:132852/141‑1 (MQ=255)                                                                                                                aaaTGATTTAAGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGATCTCCTATTGACTGACATGACTCTATATTc            >  1:204220/1‑141 (MQ=255)                                                                                                                          aGACAGTCGTTTTATCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGATCTCCTATTGACTGACATGACTCTATATTCAGTTATAGtt  >  1:278024/1‑141 (MQ=255)                                                                                                                                        |                                                                                                                               AGTTCACTTAAAGGTATGTCTTCAGATACTAAGGCACCATCAACACTCATTTTCCCAAATTCATTCATGATGACTGCTGGTTTCAGCGATTCCTTTAATGATTCAGCGAGCAAATGATTTAAGACAGTCGTTTTACCGCCACCTAAAAACCCACCTATAATAACTATTTTCATCTGCGACTCTCCTCTTATACTTCATATGAATAGCATATTCCAAGATAGATCTCCTATTGACTGACATGACTCTATATTCAGTTATAGTT  >  CP000730/2722584‑2722845

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