Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,291 T→C 6.8% intergenic (+43/‑93) 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,2910TC 6.8% 165.9 / 2.4 59intergenic (+43/‑93)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base T (37/18);  new base C (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 = 9.40e-01

CCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACA  >  Exported/2434222‑2434359
                                                                     |                                                                    
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)
                              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)
                                                                     tGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAAcaca  >  1:4886431/1‑69 (MQ=255)
                                                                     |                                                                    
CCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACA  >  Exported/2434222‑2434359

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: