Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,808,227 T→A 6.1% intergenic (+31/+22) yicI → / ← yicH predicted alpha‑glucosidase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,808,2270TA 6.1% 182.5 / 3.8 66intergenic (+31/+22)yicI/yicHpredicted alpha‑glucosidase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/18);  new base A (2/2);  total (46/20)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

AAGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGC  >  AP009048/3808161‑3808295
                                                                  |                                                                    
aaGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTc                                                                     >  1:2787541/1‑68 (MQ=255)
 aGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTAc                                                                 >  1:6236883/1‑71 (MQ=255)
  gCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTAcc                                                                >  1:6549137/1‑71 (MQ=255)
  gCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTAcc                                                                >  1:7555915/1‑71 (MQ=255)
  gCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGACTCAGACGTAcc                                                                <  1:423211/71‑1 (MQ=38)
  gCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGACTCAGACGTAcc                                                                <  1:5803287/71‑1 (MQ=38)
      cAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGc                                                            >  1:4997249/1‑71 (MQ=255)
      cAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGc                                                            >  1:6195289/1‑71 (MQ=255)
      cAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGACTCAGACGTACCCg                                                              >  1:7756496/1‑69 (MQ=38)
       aaGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGc                                                            >  1:6906619/1‑70 (MQ=255)
       aaGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGc                                                            >  1:6690556/1‑70 (MQ=255)
       aaGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGc                                                            >  1:2282972/1‑70 (MQ=255)
       aaGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGc                                                            >  1:5688892/1‑70 (MQ=255)
        aGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCga                                                          <  1:4471548/71‑1 (MQ=255)
        aGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCga                                                          <  1:1075503/71‑1 (MQ=255)
        aGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCga                                                          <  1:5915357/71‑1 (MQ=255)
          ggAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCgaga                                                        <  1:5081175/71‑1 (MQ=255)
           gAATGCGCTGACAATTACTTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCgaga                                                        <  1:870743/70‑1 (MQ=255)
                cgcTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATg                                                  >  1:4655483/1‑71 (MQ=255)
                 tcTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTa                                                                  <  1:3971572/53‑1 (MQ=255)
                  cTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGtt                                                <  1:4384156/71‑1 (MQ=255)
                  cTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGtt                                                >  1:2735854/1‑71 (MQ=255)
                    gACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTa                                              >  1:1648004/1‑71 (MQ=255)
                      cAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTAcc                                            <  1:5322553/71‑1 (MQ=255)
                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTc                                                                     >  1:5211141/1‑45 (MQ=255)
                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTc                                                                     >  1:1384909/1‑45 (MQ=255)
                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTc                                                                     >  1:2180168/1‑45 (MQ=255)
                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGt                                                                   <  1:1969392/47‑1 (MQ=255)
                       aaTTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGt                                                                   <  1:5859782/47‑1 (MQ=255)
                           aCGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCAcc                                       >  1:3899586/1‑71 (MQ=255)
                            cGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCgaga                                                        >  1:4620975/1‑53 (MQ=255)
                            cGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCg                                                           >  1:6084815/1‑50 (MQ=255)
                              ttGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATg                                                  <  1:7910088/57‑1 (MQ=255)
                                gTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTc                                  >  1:1526301/1‑71 (MQ=255)
                                gTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTc                                  >  1:3326855/1‑71 (MQ=255)
                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGt                                                 >  1:3785333/1‑53 (MQ=255)
                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCg                                <  1:6553923/70‑1 (MQ=255)
                                   aTCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGt                               >  1:163905/1‑71 (MQ=255)
                                      tATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGa                            >  1:5758721/1‑71 (MQ=255)
                                       aTTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCAc                                        >  1:1109696/1‑58 (MQ=255)
                                       aTTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCAc                                        >  1:1902647/1‑58 (MQ=255)
                                        ttGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGAt                           <  1:2234211/70‑1 (MQ=255)
                                            tGGAATGCCGGGTACGTCTGTGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAAc                      >  1:5443724/1‑71 (MQ=255)
                                            tGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAAc                      <  1:127729/71‑1 (MQ=255)
                                            tGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAAc                      >  1:2969868/1‑71 (MQ=255)
                                                aTGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCa                  >  1:4914250/1‑71 (MQ=255)
                                                 tGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTAc                                             <  1:3042175/43‑1 (MQ=255)
                                                 tGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTcc                                 >  1:6923110/1‑55 (MQ=255)
                                                 tGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCAc                 >  1:4945960/1‑71 (MQ=255)
                                                 tGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCAc                 >  1:2299821/1‑71 (MQ=255)
                                                   ccGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCg               >  1:6714431/1‑71 (MQ=255)
                                                     gggTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCa                                         >  1:1241755/1‑43 (MQ=255)
                                                       gTACGACTCAGACGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCt             >  1:4718005/1‑69 (MQ=38)
                                                         aCGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCg               <  1:2853359/65‑1 (MQ=255)
                                                         aCGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          >  1:2541920/1‑70 (MQ=255)
                                                         aCGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGa          >  1:2294425/1‑70 (MQ=255)
                                                          cGTCTGATTCGTACCCGGCGAGATTGATGTTTACCGc                                          <  1:4281589/37‑1 (MQ=255)
                                                           gTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAAc       >  1:2337390/1‑71 (MQ=255)
                                                            tCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCa                                         >  1:533082/1‑36 (MQ=255)
                                                            tCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAAc       <  1:3281051/70‑1 (MQ=255)
                                                             cTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGa     >  1:477164/1‑71 (MQ=255)
                                                              tGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGAt    >  1:4295483/1‑71 (MQ=255)
                                                               gAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCt             >  1:1956894/1‑61 (MQ=255)
                                                               gAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATg   >  1:6420432/1‑71 (MQ=255)
                                                               gAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATg   >  1:5053049/1‑71 (MQ=255)
                                                                aGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGc  >  1:4590554/1‑71 (MQ=255)
                                                                  tCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACg      <  1:6606756/65‑1 (MQ=255)
                                                                  tCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGc  >  1:1834425/1‑69 (MQ=255)
                                                                  |                                                                    
AAGCCTCAAGGGAATGCGCTGACAATTACGTTGTAATCTATTGATGGAATGCCGGGTACGTCTGAGTCGTACCCGGCGAGATTGATGTTTACCGCACCGGCTCCGTCGATGTAACTTCACCGCTGGAAACGATGC  >  AP009048/3808161‑3808295

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: