Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,683,257 G→A 17.1% D190N (GAT→AAT)  folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,683,2570GA17.1% 82.0 / 10.2 35D190N (GAT→AAT) folMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base A (2/4);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

GCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAG  >  AP009048/1683188‑1683323
                                                                     |                                                                  
gCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatg                                                                    <  1:1885887/70‑1 (MQ=255)
 cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  <  1:857023/71‑1 (MQ=255)
 cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  <  1:2608846/71‑1 (MQ=255)
  ccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatg                                                                 <  1:4580393/71‑1 (MQ=255)
   cGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatg                                                                 <  1:3006128/70‑1 (MQ=255)
    gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat                                                               <  1:4427183/71‑1 (MQ=255)
    gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat                                                               >  1:331470/1‑71 (MQ=255)
     cAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgatg                                                              >  1:2551219/1‑71 (MQ=255)
           ggCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGaa                                                         >  1:4109880/1‑70 (MQ=255)
            gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAata                                                       <  1:186745/71‑1 (MQ=255)
                       gaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAgg                                            >  1:1552299/1‑71 (MQ=255)
                        agagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGgc                                           >  1:2865452/3‑71 (MQ=255)
                         aagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga                                                                >  1:1595426/1‑49 (MQ=255)
                          agtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga                                                                >  1:369516/1‑48 (MQ=255)
                            tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg                                                           >  1:715910/1‑51 (MQ=255)
                             gaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACataat                                                                  >  1:2005845/1‑43 (MQ=255)
                                 tCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATa                                   >  1:1058434/1‑70 (MQ=255)
                                 tCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATa                                   >  1:1808092/1‑70 (MQ=255)
                                  cTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATaaa                                 >  1:1777249/1‑71 (MQ=255)
                                  cTATTGCGCCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGa                                                   >  1:1709589/1‑53 (MQ=255)
                                   tATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga                                                                <  1:399819/39‑1 (MQ=255)
                                   tATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATaa                                  <  1:3914995/69‑1 (MQ=255)
                                     ttGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAat                                                        <  1:4138624/45‑1 (MQ=255)
                                     ttGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATc                               >  1:2707383/1‑70 (MQ=255)
                                        cgcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAata                                                       <  1:4125130/43‑1 (MQ=255)
                                        cgcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACt                            <  1:3280493/70‑1 (MQ=255)
                                        cgcCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           <  1:2433246/71‑1 (MQ=255)
                                          ccATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgat                         <  1:444256/71‑1 (MQ=255)
                                            aTCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatga                       <  1:1356547/71‑1 (MQ=255)
                                            aTCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatga                       <  1:246031/71‑1 (MQ=255)
                                             tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatga                       >  1:734101/1‑70 (MQ=255)
                                             tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:2038623/71‑1 (MQ=255)
                                             tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:4316036/71‑1 (MQ=255)
                                              cGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATaaa                                 >  1:3077675/1‑59 (MQ=255)
                                               gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaa                     <  1:200390/70‑1 (MQ=255)
                                                   aTCCTGTTTAATGAACATGATGATGCCGAATATCGa                                                   <  1:3835566/36‑1 (MQ=255)
                                                     ccTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAAt                                <  1:3949745/53‑1 (MQ=255)
                                                     ccTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCgcgc              >  1:856424/1‑71 (MQ=255)
                                                         tttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATa                                   <  1:1400436/46‑1 (MQ=255)
                                                         tttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTg           <  1:1034170/70‑1 (MQ=255)
                                                          ttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCgcgc              >  1:1616250/1‑66 (MQ=255)
                                                             aTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCgag      <  1:27150/71‑1 (MQ=255)
                                                              tGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCgaga     <  1:2875475/71‑1 (MQ=255)
                                                                 aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCga       >  1:2383279/1‑66 (MQ=255)
                                                                 aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaag  <  1:2094050/71‑1 (MQ=255)
                                                                 aCATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaag  <  1:3429331/71‑1 (MQ=255)
                                                                  cATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaag  >  1:1603077/1‑70 (MQ=255)
                                                                    tgatgatgCCGAATATCGACAACAGGCGCTGAATa                                   <  1:3339483/35‑1 (MQ=255)
                                                                     gatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCt            <  1:2088166/57‑1 (MQ=255)
                                                                     |                                                                  
GCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAG  >  AP009048/1683188‑1683323

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