Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,222,514 G→C 6.3% G409A (GGC→GCC)  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,5140GC 6.3% 390.4 / 10.5 126G409A (GGC→GCC) ebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (58/60);  new base C (4/4);  total (62/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 143 total aligned reads displayed.

TGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATC  >  AP009048/3222445‑3222584
                                                                     |                                                                      
  tgtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                                     <  1:2110493/71‑1 (MQ=255)
   gtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGc                                                                       <  1:3496776/68‑1 (MQ=255)
   gtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                                     >  1:2014113/1‑70 (MQ=255)
     gATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                                     >  1:4354423/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:154249/71‑1 (MQ=255)
       tatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTa                                                                <  1:6407565/71‑1 (MQ=255)
         tCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGc                                                              >  1:1919479/1‑71 (MQ=255)
           tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                                                                   >  1:6177278/1‑64 (MQ=255)
              ggCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCg                                                            <  1:2161734/68‑1 (MQ=255)
                ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATATCGCCGACAtt                                                                 <  1:5182558/61‑1 (MQ=38)
                ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATATCGCCGACAtt                                                                 <  1:5405375/61‑1 (MQ=38)
                ccTGTTTGTGATGGCGGAAACCGAAGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTa                                                       >  1:4211744/1‑71 (MQ=255)
                  tGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                                     >  1:135184/1‑71 (MQ=255)
                   ggttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGa                                                                     <  1:286493/52‑1 (MQ=255)
                   gTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGcc                                                             <  1:1693253/62‑1 (MQ=255)
                   gTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAc                                                      >  1:5787070/1‑69 (MQ=255)
                    tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGat                                                                    <  1:4186610/54‑1 (MQ=255)
                    tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           <  1:4814421/63‑1 (MQ=255)
                    tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           <  1:5840153/63‑1 (MQ=255)
                    tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                                     <  1:2382873/69‑1 (MQ=255)
                    tttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                                     <  1:2884511/69‑1 (MQ=255)
                     ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                                  >  1:3576333/1‑71 (MQ=255)
                     ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                                  >  1:5561100/1‑71 (MQ=255)
                      tgtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                                  >  1:2656534/1‑70 (MQ=255)
                      tgtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                                  >  1:6492658/1‑70 (MQ=255)
                       gtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgacg                                                 <  1:4163481/70‑1 (MQ=255)
                       gtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcga                                                   >  1:3785725/1‑68 (MQ=255)
                         gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                                        <  1:6101283/45‑1 (MQ=255)
                         gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                                        <  1:1239818/45‑1 (MQ=255)
                          aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc                                             >  1:7249234/1‑71 (MQ=255)
                          aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc                                             >  1:1119823/1‑71 (MQ=255)
                          aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc                                             >  1:4637556/1‑71 (MQ=255)
                          aTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc                                              <  1:4751617/70‑1 (MQ=255)
                           tGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATAtt                                                                 <  1:7842407/50‑1 (MQ=255)
                            ggCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCg                                            >  1:900684/1‑70 (MQ=255)
                             gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           <  1:4804947/54‑1 (MQ=255)
                             gCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac                                                  >  1:697503/1‑63 (MQ=255)
                              cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAg                                         <  1:6373049/71‑1 (MQ=255)
                              cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAg                                         <  1:5132735/71‑1 (MQ=255)
                              cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAg                                         <  1:7692023/71‑1 (MQ=255)
                               ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc                                                     >  1:2537456/1‑58 (MQ=255)
                               ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                                        >  1:6420729/1‑71 (MQ=255)
                               ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                                        >  1:6352803/1‑71 (MQ=255)
                               ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                                        >  1:919776/1‑71 (MQ=255)
                                 aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                                        >  1:2402557/1‑37 (MQ=255)
                                 aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                                        >  1:7216908/1‑37 (MQ=255)
                                 aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                                                        >  1:6821925/1‑37 (MQ=255)
                                 aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                                                                   >  1:3658354/1‑42 (MQ=255)
                                 aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                                                                   >  1:307642/1‑42 (MQ=255)
                                 aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcga                                                   >  1:6674146/1‑58 (MQ=255)
                                     cGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                                  <  1:6428472/71‑1 (MQ=255)
                                     cGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                                  <  1:5888691/71‑1 (MQ=255)
                                       aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                                  >  1:1267208/1‑69 (MQ=255)
                                       aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                                  >  1:5875574/1‑69 (MQ=255)
                                       aCGTCGAATCGCACGGCTTTGCTAATATCGCCGacat                                                                  >  1:6535172/1‑37 (MQ=21)
                                       aCGTCGAATCGCACGGCTTTGCTAATATCGCCGacat                                                                  >  1:4893642/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)
                                          tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCg                                            <  1:861106/56‑1 (MQ=255)
                                          tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTc                             <  1:7580434/71‑1 (MQ=255)
                                            gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           >  1:536200/1‑39 (MQ=255)
                                            gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           >  1:4023159/1‑39 (MQ=255)
                                            gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           >  1:5113229/1‑39 (MQ=255)
                                            gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                                                           >  1:7204700/1‑39 (MQ=255)
                                            gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAg                               <  1:2305156/67‑1 (MQ=255)
                                            gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAg                               <  1:4487293/67‑1 (MQ=255)
                                             aaTCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAg                               <  1:2295913/66‑1 (MQ=255)
                                               tCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGt                        >  1:7772111/1‑71 (MQ=255)
                                                 gCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGt                                        >  1:4269112/1‑53 (MQ=255)
                                                      gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgcgc                 <  1:2559231/71‑1 (MQ=255)
                                                         ttGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCAtt              >  1:2485572/1‑71 (MQ=255)
                                                          tGCTAATGTCGGCGATATTAGCCGTATTACcgacga                                                <  1:7473972/36‑1 (MQ=255)
                                                            cTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                                  <  1:1208167/48‑1 (MQ=255)
                                                            cTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGaaa                                  <  1:3415549/48‑1 (MQ=255)
                                                            cTAATATCGCCGACATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGtt           >  1:4260188/1‑71 (MQ=39)
                                                            cTAATATCGCCGACATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGt            >  1:3546129/1‑70 (MQ=39)
                                                               aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCa                                          <  1:7468206/37‑1 (MQ=255)
                                                               aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCt                            >  1:6229095/1‑51 (MQ=255)
                                                               aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACg                         >  1:2775568/1‑54 (MQ=255)
                                                               aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACg                         >  1:3144124/1‑54 (MQ=255)
                                                                tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGa                     <  1:7601017/57‑1 (MQ=255)
                                                                  tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCg                      <  1:5103454/54‑1 (MQ=255)
                                                                  tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgc                   >  1:2259548/1‑57 (MQ=255)
                                                                  tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgc                   >  1:423104/1‑57 (MQ=255)
                                                                  tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGcc       >  1:1031881/1‑69 (MQ=255)
                                                                   cGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgcg                  <  1:2259330/57‑1 (MQ=255)
                                                                   cGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgcg                  <  1:1892721/57‑1 (MQ=255)
                                                                   cGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAgcg                  <  1:1082703/57‑1 (MQ=255)
                                                                    ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTAc                          >  1:3558813/1‑48 (MQ=255)
                                                                    ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTAc                          >  1:6158217/1‑48 (MQ=255)
                                                                    ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTg             <  1:3210260/61‑1 (MQ=255)
                                                                     gCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTa                           >  1:3266410/1‑46 (MQ=255)
                                                                     gCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTa                           >  1:806448/1‑46 (MQ=255)
                                                                     gCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTa                           >  1:3452776/1‑46 (MQ=255)
                                                                     gCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTa                           >  1:6249435/1‑46 (MQ=255)
                                                                     gCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCa                >  1:5493106/1‑57 (MQ=255)
                                                                     |                                                                      
TGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATC  >  AP009048/3222445‑3222584

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: