Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,311,972 A→T 7.3% L169Q (CTG→CAG)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,311,9720AT 7.3% 146.4 / 2.1 53L169Q (CTG→CAG) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base T (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGC  >  AP009048/4311902‑4312042
                                                                      |                                                                      
cgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCa                                                                        <  1:6010360/71‑1 (MQ=255)
cgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCa                                                                        <  1:3342307/71‑1 (MQ=255)
  ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                      <  1:5632194/71‑1 (MQ=255)
     ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACTCGCAGTGTTGGGTCATATCTCCCAgggg                                                                    >  1:7138180/1‑70 (MQ=255)
       ccAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                  <  1:184821/70‑1 (MQ=255)
       ccAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                  >  1:4293080/1‑70 (MQ=255)
           agCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                             <  1:6280284/71‑1 (MQ=255)
            gCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                             <  1:7297246/70‑1 (MQ=255)
               aTTGGTTTCCAGGCACCCAGGATTGCCACACGTGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                               <  1:4253313/65‑1 (MQ=255)
               aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                      <  1:5939205/58‑1 (MQ=255)
               aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                           >  1:5825961/1‑69 (MQ=255)
                  ggTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                       >  1:5219794/1‑54 (MQ=255)
                  ggTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                       >  1:2603593/1‑70 (MQ=255)
                  ggTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                       >  1:7829215/1‑70 (MQ=255)
                   gTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGtt                                                      <  1:4838508/70‑1 (MQ=255)
                      tCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                  >  1:6118296/1‑71 (MQ=255)
                       ccAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                  <  1:279144/54‑1 (MQ=255)
                        cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGata                                                                >  1:4507095/1‑55 (MQ=25)
                        cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGata                                                                >  1:3711675/1‑55 (MQ=25)
                         aGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                               <  1:7683898/71‑1 (MQ=255)
                          ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                       >  1:498656/1‑46 (MQ=255)
                          ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                       >  1:44297/1‑46 (MQ=255)
                          ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                       >  1:1196665/1‑46 (MQ=255)
                           gCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                             >  1:7637531/1‑55 (MQ=255)
                           gCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCca                                             >  1:5636280/1‑71 (MQ=255)
                              cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                       <  1:1905303/58‑1 (MQ=255)
                              cccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                       <  1:2556360/58‑1 (MQ=255)
                               ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                               <  1:3884894/49‑1 (MQ=255)
                               ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                               <  1:5362441/49‑1 (MQ=255)
                               ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                               <  1:7029045/49‑1 (MQ=255)
                                cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                      <  1:2029082/41‑1 (MQ=255)
                                cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCg                                        <  1:2491336/71‑1 (MQ=255)
                                cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCg                                        <  1:3964394/71‑1 (MQ=255)
                                 aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                       <  1:3027843/71‑1 (MQ=255)
                                 aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                       <  1:5061556/71‑1 (MQ=255)
                                  ggATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                      <  1:1288495/71‑1 (MQ=255)
                                   gATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                     <  1:325149/71‑1 (MQ=255)
                                   gATTGACACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACCCCgtgt                                     >  1:2231038/1‑71 (MQ=255)
                                           cccgcgCAGTGTTGTGTCATATCTCCCTGGGTGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGt                             <  1:1675707/69‑1 (MQ=255)
                                                    tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                    <  1:4669530/71‑1 (MQ=255)
                                                     gtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCa                   >  1:1897395/1‑71 (MQ=255)
                                                           tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGc                                   <  1:6773618/49‑1 (MQ=255)
                                                             atatCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCCCACCGTGTGCACCCGTCCACGGCGCACCGTTCa            >  1:2134283/1‑70 (MQ=255)
                                                             atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa            <  1:1622652/70‑1 (MQ=255)
                                                             atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt           >  1:7702335/1‑71 (MQ=255)
                                                               atCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCgtg                                      <  1:1005745/42‑1 (MQ=18)
                                                                tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt           >  1:7939607/1‑68 (MQ=255)
                                                                tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt           >  1:7907842/1‑68 (MQ=255)
                                                                tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt           >  1:5975534/1‑68 (MQ=255)
                                                                 ctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAc                  >  1:3162310/1‑60 (MQ=255)
                                                                   cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAc                                 >  1:4009314/1‑43 (MQ=255)
                                                                   cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAc                                 >  1:4426371/1‑43 (MQ=255)
                                                                   cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAc                                 >  1:6215562/1‑43 (MQ=255)
                                                                   cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt           >  1:7878764/1‑65 (MQ=255)
                                                                    ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACt    >  1:499503/1‑71 (MQ=255)
                                                                      aGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGc  >  1:7601632/1‑71 (MQ=255)
                                                                      aGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGc  >  1:6489896/1‑71 (MQ=255)
                                                                      |                                                                      
CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGC  >  AP009048/4311902‑4312042

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: