Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,462,102 A→G 9.6% F53L (TTC→CTC)  sixA ← phosphohistidine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,462,1020AG 9.6% 138.7 / 7.1 52F53L (TTC→CTC) sixAphosphohistidine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (20/27);  new base G (2/3);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

ACTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGA  >  AP009048/2462035‑2462168
                                                                   |                                                                  
aCTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACg                                                                  <  1:3913573/70‑1 (MQ=255)
 cTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACggg                                                                <  1:2446052/71‑1 (MQ=255)
       cACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTc                                                             <  1:1470703/68‑1 (MQ=255)
       cACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTc                                                             <  1:2066408/68‑1 (MQ=255)
       cACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCa                                                            >  1:2785513/1‑69 (MQ=255)
       cACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCAc                                                           <  1:133742/70‑1 (MQ=255)
        aCTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCAc                                                           >  1:1915373/1‑69 (MQ=255)
          tGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGaa                                                                    >  1:560297/1‑58 (MQ=255)
          tGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAAc                                                                   <  1:3308792/59‑1 (MQ=255)
            gAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAAc                                                     >  1:2081502/1‑71 (MQ=255)
               ggCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAAcgcg                                                  <  1:2852745/71‑1 (MQ=255)
               ggCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAAcgcg                                                  <  1:1650013/71‑1 (MQ=255)
                gCAGGTTCAGACAATCCCCTACTTCTTCTAGTGTTTGCTCGGATCGCAGGAACGGGCt                                                              >  1:151591/1‑58 (MQ=255)
                gCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAAc                                                                   <  1:2590642/53‑1 (MQ=255)
                   ggTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACg                                                                  <  1:118091/51‑1 (MQ=255)
                    gTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGaa                                                                    >  1:2547885/1‑48 (MQ=255)
                    gTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGa                                             <  1:2508647/71‑1 (MQ=255)
                     ttCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGAt                                            <  1:315102/71‑1 (MQ=255)
                      tCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCa                                                         >  1:2407252/1‑57 (MQ=255)
                      tCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGaa                                                      >  1:4185591/1‑60 (MQ=255)
                      tCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGa                                             >  1:3647027/1‑69 (MQ=255)
                       cAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGAtt                                           >  1:3373261/1‑70 (MQ=255)
                       cAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGAttt                                          >  1:2905200/1‑71 (MQ=255)
                       cAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGAtt                                           <  1:4465886/70‑1 (MQ=255)
                         gACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACg                                                                  <  1:676333/45‑1 (MQ=255)
                          aCAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTc                                         >  1:4270324/1‑69 (MQ=255)
                           cAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCAc                                      >  1:3275331/1‑71 (MQ=255)
                            aaTCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACt                                     <  1:3158449/71‑1 (MQ=255)
                              tCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACttt                                   >  1:353733/1‑71 (MQ=255)
                                cccTACTTCTTCCAGTGTTTGCTCGGCTCTCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACtttt                                  <  1:973434/70‑1 (MQ=255)
                                cccTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACttttt                                 <  1:3467795/71‑1 (MQ=255)
                                  cTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAg                                                        <  1:4084623/46‑1 (MQ=255)
                                   tACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTc                                               >  1:1878604/1‑54 (MQ=255)
                                   tACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGa                               >  1:3005269/1‑70 (MQ=255)
                                     cttcttCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGAcc                             <  1:4332385/70‑1 (MQ=255)
                                      ttcttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTc                                         >  1:4462448/1‑57 (MQ=255)
                                        cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCttt                          >  1:908294/1‑70 (MQ=255)
                                        cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCttt                          <  1:3874238/70‑1 (MQ=255)
                                        cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGATCACCAGAACGCGTTCGATTTc                                         >  1:3774529/1‑55 (MQ=255)
                                           ccAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCttt                          >  1:2330754/1‑67 (MQ=255)
                                           ccAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAg                       <  1:3794744/70‑1 (MQ=255)
                                                gtTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGa                                             >  1:1589157/1‑43 (MQ=255)
                                                gtTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGa                                             >  1:2883792/1‑43 (MQ=255)
                                                       cggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTg                                <  1:2247140/49‑1 (MQ=255)
                                                        ggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCa         <  1:2022500/71‑1 (MQ=255)
                                                          ctcgCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACttttt                                 <  1:1937919/45‑1 (MQ=255)
                                                           tcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCttt                          <  1:1757526/51‑1 (MQ=255)
                                                             gCAGGAACGGGCTCACCAGAACGCGTTCGATTTGCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCga    <  1:3449388/71‑1 (MQ=255)
                                                             gCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGc                      <  1:3576664/53‑1 (MQ=255)
                                                             gCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGc                      <  1:244717/53‑1 (MQ=255)
                                                               aGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGc      >  1:4325581/1‑67 (MQ=255)
                                                               aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCt                            <  1:3652245/45‑1 (MQ=255)
                                                               aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGc      >  1:3582332/1‑67 (MQ=255)
                                                               aGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCgata  <  1:1976507/71‑3 (MQ=255)
                                                                  aaCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCAt           >  1:2017486/1‑59 (MQ=255)
                                                                   |                                                                  
ACTTCTGCACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGA  >  AP009048/2462035‑2462168

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