Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,222,510 G→A 6.0% V408I (GTC→ATC)  ebgA → cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,222,5100GA 6.0% 432.6 / 11.1 133V408I (GTC→ATC) ebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (59/66);  new base A (4/4);  total (63/70)
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.40e-01

Only 100 of 133 total aligned reads displayed.

AACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTT  >  AP009048/3222442‑3222575
                                                                    |                                                                 
aaCTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTc                                                                 <  1:469384/71‑1 (MQ=255)
 aCTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                >  1:2968538/1‑71 (MQ=255)
 aCTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:2148082/71‑1 (MQ=255)
 aCTGTGTGATATCTACGGCCTATTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:7507135/71‑1 (MQ=255)
   tgtgtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGc                                                              >  1:2399760/1‑71 (MQ=255)
     tgtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGc                                                              <  1:1865307/69‑1 (MQ=255)
     tgtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                            <  1:2110493/71‑1 (MQ=255)
      gtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGc                                                              <  1:3496776/68‑1 (MQ=255)
        gATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                            >  1:4354423/1‑68 (MQ=255)
        gATATCTACGACCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                            >  1:6284936/1‑68 (MQ=255)
         atatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGatat                                                         >  1:1720339/1‑70 (MQ=255)
         atatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATAtt                                                        <  1:2884763/71‑1 (MQ=255)
          tatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGatat                                                         >  1:3923442/1‑69 (MQ=255)
          tatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTa                                                       <  1:6407565/71‑1 (MQ=255)
          tatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTa                                                       <  1:154249/71‑1 (MQ=255)
              tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                                                          >  1:6177278/1‑64 (MQ=255)
              tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCg                                                   >  1:4811664/1‑71 (MQ=255)
                cGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGcc                                                    <  1:6391694/68‑1 (MQ=255)
                cGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCg                                                   <  1:2121619/69‑1 (MQ=255)
                  gCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTAtt                                               <  1:5394092/71‑1 (MQ=255)
                   ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATATCGCCGACAtt                                                        <  1:5182558/61‑1 (MQ=38)
                   ccTGTTTGTGATGGCGGAAACCGAAGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTa                                              >  1:4211744/1‑71 (MQ=255)
                     tGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                            >  1:135184/1‑71 (MQ=255)
                      gTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGcc                                                    <  1:1693253/62‑1 (MQ=255)
                      gTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGc                                                     <  1:6437200/61‑1 (MQ=255)
                      gTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAc                                             >  1:5787070/1‑69 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:7161311/49‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:5997816/49‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:1954035/49‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:1941551/49‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                                                <  1:7918059/49‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGat                                                           <  1:4186610/54‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                            <  1:2884511/69‑1 (MQ=255)
                       tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                            <  1:2382873/69‑1 (MQ=255)
                        ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                         >  1:5561100/1‑71 (MQ=255)
                        ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                         >  1:3576333/1‑71 (MQ=255)
                         tgtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                         >  1:6492658/1‑70 (MQ=255)
                         tgtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                         >  1:2656534/1‑70 (MQ=255)
                          gtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgacga                                       >  1:3179019/1‑71 (MQ=255)
                          gtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgacg                                        <  1:4163481/70‑1 (MQ=255)
                          gtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcga                                          >  1:3785725/1‑68 (MQ=255)
                           tgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgacga                                       <  1:5858276/70‑1 (MQ=255)
                            gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                               <  1:1239818/45‑1 (MQ=255)
                            gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                               <  1:6101283/45‑1 (MQ=255)
                            gATGGCGGAAACCGACGGCGAATCGCACGGCTTTGCTAATGTCgg                                                               <  1:7895175/45‑1 (MQ=255)
                             aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc                                    >  1:1119823/1‑71 (MQ=255)
                             aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc                                    >  1:7249234/1‑71 (MQ=255)
                             aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc                                    >  1:4637556/1‑71 (MQ=255)
                             aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc                                     <  1:3559074/70‑1 (MQ=255)
                             aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc                                     <  1:4751617/70‑1 (MQ=255)
                              tGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATAtt                                                        <  1:7842407/50‑1 (MQ=255)
                              tGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCg                                   >  1:6948956/1‑71 (MQ=255)
                               ggCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCg                                   >  1:900684/1‑70 (MQ=255)
                                gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                  <  1:4804947/54‑1 (MQ=255)
                                gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGc                                  <  1:228498/70‑1 (MQ=255)
                                gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTGGCGATATTAGCCGTATTACCGACGATCCGCa                                 >  1:5630731/1‑70 (MQ=255)
                                 cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCa                                 >  1:862584/1‑70 (MQ=255)
                                 cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAg                                <  1:7692023/71‑1 (MQ=255)
                                 cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAg                                <  1:5132735/71‑1 (MQ=255)
                                  ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                               >  1:6352803/1‑71 (MQ=255)
                                  ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                               >  1:6420729/1‑71 (MQ=255)
                                  ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                               >  1:919776/1‑71 (MQ=255)
                                   gAAACCGACGTCGAATCGCACGGCTTTGCTCATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTg                              <  1:7741713/71‑1 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                               >  1:6821925/1‑37 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                               >  1:7216908/1‑37 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                                                          >  1:3658354/1‑42 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                                                          >  1:307642/1‑42 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTAtt                                               >  1:3521597/1‑53 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcga                                          >  1:6674146/1‑58 (MQ=255)
                                        cGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                         <  1:5888691/71‑1 (MQ=255)
                                          aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                         >  1:5875574/1‑69 (MQ=255)
                                          aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                         >  1:1267208/1‑69 (MQ=255)
                                          aCGTCGAATCGCACGGCTTTGCTAATATCGCCGacat                                                         >  1:4893642/1‑37 (MQ=21)
                                          aCGTCGAATCGCACGGCTTTGCTAATATCGCCGacat                                                         >  1:6535172/1‑37 (MQ=21)
                                            gTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                            >  1:2354259/1‑48 (MQ=255)
                                            gTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGt                     <  1:7048147/71‑1 (MQ=255)
                                            gTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGt                     >  1:2414969/1‑71 (MQ=255)
                                             tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGc                                  <  1:2508183/57‑1 (MQ=255)
                                             tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTc                    <  1:1931433/71‑1 (MQ=255)
                                             tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTc                    <  1:1175409/71‑1 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                  >  1:7204700/1‑39 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                  >  1:5113229/1‑39 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                  >  1:536200/1‑39 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                  >  1:4023159/1‑39 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAg                      <  1:4487293/67‑1 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAg                      <  1:2305156/67‑1 (MQ=255)
                                                aaTCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAg                      <  1:2295913/66‑1 (MQ=255)
                                                    gCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                               >  1:4269112/1‑53 (MQ=255)
                                                         gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgcgc        <  1:2559231/71‑1 (MQ=255)
                                                          cTTTGCTAATGTCGGCGATATTAGACGTATTACCGACGATCCGCAGTGGGAAAAAGTCTAAGTCGAgcg         >  1:7467087/1‑69 (MQ=255)
                                                            ttGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCAtt     >  1:2485572/1‑71 (MQ=255)
                                                             tGCTAATGTCGGCGATATTAGCCGTATTACcgacga                                       <  1:7473972/36‑1 (MQ=255)
                                                               cTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                         <  1:3415549/48‑1 (MQ=255)
                                                               cTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                         <  1:1208167/48‑1 (MQ=255)
                                                               cTAATATCGCCGACATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGtt  >  1:4260188/1‑71 (MQ=39)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCt                   >  1:6229095/1‑51 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACg                >  1:2775568/1‑54 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACg                >  1:3144124/1‑54 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGt   <  1:753140/67‑1 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGa            <  1:7601017/57‑1 (MQ=255)
                                                                    |                                                                 
AACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTT  >  AP009048/3222442‑3222575

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: