Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,806,471 G→A 13.6% P198P (CCG→CCA yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,806,4710GA13.6% 181.0 / 17.1 66P198P (CCG→CCAyicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (3/6);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GGCCAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTGG  >  AP009048/3806403‑3806538
                                                                    |                                                                   
ggCCAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGt                                                                    >  1:3409070/1‑70 (MQ=255)
  ccAGACGGTAGAGACCTGGAACCGGTACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCt                                                                 <  1:2820921/71‑1 (MQ=255)
  ccAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCt                                                                 <  1:1423655/71‑1 (MQ=255)
    aGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTa                                                                >  1:4000463/1‑70 (MQ=255)
    aGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTa                                                                >  1:2923996/1‑70 (MQ=255)
      aCGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCATTCTACAt                                                             <  1:917456/71‑1 (MQ=255)
       cGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACAt                                                             <  1:2569680/70‑1 (MQ=255)
          tAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACt                                                         <  1:1307357/71‑1 (MQ=255)
          tAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCATTCTACATGACt                                                         <  1:2949057/71‑1 (MQ=255)
              gaCCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAAcc                                                     >  1:2488032/1‑71 (MQ=255)
               aCCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAAc                                                      <  1:2057281/69‑1 (MQ=255)
               aCCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCg                                                    >  1:970264/1‑71 (MQ=255)
               aCCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCATTCTACATGACTAAcc                                                     <  1:4062959/70‑1 (MQ=255)
                 cTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTAc                                                               <  1:4291075/58‑1 (MQ=255)
                   ggAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGt                                                >  1:1864433/1‑71 (MQ=255)
                   ggAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGt                                                >  1:2094534/1‑71 (MQ=255)
                      aCCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTAt                                             <  1:3428702/71‑1 (MQ=255)
                       ccGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATg                                            >  1:361555/1‑71 (MQ=255)
                       ccGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATg                                            >  1:4252650/1‑71 (MQ=255)
                         gggACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACt                                                         >  1:279977/1‑56 (MQ=255)
                          ggACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGt                                                >  1:1366144/1‑64 (MQ=255)
                               gcggcACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACa                                                              >  1:1541747/1‑45 (MQ=255)
                                cggcACAAGTACTGAAGAGGCGTATAAAAATATCCCGTTCTACATGACTAAcc                                                     <  1:1671158/53‑1 (MQ=255)
                                cggcACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCt                                                                 <  1:2490273/41‑1 (MQ=255)
                                  gcACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTc                                 <  1:1519970/71‑1 (MQ=255)
                                  gcACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTc                                 <  1:1117383/71‑1 (MQ=255)
                                   cacaAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAAc                                                      <  1:722912/49‑1 (MQ=255)
                                   cacaAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCa                                <  1:4395032/71‑1 (MQ=255)
                                    acaAGTACTGAACAGTCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATg                                            <  1:4146706/58‑1 (MQ=255)
                                    acaAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTgg                                                 <  1:2639386/53‑1 (MQ=255)
                                    acaAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCaa                               >  1:366284/1‑71 (MQ=255)
                                       agcaCTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAatca                            >  1:3256648/4‑71 (MQ=255)
                                       aGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGtt                                               <  1:4502033/52‑1 (MQ=255)
                                        gTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTAt                                             <  1:1409099/53‑1 (MQ=255)
                                         tACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTa                                              <  1:2340589/51‑1 (MQ=255)
                                            tGAACAGGCGTATAAAAATATCCCGTTCTACATGACTa                                                        <  1:1965907/38‑1 (MQ=255)
                                             gAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATccc                        >  1:1232308/1‑69 (MQ=255)
                                              aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATgggg                                         >  1:2959666/1‑51 (MQ=255)
                                              aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATgggg                                         >  1:715237/1‑51 (MQ=255)
                                              aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                     >  1:1723469/1‑71 (MQ=255)
                                              aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                     <  1:1156700/71‑1 (MQ=255)
                                              aaCAGGCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                     >  1:3036984/1‑71 (MQ=255)
                                                cAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTAc                                      >  1:928055/1‑52 (MQ=255)
                                                cAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTAc                                      >  1:1331393/1‑52 (MQ=255)
                                                 aGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtgt                  >  1:1598920/1‑71 (MQ=255)
                                                 aGGCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtgt                  <  1:977145/71‑1 (MQ=255)
                                                    cGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGtc               <  1:817637/71‑1 (MQ=255)
                                                      tataAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGtctc             <  1:4207596/71‑1 (MQ=255)
                                                      tataAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGtctc             <  1:3308065/71‑1 (MQ=255)
                                                        taAAAATATCCCGTTCTACATGACTAACCGTGGTTa                                              <  1:3424266/36‑1 (MQ=255)
                                                        taAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTctct            <  1:4214598/70‑1 (MQ=255)
                                                        taAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCtt           >  1:3584495/1‑71 (MQ=255)
                                                        taAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCtt           >  1:3988732/1‑71 (MQ=255)
                                                         aaaaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATcccc                       <  1:4223166/58‑1 (MQ=255)
                                                          aaaaTATCCCGTTCTACATGACTAACCGTGGTTATgg                                           >  1:315136/1‑37 (MQ=255)
                                                          aaaaTATCCCGTTCTACATGACTAACCGTGGTTATg                                            >  1:1733543/1‑36 (MQ=255)
                                                          aaaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCtttt         <  1:3932069/71‑1 (MQ=255)
                                                          aaaaTATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCtttt         >  1:834914/1‑71 (MQ=255)
                                                           aaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAatcatc                          >  1:2910050/1‑53 (MQ=255)
                                                           aaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtgt                  <  1:1159404/61‑1 (MQ=255)
                                                           aaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCtttt         >  1:2079431/1‑70 (MQ=255)
                                                           aaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTg        <  1:1368113/71‑1 (MQ=255)
                                                             atatCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTg        >  1:3014332/1‑69 (MQ=255)
                                                                tCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGt    <  1:502574/70‑1 (MQ=255)
                                                                 cccATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTgg  >  1:69136/1‑71 (MQ=255)
                                                                   cGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCttt          <  1:2202604/61‑1 (MQ=255)
                                                                    |                                                                   
GGCCAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTGG  >  AP009048/3806403‑3806538

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: