Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	796,704	C→T	Q51* (CAA→TAA)	modC →	molybdate ABC transporter ATPase; chlorate resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	796,704	0	C	T	100.0%	172.3 / NA	50	Q51* (CAA→TAA)	modC	molybdate ABC transporter ATPase; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/26);  total (24/26)

CGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATG  >  NC_000913/796612‑796794                                                                                             |                                                                                           cGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAgg                                                                                     >  1:350303/1‑100 (MQ=255) cGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAg                                                                                      <  2:752250/99‑1 (MQ=255)   cTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACAGCAGTGGACTGACGCGCCCGTAAAAAGGGc                                                                                   <  2:494693/100‑1 (MQ=255)      cttccgatctgCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGa                                                                                <  2:633008/90‑1 (MQ=255)      cccGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGa                                                                                >  1:780185/1‑100 (MQ=255)              tGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCTGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGAATGTCCTc                                                                        <  2:684884/100‑1 (MQ=255)                gCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGaagatcggaag                                                                      >  1:633012/1‑90 (MQ=255)                gCATCACTGCTATCTTTGGCGTCTCCGGCGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCaa                                                                      >  2:127793/1‑100 (MQ=255)                  aTCACTGCTATCTTTGGCGTCGCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATg                                                                    <  1:1035747/100‑1 (MQ=255)                    cACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATggg                                                                  >  1:1038839/1‑100 (MQ=255)                        gCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCggg                                                              <  1:3050/100‑1 (MQ=255)                            tcgttGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGCACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACt                                                          <  2:491808/97‑1 (MQ=255)                             cTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTa                                                         >  1:1185172/1‑100 (MQ=255)                                tGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAat                                                      <  1:291781/100‑1 (MQ=255)                                tGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGCACTAAat                                                      >  2:1288128/1‑100 (MQ=255)                                 ggCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatg                                                     <  2:158487/100‑1 (MQ=255)                                  gCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatga                                                    <  2:541021/100‑1 (MQ=255)                                  gCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatga                                                    <  1:531560/100‑1 (MQ=255)                                   cGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatgat                                                   >  1:738738/1‑100 (MQ=255)                                    gTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatgatg                                                  >  1:1119581/1‑100 (MQ=255)                                     tctcCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGc                                                 <  1:1004792/100‑1 (MQ=255)                                      ctcCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGcc                                                <  2:463990/100‑1 (MQ=255)                                       tcCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCg                                               >  2:1232715/1‑100 (MQ=255)                                        ccGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGa                                              <  1:803424/100‑1 (MQ=255)                                             gCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAg                                         >  1:997213/1‑100 (MQ=255)                                              ccGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAgg                                        <  1:49969/100‑1 (MQ=255)                                              ccGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAgg                                        <  1:1166815/100‑1 (MQ=255)                                              ccGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAgg                                        <  2:1001041/100‑1 (MQ=255)                                                ggAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTa                                      >  1:105696/1‑100 (MQ=255)                                                      aCTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGcc                                >  1:699812/1‑100 (MQ=255)                                                        ttCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTg                              <  2:699808/100‑1 (MQ=255)                                                        ttCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATAAAGAAGAAAAAGGTATCTGCCTg                              <  2:997204/100‑1 (MQ=255)                                                         tCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGa                             >  2:484123/1‑100 (MQ=255)                                                            cTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGc                          <  2:400769/100‑1 (MQ=255)                                                             tGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGcc                         >  2:1270894/1‑100 (MQ=255)                                                             tGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGcc                         <  2:10674/100‑1 (MQ=255)                                                             tGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGcc                         >  1:720743/1‑100 (MQ=255)                                                             tGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGcc                         <  1:732424/100‑1 (MQ=255)                                                              gATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg                        >  2:198897/1‑100 (MQ=255)                                                              gATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg                        >  1:471893/1‑100 (MQ=255)                                                                  aacgagaacAGTGGACTGACGCGCCCGTAAAAAGTGCGGATTGTTTTCAATGGGCGGGTAATAAAAAAAAAAAAAAAAGGTATCTGCCTGACGcc                         <  2:1151767/87‑1 (MQ=255)                                                                  aaCGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGaaa                    >  1:482072/1‑100 (MQ=255)                                                                   aCGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGaaaa                   >  2:1284259/1‑100 (MQ=255)                                                                         tgaGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTcg             >  1:1205605/3‑100 (MQ=255)                                                                         tCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTcg             >  2:402432/1‑100 (MQ=255)                                                                           aGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTcgcg           >  2:372292/1‑100 (MQ=255)                                                                                aCTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGc      <  1:688357/100‑1 (MQ=255)                                                                                aCTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGc      <  1:447766/100‑1 (MQ=255)                                                                                   gACCCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAAAGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTa    <  2:102197/99‑1 (MQ=255)                                                                                    aCGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATg  >  1:1024271/1‑100 (MQ=255)                                                                                             |                                                                                           CGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCGTCGCGTTGGCTATG  >  NC_000913/796612‑796794

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