Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,424 T→C 19.3% I150T (ATT→ACT)  folM → dihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,4240TC19.3% 134.0 / 28.2 68I150T (ATT→ACT) folMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (32/22);  new base C (7/6);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCT  >  NC_000913/1681325‑1681523
                                                                                                   |                                                                                                   
aCCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACatat                                                                                                     <  2:2256704/100‑1 (MQ=255)
  cATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTg                                                                                                   >  1:2465932/1‑100 (MQ=255)
    tGCGCTTGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCg                                                                                                 <  1:668489/100‑1 (MQ=255)
    tGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCg                                                                                                 >  1:108359/1‑100 (MQ=255)
    tGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCg                                                                                                 <  1:1911243/100‑1 (MQ=255)
       gcTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAt                                                                                              >  2:2270565/1‑100 (MQ=255)
       gcTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAt                                                                                              >  2:1957730/1‑100 (MQ=255)
       gcTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAt                                                                                              <  2:1550519/100‑1 (MQ=255)
       gatGGAAAGATTAATGCGTGGGCACGGACACGCCGCCAGCGAGATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAt                                                                                              >  2:1753271/3‑100 (MQ=255)
         tGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTAtgc                                                                                            >  1:2082500/1‑100 (MQ=255)
            aaaGATTACTGCGGGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAtgctgc                                                                                         >  2:1760215/1‑100 (MQ=255)
            aaaGATTACTGAGTGGGCACGGTCACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTAtgct                                                                                           <  1:1420827/98‑1 (MQ=255)
             aaGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTgca                                                                                        <  1:1147926/100‑1 (MQ=255)
             aaGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTgca                                                                                        <  1:225545/100‑1 (MQ=255)
               gATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTgcaag                                                                                      >  2:1881887/1‑100 (MQ=255)
               gATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTgcaag                                                                                      <  1:2402611/100‑1 (MQ=255)
                aTTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTgcaagc                                                                                     <  2:2547673/100‑1 (MQ=255)
                    cTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAg                                                                                 >  1:1393188/1‑100 (MQ=255)
                    cTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAg                                                                                 >  1:1440609/1‑100 (MQ=255)
                    cTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTATGCTGCAAGCAAAg                                                                                 >  1:2128299/1‑100 (MQ=255)
                       cGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCgg                                                                              <  1:944672/100‑1 (MQ=255)
                        gTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGc                                                                             <  2:1909608/100‑1 (MQ=255)
                         tGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                            >  2:1515700/1‑100 (MQ=255)
                         tGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                            >  1:2540204/1‑100 (MQ=255)
                          gggCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGCGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCa                                                                            >  1:2028480/1‑99 (MQ=255)
                           ggCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGGGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACt                                                                          >  2:642216/1‑100 (MQ=255)
                            gCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTg                                                                         >  1:831650/1‑100 (MQ=255)
                              aCGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGa                                                                       >  1:1793562/1‑100 (MQ=255)
                                  acacGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGataat                                                                   >  1:2585711/1‑100 (MQ=255)
                                  acacGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGataat                                                                   <  2:1018211/100‑1 (MQ=255)
                                   cacGCCGCCATCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGataata                                                                  <  1:1694697/100‑1 (MQ=255)
                                   cacGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGataata                                                                  <  2:59448/100‑1 (MQ=255)
                                    acGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAatat                                                                 <  1:1961930/100‑1 (MQ=255)
                                     cgccgcCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATg                                                                >  1:406873/1‑100 (MQ=255)
                                     cgccgcCAGCGATATCATGCACTTTACCGAATATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAatat                                                                 <  2:2152376/99‑1 (MQ=255)
                                       ccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGAc                                                              <  2:267265/100‑1 (MQ=255)
                                        cgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGAcc                                                             <  2:2270622/100‑1 (MQ=255)
                                        cgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGAAGCGACAAACATATTGCGTATGCTGCAAGCAACGCGGCACTGGATAATATGAc                                                              >  1:1338074/1‑99 (MQ=255)
                                                  atCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTg                                                   <  1:839508/100‑1 (MQ=255)
                                                   tCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGc                                                  >  1:1422418/1‑100 (MQ=255)
                                                    cATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGcc                                                 >  2:921060/1‑100 (MQ=255)
                                                    cATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGcc                                                 >  2:2275870/1‑100 (MQ=255)
                                                    cATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGcc                                                 <  1:848815/100‑1 (MQ=255)
                                                       tCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGc                                              >  2:1296732/1‑100 (MQ=255)
                                                        cACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCCCTGGATAATATGACCCGCTCGTTTGCCCGCa                                             >  2:480908/1‑100 (MQ=255)
                                                            ttACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCt                                         >  2:1979377/1‑100 (MQ=255)
                                                               ccGATTATGTGGTGGAGCGCGGTAGCGTCAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGc                                      <  1:1785715/100‑1 (MQ=255)
                                                               acGATTATCCGGTGGCCCCCGGTATCGACAAACATATTGAGAATGCCGCAAGCAAAGCGGGACTGGATAATATGACCCGCTCGTTTGCCCGCAAGATg                                        >  1:390877/2‑98 (MQ=255)
                                                                cGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCa                                     >  2:1327396/1‑100 (MQ=255)
                                                                 gATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCAc                                    >  2:1221079/1‑100 (MQ=255)
                                                                  aTTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCAcc                                   <  1:2306641/100‑1 (MQ=255)
                                                                   ttccgatctggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCAc                                    <  2:506888/90‑1 (MQ=255)
                                                                      tgtgGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaa                               >  1:146022/1‑100 (MQ=255)
                                                                        tggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagt                             >  1:1473205/1‑100 (MQ=255)
                                                                           tggAGCGCGGTAGCGACAAACATAGTGCGTATGCTGCAAGCAAAGCGGCACTGGATAAGATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaa                          <  1:2009030/100‑1 (MQ=255)
                                                                           tggAGCGCGGTAGCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaa                          >  2:1134217/1‑100 (MQ=255)
                                                                            ggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaa                         >  1:1114900/1‑100 (MQ=255)
                                                                             gAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaag                        <  2:2207493/100‑1 (MQ=255)
                                                                             gAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaag                        >  1:2138084/1‑100 (MQ=255)
                                                                                  cgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAt                   >  2:1599101/1‑100 (MQ=255)
                                                                                  cgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAt                   >  2:238610/1‑100 (MQ=255)
                                                                                  cgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAt                   <  2:1037559/100‑1 (MQ=255)
                                                                                       gCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTAt              >  2:2106111/1‑100 (MQ=255)
                                                                                       gCGACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTAt              <  1:774353/100‑1 (MQ=255)
                                                                                         gACAAACATACTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTg            >  2:531395/1‑100 (MQ=255)
                                                                                          aCAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCAGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgc           >  1:19335/1‑100 (MQ=255)
                                                                                            aaaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcgc         <  1:901883/100‑1 (MQ=255)
                                                                                            aaaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcgc         <  2:495276/100‑1 (MQ=255)
                                                                                             aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGcc        <  1:1361579/100‑1 (MQ=255)
                                                                                                  atTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGc   <  2:1128895/100‑1 (MQ=255)
                                                                                                   ttGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCt  <  2:1017843/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ACCATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCT  >  NC_000913/1681325‑1681523

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