Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,053,962 G→A 13.0% intergenic (+120/+61) 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,9620GA13.0% 121.2 / 7.7 46intergenic (+120/+61)cca/bacAfused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase/undecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (27/13);  new base A (4/2);  total (31/15)
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.59e-01

CAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAATTAAAAGA  >  Exported/2053894‑2054030
                                                                    |                                                                    
cAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGaa                                                                    >  1:4671056/1‑71 (MQ=255)
 aaCGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                                                                   >  1:2545792/1‑71 (MQ=255)
   cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGc                                                                 >  1:688297/1‑71 (MQ=255)
   cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGc                                                                 <  1:483099/71‑1 (MQ=255)
       gtttGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCtt                                                              <  1:4293478/68‑1 (MQ=11)
       gtttGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCtt                                                              <  1:3890482/68‑1 (MQ=11)
       aGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTa                                                             <  1:1025895/71‑1 (MQ=255)
        gTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTAt                                                            <  1:5202251/71‑1 (MQ=255)
          tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                                                                   >  1:2172949/1‑62 (MQ=255)
           gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGTAAACGCCTTATCCg                                                         >  1:384209/1‑71 (MQ=14)
           gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGTAAACGCCTTATCCg                                                         >  1:11775/1‑71 (MQ=14)
              ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGcc                                                      <  1:3706205/71‑1 (MQ=255)
                 cTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTAc                                                   >  1:4226647/1‑71 (MQ=255)
                   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)
                         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)
                                                  gATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTa                  >  1:2078098/1‑71 (MQ=255)
                                                        cGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCAt                           >  1:4484728/1‑56 (MQ=255)
                                                        cGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCAt                           >  1:200924/1‑56 (MQ=255)
                                                         ggATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGcc                                 <  1:599692/49‑1 (MQ=255)
                                                         ggATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCt             >  1:4384937/1‑69 (MQ=255)
                                                         ggATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTc            >  1:3151562/1‑70 (MQ=255)
                                                            tGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAAt         >  1:742562/1‑70 (MQ=255)
                                                            tGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAAt         >  1:3317704/1‑70 (MQ=255)
                                                             gcggcgTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCaa          >  1:1907140/1‑68 (MQ=255)
                                                              cggcgTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAATTa       >  1:355776/1‑70 (MQ=255)
                                                                  gTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAATTAAAAGa  <  1:5182873/71‑1 (MQ=255)
                                                                   tGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATcc                         <  1:2846034/47‑1 (MQ=255)
                                                                   tGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAATTAAAAGa  <  1:4857100/70‑1 (MQ=255)
                                                                    |                                                                    
CAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGACAACGCCAAGCATCCGACACTATTCCTCAATTAAAAGA  >  Exported/2053894‑2054030

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: