Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,053,941 G→A 9.6% intergenic (+99/+82) cca → / ← bacA fused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase/undecaprenyl pyrophosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,053,9410GA 9.6% 143.6 / 4.8 52intergenic (+99/+82)cca/bacAfused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase/undecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (28/19);  new base A (3/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

AACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCC  >  Exported/2053872‑2054007
                                                                      |                                                                  
aaCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGt                                                                   >  1:1842913/1‑71 (MQ=255)
 aCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGGAAGCGTCGCATCAGGCGtt                                                                  >  1:3875910/1‑71 (MQ=255)
 aCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGtt                                                                  <  1:129242/71‑1 (MQ=255)
  cGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGtt                                                                  >  1:3281464/1‑70 (MQ=255)
   gCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGa                                                                >  1:328791/1‑71 (MQ=255)
      ttATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGtt                                                                  <  1:1417792/66‑1 (MQ=255)
      ttATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTg                                                             >  1:5198770/1‑71 (MQ=255)
       tATCCGGCCTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCattgattg                                                             >  1:4402974/1‑71 (MQ=2)
        aTCCGGCCTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGc                                                            >  1:3497838/1‑71 (MQ=2)
        aTCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTg                                                             <  1:1559422/69‑1 (MQ=255)
            ggccTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGa                                                        >  1:5129818/4‑71 (MQ=14)
            ggCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGAt                                                       >  1:2698069/1‑71 (MQ=255)
                   aTGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGAtt                                                              <  1:3024428/58‑1 (MQ=255)
                   aTGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATgcgg                                                   <  1:1681334/69‑1 (MQ=255)
                    tGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATgcggc                                                  <  1:1625650/69‑1 (MQ=255)
                     gTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGt                                                >  1:4160627/1‑70 (MQ=255)
                       cAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGaa                                             >  1:4671056/1‑71 (MQ=255)
                        aaCGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                                            >  1:2545792/1‑71 (MQ=255)
                          cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGc                                          <  1:483099/71‑1 (MQ=255)
                          cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGc                                          >  1:688297/1‑71 (MQ=255)
                              gtttGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCtt                                       <  1:4293478/68‑1 (MQ=11)
                              gtttGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCtt                                       <  1:3890482/68‑1 (MQ=11)
                              aGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTa                                      <  1:1025895/71‑1 (MQ=255)
                               gTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCg                                                          <  1:4145213/50‑1 (MQ=255)
                               gTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTAt                                     <  1:5202251/71‑1 (MQ=255)
                                ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGCAt                                                       <  1:4536747/52‑1 (MQ=255)
                                 tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                                            >  1:2172949/1‑62 (MQ=255)
                                  gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGTAAACGCCTTATCCg                                  >  1:11775/1‑71 (MQ=14)
                                  gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGTAAACGCCTTATCCg                                  >  1:384209/1‑71 (MQ=14)
                                     ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGcc                               <  1:3706205/71‑1 (MQ=255)
                                        cTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTAc                            >  1:4226647/1‑71 (MQ=255)
                                         tGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGta                                               <  1:4217731/51‑2 (MQ=14)
                                          gATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACg                           >  1:5380770/1‑70 (MQ=255)
                                           aTAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGAGGCGTAAACGCCTTAt                                     >  1:243636/1‑59 (MQ=17)
                                           aTAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGc                          >  1:4815312/1‑70 (MQ=255)
                                             aaGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAAt                       >  1:1385256/1‑71 (MQ=255)
                                              aGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATc                      <  1:5636696/71‑1 (MQ=255)
                                               gACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCg                     >  1:5637974/1‑71 (MQ=255)
                                                aCGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATgcggc                                                  <  1:1660086/41‑1 (MQ=255)
                                                aCGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGt                    >  1:3671525/1‑71 (MQ=255)
                                                 cgcgCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGtt                   >  1:3079450/1‑71 (MQ=255)
                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGa                 >  1:946860/1‑70 (MQ=255)
                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGa                 >  1:1416237/1‑70 (MQ=255)
                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGAc                <  1:4357513/71‑1 (MQ=255)
                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGAc                <  1:5568931/71‑1 (MQ=255)
                                                     acaGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                            >  1:3716432/3‑58 (MQ=11)
                                                        gCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                                            >  1:1767830/1‑39 (MQ=255)
                                                        gCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCATATTCCGCCTACGCAATCGTTGACAACg            >  1:2247135/1‑71 (MQ=255)
                                                             gCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGcc          >  1:2436290/1‑68 (MQ=255)
                                                             gCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCaa        <  1:394944/70‑1 (MQ=255)
                                                             gCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAg       >  1:4752680/1‑71 (MQ=255)
                                                                  aGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTa                             >  1:1099773/1‑44 (MQ=255)
                                                                   ggCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATcc  <  1:1006173/70‑1 (MQ=255)
                                                                      |                                                                  
AACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCC  >  Exported/2053872‑2054007

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: