Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,176,543 C→G 5.9% L53F (TTG→TTC gatB ← galactitol‑specific enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,176,5430CG 5.9% 417.4 / 9.5 133L53F (TTG→TTCgatBgalactitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base C (71/54);  new base G (2/6);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

AATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTA  >  AP009048/2176474‑2176612
                                                                     |                                                                     
aaTAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATc                                                                       <  1:6692379/70‑1 (MQ=255)
aaTAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCa                                                                      <  1:6426851/71‑1 (MQ=255)
  tAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    <  1:5184509/71‑1 (MQ=255)
      aGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATGAAATGca                                                                >  1:3156233/1‑71 (MQ=255)
       ggCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGcac                                                               <  1:4808873/71‑1 (MQ=255)
           tGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:7724404/1‑62 (MQ=255)
           tGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:5355323/1‑62 (MQ=255)
           tGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:4878905/1‑62 (MQ=255)
           tGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:799108/1‑62 (MQ=255)
           tGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAcca                                                           <  1:3760843/71‑1 (MQ=255)
           tGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAcca                                                           <  1:2482176/71‑1 (MQ=255)
            gCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:267950/1‑61 (MQ=255)
            gCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:3934933/1‑61 (MQ=255)
            accGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:595701/2‑61 (MQ=255)
             ccGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATTTGATATGcacca                                                             <  1:3085592/67‑3 (MQ=38)
             ccGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGcacca                                                             <  1:3622762/67‑3 (MQ=255)
             ccGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGcacca                                                             <  1:2650792/67‑3 (MQ=255)
                tGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGTTGCATATCAAATGCACAccatcca                                                       >  1:7467022/1‑70 (MQ=255)
                 gAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCa                                                                      <  1:2147813/54‑1 (MQ=255)
                 gAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccat                                                          >  1:4666413/1‑66 (MQ=255)
                 gAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccat                                                          >  1:6913727/1‑66 (MQ=255)
                  aaCTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGcaca                                                              <  1:3858049/61‑1 (MQ=255)
                      aaCGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    <  1:2364073/51‑1 (MQ=255)
                      aaCGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccat                                                          <  1:2181395/61‑1 (MQ=255)
                       aCGGAATATCGCCAAAACTACGTTCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGtt                                               >  1:2391208/1‑71 (MQ=255)
                       aCGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGtt                                               >  1:3803061/1‑71 (MQ=255)
                       aCGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGtt                                               >  1:2998293/1‑71 (MQ=255)
                         ggAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTc                                             <  1:1782602/71‑1 (MQ=255)
                           aaTATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATc                                                                       >  1:4018404/1‑43 (MQ=255)
                           aaTATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATGAAATGCACACCATCCATATAGGTTTCTa                                           >  1:6068937/1‑71 (MQ=255)
                           aaTATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatc                                                         >  1:1596878/1‑57 (MQ=255)
                           aaTATCCCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTa                                           >  1:3978454/1‑71 (MQ=255)
                             tatCGCCAAAACTACGATCCACTTTGGCAGTTGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         <  1:6779555/71‑1 (MQ=255)
                             tatCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAgg                                                 <  1:2138118/63‑1 (MQ=255)
                             tatCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         >  1:1403404/1‑71 (MQ=255)
                              atCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         <  1:1396194/70‑1 (MQ=255)
                               tCGCCAAAGCTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:7219705/1‑42 (MQ=255)
                               tCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaaa                                                                    >  1:483277/1‑42 (MQ=255)
                               tCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAttt                                        >  1:1351470/1‑70 (MQ=255)
                                cGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGcaca                                                              >  1:1237747/1‑47 (MQ=255)
                                cGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGcaca                                                              >  1:4885847/1‑47 (MQ=255)
                                cGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatcc                                                        <  1:728644/53‑1 (MQ=255)
                                   cAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCGata                                                     <  1:7191104/53‑1 (MQ=255)
                                   cAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCata                                                     <  1:3880719/53‑1 (MQ=255)
                                    aaaaCTTCGATCCACTTTGGCAGTGGTGCATATCa                                                                      <  1:7738170/35‑1 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTTCATATCAAATGCACACCATCCata                                                     <  1:3511865/52‑1 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatcc                                                        >  1:1848255/1‑49 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatcc                                                        >  1:580423/1‑49 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTc                                       >  1:6787230/1‑66 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTc                                       >  1:7575842/1‑66 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCAtt                                    <  1:5717274/69‑1 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCAtt                                    <  1:67450/69‑1 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCAtt                                    <  1:7373973/69‑1 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTaa                                  <  1:4645506/71‑1 (MQ=255)
                                    aaaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAACATCCATATAGGTTTCTAt                                          >  1:4558597/1‑63 (MQ=255)
                                     aaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         >  1:3479763/1‑63 (MQ=255)
                                     aaaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         >  1:5189052/1‑63 (MQ=255)
                                      aaCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAcc                                >  1:2748553/1‑71 (MQ=255)
                                       aCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAccc                               <  1:4869040/71‑1 (MQ=255)
                                       aCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAccc                               <  1:6885756/71‑1 (MQ=255)
                                       aCTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAccc                               <  1:5613456/71‑1 (MQ=255)
                                        cTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatcca                                                       >  1:1118588/1‑46 (MQ=255)
                                        cTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatcca                                                       >  1:1199920/1‑46 (MQ=255)
                                        cTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAccatcc                                                        >  1:3803036/1‑45 (MQ=255)
                                           cGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGaca                           >  1:6882315/1‑71 (MQ=255)
                                           cGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGaca                           >  1:7430800/1‑71 (MQ=255)
                                           cGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGa                             <  1:6034330/69‑1 (MQ=255)
                                            gATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGacac                          <  1:5936408/71‑1 (MQ=255)
                                              tCCACTATGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTg                        >  1:8088/1‑71 (MQ=255)
                                                cACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTgg                       <  1:2020775/70‑1 (MQ=255)
                                                   tttGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGAtt                    >  1:4971179/1‑70 (MQ=255)
                                                    ttGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCAGACACTGGATTaa                  >  1:2434684/1‑71 (MQ=255)
                                                     tGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTaa                  <  1:3247786/70‑1 (MQ=255)
                                                       gCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAAtt                >  1:2222377/1‑70 (MQ=255)
                                                        cAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCa              <  1:5412531/71‑1 (MQ=255)
                                                         aGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTaa                                  <  1:6241930/50‑1 (MQ=255)
                                                         aGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCaa             <  1:446637/71‑1 (MQ=255)
                                                          gtggtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCaa             >  1:5101605/1‑70 (MQ=255)
                                                           tggtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTc                                       >  1:2880238/1‑43 (MQ=255)
                                                           tggtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTc                                       >  1:6532707/1‑43 (MQ=255)
                                                           tggtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTc                                       >  1:2010849/1‑43 (MQ=255)
                                                             gtgCATTTGATATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGAtt                    <  1:2094787/60‑1 (MQ=255)
                                                             gtgCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         <  1:1589934/39‑1 (MQ=255)
                                                             gtgCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                         <  1:5650010/39‑1 (MQ=255)
                                                             gtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAc                                 >  1:4452723/1‑47 (MQ=255)
                                                             gtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAc                                 >  1:2274408/1‑47 (MQ=255)
                                                              tgCATTTGATATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGa                      <  1:1721400/57‑1 (MQ=255)
                                                               gCATATCAAATGCACACCATCCATATAGGTTTCTa                                           >  1:5390745/1‑35 (MQ=255)
                                                                cATATCAAATGCACACCATCCATATAGGTTTCTAttt                                        <  1:3440135/37‑1 (MQ=255)
                                                                 atatCAAATGCACACCATCCATATAGGTTTCTATTTCATTaa                                  <  1:4974375/42‑1 (MQ=255)
                                                                  tatCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGa                             <  1:7348069/46‑1 (MQ=255)
                                                                  tatCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAatat    >  1:1325563/1‑71 (MQ=255)
                                                                  tatCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAata     >  1:2562370/1‑70 (MQ=255)
                                                                    tCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGa                             >  1:5626268/1‑44 (MQ=255)
                                                                    tCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGa                             >  1:6744966/1‑44 (MQ=255)
                                                                    tCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAAtatta  >  1:2664997/1‑71 (MQ=255)
                                                                    tCAAATGCACACCATCCATATAGGTTTATATTTCATTAACCCGa                             >  1:5252526/1‑44 (MQ=255)
                                                                     cAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCg                              >  1:5308920/1‑42 (MQ=255)
                                                                     cAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCg                              >  1:5944918/1‑42 (MQ=255)
                                                                     cAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAAc            <  1:1160869/60‑1 (MQ=255)
                                                                     |                                                                     
AATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTA  >  AP009048/2176474‑2176612

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