Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 797,155 T→C V201A (GTA→GCA)  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_000913797,1550TC93.5% 99.8 / ‑1.0 31V201A (GTA→GCA) modCmolybdate ABC transporter ATPase; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (13/16);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GCTGGCGTCACTGGATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGTACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAGCATCATCCGCATTACGCGA  >  NC_000913/797018‑797295
                                                                                                                                         |                                                                                                                                            
gCTGGCGTCACTGGATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGt                                                                                                                                 >  1:5700/1‑151 (MQ=255)
   ggCGTCACTGGATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTc                                                                                                                                >  2:510505/1‑149 (MQ=255)
      gTCACTGGATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTg                                                                                                                           <  2:3150/151‑1 (MQ=255)
      gTCACTGGATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTg                                                                                                                           >  1:280399/1‑151 (MQ=255)
            ggATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGcc                                                                                                                     >  2:357398/1‑151 (MQ=255)
                atTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGTACTGGAAAACGGTCAGGTGAAAGCCttt                                                                                                                  <  2:23508/150‑1 (MQ=255)
                           aCGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTgg                                                                                                        >  1:588762/1‑149 (MQ=255)
                             gcgAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTggagg                                                                                                     >  1:132684/1‑150 (MQ=255)
                               gAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAg                                                                                                  <  2:379090/151‑1 (MQ=255)
                                aaCTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAgt                                                                                                 >  2:97439/1‑151 (MQ=255)
                                           tATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGt                                                                                       <  1:481105/150‑1 (MQ=255)
                                            aTCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTa                                                                                      <  2:294352/150‑1 (MQ=255)
                                                         gACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATcc                                                                         <  1:452314/150‑1 (MQ=255)
                                                          acacGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGt                                                                       >  2:445119/1‑151 (MQ=255)
                                                           cacGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTg                                                                      <  1:217124/151‑1 (MQ=255)
                                                              gggAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCt                                                                   <  1:97439/151‑1 (MQ=255)
                                                              gggAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCt                                                                   >  1:387623/1‑151 (MQ=255)
                                                                   aTCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGa                                                              <  1:183195/151‑1 (MQ=255)
                                                                        cATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAagag                                                         <  2:247784/151‑1 (MQ=255)
                                                                             cGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGcaaca                                                    <  1:133606/151‑1 (MQ=255)
                                                                                      aTGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGTACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCAt                                           >  2:544887/1‑151 (MQ=255)
                                                                                          cAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCt                                        <  1:548962/150‑1 (MQ=255)
                                                                                                  cGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGAc                               >  2:27061/1‑151 (MQ=255)
                                                                                                   gCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACg                              <  2:473494/151‑1 (MQ=255)
                                                                                                        aTGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCt                         <  1:452293/151‑1 (MQ=255)
                                                                                                         tGAGATCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCtt                        >  2:26473/1‑151 (MQ=255)
                                                                                                            gaTCCTCCATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAg                     <  2:19034/151‑1 (MQ=255)
                                                                                                               cctccATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAGcat                  >  1:321677/1‑151 (MQ=255)
                                                                                                                 tccATCTGGCAGACAGAGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAGcatc                 >  2:546387/1‑150 (MQ=255)
                                                                                                                             acagaGTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAGCATCATCCGCATTAcgc    <  2:68532/151‑1 (MQ=255)
                                                                                                                               agagTGATGGCACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAGCATCATCCGCATTACGCGa  <  1:164556/151‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
GCTGGCGTCACTGGATATTCCGCGTAAACGCGAACTGTTGCCTTATCTGCAACGGCTGACACGGGAAATCAACATTCCGATGTTGTATGTCAGCCATTCGCTGGATGAGATCCTCCATCTGGCAGACAGAGTGATGGTACTGGAAAACGGTCAGGTGAAAGCCTTTGGCGCGCTGGAGGAAGTGTGGGGCAGTAGCGTGATGAATCCGTGGCTGCCGAAAGAGCAACAAAGTAGCATTCTGAAAGTGACGGTGCTTGAGCATCATCCGCATTACGCGA  >  NC_000913/797018‑797295

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