Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 796,681 C→A A43D (GCC→GAC)  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_000913796,6810CA98.2% 184.8 / ‑3.4 55A43D (GCC→GAC) modCmolybdate ABC transporter ATPase; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (28/26);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

GCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAA  >  NC_000913/796585‑796777
                                                                                                |                                                                                                
gCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACAt                                                                                               >  1:1265745/1‑100 (MQ=255)
    ccATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGt                                                                                           >  1:541025/1‑100 (MQ=255)
    ccATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGt                                                                                           >  2:732419/1‑100 (MQ=255)
    ccATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGt                                                                                           >  1:84505/1‑100 (MQ=255)
     cATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTg                                                                                          <  1:433053/100‑1 (MQ=255)
        tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGAc                                                                                       <  1:1253851/100‑1 (MQ=255)
        tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGAc                                                                                       <  1:297893/100‑1 (MQ=255)
          ccTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTg                                                                                     <  2:457675/100‑1 (MQ=255)
                 aTTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGcc                                                                              >  1:158489/1‑100 (MQ=255)
                  ttAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGccc                                                                             <  2:272124/100‑1 (MQ=255)
                    aaTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCgt                                                                           >  2:667342/1‑99 (MQ=255)
                     aTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCgta                                                                          >  1:1186929/1‑98 (MQ=255)
                     aTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCgta                                                                          >  1:517399/1‑98 (MQ=255)
                           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                                   <  1:531560/100‑1 (MQ=255)
                                                            gCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatga                                   <  2:541021/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                       <  2:1001041/100‑1 (MQ=255)
                                                                        ccGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAgg                       <  1:1166815/100‑1 (MQ=255)
                                                                        ccGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAgg                       <  1:49969/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:732424/100‑1 (MQ=255)
                                                                                       tGATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGcc        >  1:720743/1‑100 (MQ=255)
                                                                                        gATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg       >  2:198897/1‑100 (MQ=255)
                                                                                        gATTAACGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg       >  1:471893/1‑100 (MQ=255)
                                                                                            aaCGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGaaa   >  1:482072/1‑100 (MQ=255)
                                                                                             aCGACATCAGTGGACTGACGCGCCCGTAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGaaaa  >  2:1284259/1‑100 (MQ=255)
                                                                                                |                                                                                                
GCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAA  >  NC_000913/796585‑796777

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: