Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,290 A→G 6.8% intergenic (+42/‑94) rhlB → / → gpp ATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,434,2900AG 6.8% 164.7 / 2.4 59intergenic (+42/‑94)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (37/18);  new base G (2/2);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACA  >  Exported/2434222‑2434357
                                                                    |                                                                   
ccGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATg                                                                   >  1:1212777/1‑71 (MQ=255)
ccGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTCTATCAAGATg                                                                   >  1:3154847/1‑71 (MQ=255)
       aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa                                                               >  1:2752797/1‑68 (MQ=255)
       aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa                                                               >  1:4905170/1‑68 (MQ=255)
           tcgtcgtTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatg                                                        <  1:2014098/71‑1 (MQ=255)
                 ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatga                                                       <  1:1491227/66‑1 (MQ=255)
                 ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatga                                                       <  1:1541547/66‑1 (MQ=255)
                 ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATccc                                                  >  1:2637809/1‑71 (MQ=255)
                  tCAGGTTAAATAAAAAATGTCATACTCAAAATAATTCAAACTGTATCAAGATGGCAAGAgaatga                                                       <  1:4206878/65‑1 (MQ=255)
                      gTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaat                                                         <  1:1229138/59‑1 (MQ=255)
                       ttAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaa                                                          >  1:4194020/1‑57 (MQ=255)
                        taaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGc                                           >  1:4197692/1‑71 (MQ=255)
                        taaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGc                                           >  1:3528107/1‑71 (MQ=255)
                        taaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCACGATGAGc                                           >  1:5042383/1‑71 (MQ=255)
                         aaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatgaat                                                     >  1:1404664/1‑60 (MQ=255)
                         aaataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCt                                          >  1:4493462/1‑71 (MQ=255)
                           ataaaAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTa                                        <  1:2318112/71‑1 (MQ=255)
                              aaaaaTGTCATACTCAAAATCATTCAAACTGTATCAAGa                                                                     >  1:175193/1‑39 (MQ=255)
                              aaaaaTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTg                                     >  1:2179355/1‑71 (MQ=255)
                                  aTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                  >  1:2023369/1‑70 (MQ=255)
                                     tCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtagt                                 <  1:1230803/68‑1 (MQ=255)
                                      taTACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCt                                          <  1:997666/57‑1 (MQ=9)
                                      taTACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGt                             >  1:4293028/2‑71 (MQ=9)
                                       aTACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaa                                                          <  1:32071/41‑1 (MQ=255)
                                       aTACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGt                             >  1:5601956/1‑70 (MQ=255)
                                        tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa                           >  1:195662/1‑71 (MQ=255)
                                        tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa                           >  1:3953088/1‑71 (MQ=255)
                                        tACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGa                           <  1:4024673/71‑1 (MQ=9)
                                          cTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAg                                            >  1:3876870/1‑52 (MQ=255)
                                           acAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTc                        >  1:2925141/2‑71 (MQ=12)
                                            cAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatgaa                                                      >  1:5197303/1‑40 (MQ=255)
                                             aaaaTCATTCAAACTGTATCAAGATGGCAAGAgaa                                                          >  1:5102755/1‑35 (MQ=255)
                                             aaaaTCATTCAAACTGTATCAAGATGGCAAGAgaa                                                          >  1:3690673/1‑35 (MQ=255)
                                             aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                       >  1:3893752/1‑70 (MQ=255)
                                             aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                       <  1:4193246/70‑1 (MQ=255)
                                             aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg                      >  1:5458564/1‑71 (MQ=255)
                                             aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg                      >  1:2517386/1‑71 (MQ=255)
                                              aaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCg                                                 <  1:3036918/43‑1 (MQ=255)
                                              aaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCt                                          <  1:1694394/50‑1 (MQ=255)
                                                aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGAt                                               >  1:4469841/1‑43 (MQ=255)
                                                aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt                         >  1:4116376/1‑65 (MQ=255)
                                                aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                    >  1:5154368/1‑70 (MQ=255)
                                                aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                    >  1:192398/1‑70 (MQ=255)
                                                   aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg                            >  1:2951237/1‑59 (MQ=255)
                                                   aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg                            >  1:5533870/1‑59 (MQ=255)
                                                   aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCACGATGAGCTTACTGTAGTAAGTg                            >  1:3781978/1‑59 (MQ=255)
                                                      cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                  >  1:5087141/1‑50 (MQ=255)
                                                      cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                  >  1:1442162/1‑50 (MQ=255)
                                                        aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                        <  1:4275211/58‑1 (MQ=255)
                                                        aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                        <  1:5226499/58‑1 (MQ=255)
                                                        aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                        <  1:2491927/58‑1 (MQ=255)
                                                             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt                         <  1:55517/52‑1 (MQ=255)
                                                             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg        <  1:1632874/69‑1 (MQ=255)
                                                             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg        <  1:3800996/69‑1 (MQ=255)
                                                             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGc       >  1:1956882/1‑70 (MQ=255)
                                                               tCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa    >  1:398427/1‑71 (MQ=255)
                                                                 aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca  >  1:2247963/1‑71 (MQ=255)
                                                                 aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAac   <  1:80674/70‑1 (MQ=255)
                                                                  aGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca  >  1:4882327/1‑70 (MQ=255)
                                                                   gATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGaa             >  1:3420785/1‑58 (MQ=255)
                                                                    |                                                                   
CCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACA  >  Exported/2434222‑2434357

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: