Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,264,702 G→A 6.6% R15C (CGC→TGC)  fruA ← fused fructose‑specific PTS enzyme IIB'BC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,264,7020GA 6.6% 238.4 / 4.1 75R15C (CGC→TGC) fruAfused fructose‑specific PTS enzyme IIB'BC components
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base A (2/3);  total (37/38)
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

GATGATTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTAAAA  >  AP009048/2264640‑2264770
                                                              |                                                                    
gatgatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGa                                                               <  1:4141775/70‑1 (MQ=255)
gatgatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGAc                                                              >  1:3653201/1‑71 (MQ=255)
gatgatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGAc                                                              <  1:7439794/71‑1 (MQ=255)
 atgatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGAcc                                                             >  1:4492506/1‑71 (MQ=255)
  tgatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGAcc                                                             >  1:6100317/1‑70 (MQ=255)
   gatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGAcc                                                             >  1:5902456/1‑69 (MQ=255)
   gatTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGAcc                                                             >  1:2441641/1‑69 (MQ=255)
     tttCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCgaga                                                         <  1:5473554/71‑1 (MQ=255)
        ccAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGAt                                                        >  1:4311285/1‑69 (MQ=255)
         cAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAg                                                     >  1:4729483/1‑71 (MQ=255)
          aGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGc                                                    <  1:4969249/71‑1 (MQ=255)
            ttttGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCg                                                   >  1:1980398/1‑70 (MQ=255)
               tGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCaa                                               >  1:1663509/1‑71 (MQ=255)
               tGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCa                                                <  1:3316384/70‑1 (MQ=255)
                gCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCaa                                               >  1:1603696/1‑70 (MQ=255)
                 cTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCaata                                             <  1:1693363/71‑1 (MQ=255)
                 cTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCaata                                             <  1:1877984/71‑1 (MQ=255)
                     tCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATc                                          <  1:1239412/70‑1 (MQ=255)
                     tCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATc                                          <  1:170256/70‑1 (MQ=255)
                     tCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATc                                          <  1:5217976/70‑1 (MQ=255)
                          ccgccgcGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGCCCGAGATTAGCGTCAATAATCagcagc                                    >  1:6771098/1‑71 (MQ=255)
                           cgccgcGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCg                                   <  1:93309/71‑1 (MQ=255)
                             ccgcGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGtt                                 >  1:2635851/1‑71 (MQ=255)
                               gcgcCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGttt                                >  1:1563770/1‑70 (MQ=255)
                                 gcCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCa                             <  1:3221659/71‑1 (MQ=255)
                                  cccAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCa                                                <  1:7123068/51‑1 (MQ=255)
                                     agcagGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTAt                         >  1:5985691/1‑71 (MQ=255)
                                     agcagGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTAt                         <  1:3186630/71‑1 (MQ=255)
                                     agcagGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTAt                         <  1:690113/71‑1 (MQ=255)
                                      gcagGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTAt                         >  1:5691386/1‑70 (MQ=255)
                                        agGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTc                              >  1:5042578/1‑63 (MQ=255)
                                         gggTCTTCGCCATATAGGCGCGTGCCTGACCgaga                                                         >  1:5427990/1‑35 (MQ=255)
                                         gggTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGtttt                               <  1:6771731/61‑1 (MQ=255)
                                         gggTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCct                     >  1:3560422/1‑71 (MQ=255)
                                         gggTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCct                     >  1:1630486/1‑71 (MQ=255)
                                          ggTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGttt                                <  1:7921804/59‑1 (MQ=255)
                                          ggTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCctc                    >  1:2599440/1‑71 (MQ=255)
                                            tCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGc                                                    >  1:1041301/1‑37 (MQ=255)
                                            tCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGc                                                    >  1:966661/1‑37 (MQ=255)
                                            tCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGc                                                    >  1:7182977/1‑37 (MQ=255)
                                            tCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGc                                                    >  1:1060389/1‑37 (MQ=255)
                                              ttCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCaataat                                           >  1:1338264/1‑44 (MQ=255)
                                              ttCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATca                                         >  1:161799/1‑46 (MQ=255)
                                              ttCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATc                                          <  1:110724/45‑1 (MQ=255)
                                              ttCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGttt                                <  1:2213014/55‑1 (MQ=255)
                                              ttCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGc                       >  1:1986659/1‑64 (MQ=255)
                                              ttCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGc                       >  1:2734940/1‑64 (MQ=255)
                                                cGCCATCTAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTAATTATGCCTCTCctgc              >  1:1021289/1‑71 (MQ=255)
                                                cGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTcc                 <  1:6795009/68‑1 (MQ=255)
                                                cGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCctgc              >  1:2939422/1‑71 (MQ=255)
                                                cGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCctg               <  1:2846680/70‑1 (MQ=255)
                                                 gCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTcc                 <  1:5485598/67‑1 (MQ=255)
                                                 gCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTcc                 <  1:7844130/67‑1 (MQ=255)
                                                   cATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCctg               <  1:2727164/67‑1 (MQ=255)
                                                    atatAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGt           >  1:7714674/1‑70 (MQ=255)
                                                    atatAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGt           >  1:3499175/1‑70 (MQ=255)
                                                    atatAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATCATGCCTCTCCTGCTGt           >  1:5404315/1‑70 (MQ=255)
                                                     tataGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTa                          <  1:6143712/54‑1 (MQ=255)
                                                     tataGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGt           >  1:6015544/1‑69 (MQ=255)
                                                     tataGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCa         <  1:4717328/71‑1 (MQ=255)
                                                     tataGGCGCATGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGttt                                <  1:3983888/48‑1 (MQ=255)
                                                     tataGGCGCATGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGttt                                <  1:3611015/48‑1 (MQ=255)
                                                     tataGGCGCATGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGttt                                <  1:3486722/48‑1 (MQ=255)
                                                      ataGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAg        >  1:5902590/1‑71 (MQ=255)
                                                       taGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCAt                            <  1:6109418/50‑1 (MQ=255)
                                                       taGGCGCATGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTc                              >  1:1419169/1‑48 (MQ=255)
                                                       taGGCGCATGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTc                              >  1:7774525/1‑48 (MQ=255)
                                                        aGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCAt                            <  1:7205686/49‑1 (MQ=255)
                                                        aGGCGCGCGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCAt                            <  1:5615971/49‑1 (MQ=255)
                                                         ggCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTa     <  1:3397454/71‑1 (MQ=255)
                                                         ggCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTa     <  1:1472563/71‑1 (MQ=255)
                                                         ggCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTa     <  1:3599954/71‑1 (MQ=255)
                                                          gcgcGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCctgctg            <  1:7753159/63‑1 (MQ=255)
                                                          gcgcGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTc          <  1:7177319/65‑1 (MQ=255)
                                                          gcgcGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTaa    <  1:432147/71‑1 (MQ=255)
                                                            gcgTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAg        >  1:4971842/1‑65 (MQ=255)
                                                            gcgTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAg        >  1:1768344/1‑65 (MQ=255)
                                                            gcgTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAg        >  1:6638726/1‑65 (MQ=255)
                                                             cgTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCAt                            <  1:5316446/44‑1 (MQ=255)
                                                              gTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTaaaa  >  1:7668538/1‑69 (MQ=255)
                                                              |                                                                    
GATGATTTCCAGTTTTGCTTTTCGCGCCGCCGCGCCCAGCAGGGTCTTCGCCATATAGGCGCGTGCCTGACCGAGATTAGCGTCAATAATCAGCAGCGTTTTCATTATGCCTCTCCTGCTGTCAGTTAAAA  >  AP009048/2264640‑2264770

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: