Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,298 A→T 8.6% intergenic (+50/‑86) 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,2980AT 8.6% 159.2 / 6.3 58intergenic (+50/‑86)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (34/19);  new base T (2/3);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAG  >  Exported/2434233‑2434366
                                                                 |                                                                    
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)
             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:1442162/1‑50 (MQ=255)
                                           cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                           >  1:5087141/1‑50 (MQ=255)
                                             aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                 <  1:4275211/58‑1 (MQ=255)
                                             aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                 <  1:2491927/58‑1 (MQ=255)
                                             aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                 <  1:5226499/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)
                                                          cggCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCaga       <  1:3762691/70‑1 (MQ=2)
                                                            gCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga     <  1:769099/71‑1 (MQ=255)
                                                               agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg  <  1:5298159/71‑1 (MQ=255)
                                                               agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg  <  1:2893665/71‑1 (MQ=255)
                                                               agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg  <  1:4304180/71‑1 (MQ=255)
                                                                gagaATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa             >  1:279177/1‑59 (MQ=255)
                                                                 |                                                                    
TCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAG  >  Exported/2434233‑2434366

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: