Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,279 A→G 5.5% intergenic (+31/‑105) 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,2790AG 5.5% 206.0 / 2.0 73intergenic (+31/‑105)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (41/28);  new base G (2/2);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TACTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGT  >  Exported/2434212‑2434339
                                                                   |                                                            
tACTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTGAAAATCATTCAAACTg                                                           >  1:687347/1‑71 (MQ=255)
tACTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTg                                                           >  1:3867163/1‑71 (MQ=255)
tACTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTg                                                           >  1:838508/1‑71 (MQ=255)
  cTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTa                                                         >  1:1126626/1‑71 (MQ=255)
  cTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTa                                                         >  1:3085494/1‑71 (MQ=255)
     gcgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATc                                                       >  1:4470732/1‑70 (MQ=255)
     gcgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATc                                                       <  1:4692085/70‑1 (MQ=255)
     gcgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCa                                                      <  1:2443418/71‑1 (MQ=255)
      cgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:5626912/71‑1 (MQ=255)
      cgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:3970646/71‑1 (MQ=255)
      cgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:331616/71‑1 (MQ=255)
          ccGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATg                                                 >  1:1212777/1‑71 (MQ=255)
          ccGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTCTATCAAGATg                                                 >  1:3154847/1‑71 (MQ=255)
              gTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATc                                                       >  1:2390883/1‑61 (MQ=255)
                 aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:1692950/60‑1 (MQ=255)
                 aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:4820784/60‑1 (MQ=255)
                 aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa                                             >  1:2752797/1‑68 (MQ=255)
                 aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa                                             >  1:4905170/1‑68 (MQ=255)
                  tcgtcgtcgTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:349124/59‑1 (MQ=255)
                  tcgtcgtcgTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:883947/59‑1 (MQ=255)
                    gtcgtcgtTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:3866174/57‑1 (MQ=255)
                    gtcgtcgtTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:4599190/57‑1 (MQ=255)
                     tcgtcgtTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatg                                      <  1:2014098/71‑1 (MQ=255)
                           ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatga                                     <  1:1541547/66‑1 (MQ=255)
                           ttCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAgaatga                                     <  1:1491227/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)
                                  taaataaaAAATGTCATACTCAAAATCATTCAAACTGTa                                                         >  1:1042971/1‑39 (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)
                                      taaaAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:3104090/39‑1 (MQ=255)
                                      taaaAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:823176/39‑1 (MQ=255)
                                      taaaAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                     <  1:4559240/39‑1 (MQ=255)
                                        aaaaaTGTCATACTCAAAATCATTCAAACTGTATCa                                                      <  1:3731587/36‑1 (MQ=255)
                                        aaaaaTGTCATACTCAAAATCATTCAAACTGTATCa                                                      <  1:2357833/36‑1 (MQ=255)
                                        aaaaaTGTCATACTCAAAATCATTCAAACTGTATCAAGa                                                   >  1:175193/1‑39 (MQ=255)
                                        aaaaaTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTg                   >  1:2179355/1‑71 (MQ=255)
                                         aaaaTGTCATACTCAAAATCATTCAAACTGTATCAAg                                                    <  1:175639/37‑1 (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:3953088/1‑71 (MQ=255)
                                                  tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa         >  1:195662/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:3690673/1‑35 (MQ=255)
                                                       aaaaTCATTCAAACTGTATCAAGATGGCAAGAgaa                                        >  1:5102755/1‑35 (MQ=255)
                                                       aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa     >  1:3893752/1‑70 (MQ=255)
                                                       aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa     <  1:4193246/70‑1 (MQ=255)
                                                       aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg    >  1:2517386/1‑71 (MQ=255)
                                                       aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg    >  1:5458564/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:5533870/1‑59 (MQ=255)
                                                             aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg          >  1:2951237/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:5226499/58‑1 (MQ=255)
                                                                  aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc      <  1:2491927/58‑1 (MQ=255)
                                                                  aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc      <  1:4275211/58‑1 (MQ=255)
                                                                   |                                                            
TACTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGT  >  Exported/2434212‑2434339

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: