Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,176,396 C→T 10.0% M6I (ATG→ATA gatC ← galactitol‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,176,3960CT10.0% 149.7 / 6.3 50M6I (ATG→ATAgatCgalactitol‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (2/3);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGA  >  AP009048/2176327‑2176466
                                                                     |                                                                      
tAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCa                                                                       >  1:322469/1‑71 (MQ=255)
   tatTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCAt                                                                      <  1:3287615/69‑1 (MQ=255)
     ttttAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGAc                                                                   <  1:1817500/70‑1 (MQ=255)
        tAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCt                                                               >  1:2319318/1‑71 (MQ=255)
         aGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCt                                                               <  1:3304500/70‑1 (MQ=255)
          gAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCt                                                               <  1:862076/69‑1 (MQ=255)
           aaaaaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaa                                                            <  1:1389348/71‑1 (MQ=255)
            aaaaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaa                                                            >  1:2563151/1‑70 (MQ=255)
            aaaaaTAATGATGACAATCGGCACCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaa                                                            <  1:1129427/70‑1 (MQ=255)
             aaaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaaaa                                                          >  1:2993428/1‑71 (MQ=255)
               aataatGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACa                                                        <  1:3539069/71‑1 (MQ=255)
               aataatGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACa                                                        >  1:3524238/1‑71 (MQ=255)
                ataatGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACAt                                                       <  1:1372606/71‑1 (MQ=255)
                ataatGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGAc                                                                   >  1:3929064/1‑59 (MQ=255)
                 taatGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATa                                                      <  1:4429274/71‑1 (MQ=255)
                  aatGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAg                                                     >  1:3819795/1‑71 (MQ=255)
                   atgatgACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAgg                                                    >  1:2423222/1‑71 (MQ=255)
                   atgatgACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAg                                                     <  1:416897/70‑1 (MQ=255)
                           aaTCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTg                                                              <  1:2003615/53‑1 (MQ=255)
                             tCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAg                                                     >  1:3842335/1‑60 (MQ=255)
                              cGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACAt                                                       >  1:26255/1‑57 (MQ=255)
                              cGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACAt                                                       >  1:3419062/1‑57 (MQ=255)
                               ggCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTAt                                                                      <  1:1724402/41‑1 (MQ=255)
                                gcagcaTCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAg                                       <  1:2148858/71‑1 (MQ=255)
                                    caTCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACAt                                                       <  1:3663933/51‑1 (MQ=255)
                                    caTCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAg                                   <  1:3538573/71‑1 (MQ=255)
                                    caTCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAg                                   <  1:789394/71‑1 (MQ=255)
                                     aTCACCGTAGGGCCTAGGTCGAGAATATAACGTATGACTTCTGaaa                                                           <  1:1185530/46‑1 (MQ=255)
                                      tCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACt                                                                  <  1:4563643/38‑1 (MQ=255)
                                      tCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGt                                  >  1:2628087/1‑70 (MQ=255)
                                       cACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaaaa                                                          <  1:3708812/45‑1 (MQ=255)
                                       cACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGt                                  >  1:2472390/1‑69 (MQ=255)
                                        aCCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACAt                                                       >  1:1239529/1‑47 (MQ=255)
                                         ccGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGa                                      >  1:3664190/1‑63 (MQ=255)
                                              gggCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt                         <  1:2647448/71‑1 (MQ=255)
                                                  cGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAg                                       <  1:4520362/53‑1 (MQ=255)
                                                  cGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattttattt                     <  1:1291244/71‑1 (MQ=255)
                                                  cGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattttattt                     >  1:1621200/1‑71 (MQ=255)
                                                     ggTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTaa                                        <  1:1801756/49‑1 (MQ=255)
                                                      gTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGt                  >  1:3639063/1‑70 (MQ=255)
                                                        cGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAAt               >  1:3856105/1‑71 (MQ=255)
                                                         gagaATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATg              >  1:1016207/1‑71 (MQ=255)
                                                           gaATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCa                                >  1:3784886/1‑51 (MQ=255)
                                                           gaATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCt            >  1:3448087/1‑71 (MQ=255)
                                                           gaATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCt            >  1:2985448/1‑71 (MQ=255)
                                                           gaATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCt            >  1:949859/1‑71 (MQ=255)
                                                                taACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattttattt                     >  1:2074690/1‑57 (MQ=255)
                                                                taACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattttatt                      >  1:3487657/1‑56 (MQ=255)
                                                                  aCGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattt                          <  1:15054/50‑1 (MQ=255)
                                                                    gCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTcgataccg   <  1:996981/71‑1 (MQ=255)
                                                                     cATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTcgataccga  <  1:579307/71‑1 (MQ=255)
                                                                     |                                                                      
TAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGA  >  AP009048/2176327‑2176466

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: