Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,322 G→C 15.8% intergenic (+74/‑62) 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,3220GC15.8% 212.6 / 26.2 82intergenic (+74/‑62)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base C (6/7);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

AAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTA  >  Exported/2434252‑2434384
                                                                      |                                                              
aaaaaTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTg                                                                >  1:2179355/1‑71 (MQ=255)
    aTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                                             >  1:2023369/1‑70 (MQ=255)
       tCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtagt                                                            <  1:1230803/68‑1 (MQ=255)
        taTACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGt                                                        >  1:4293028/2‑71 (MQ=9)
         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)
             acAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTc                                                   >  1:2925141/2‑71 (MQ=12)
               aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                                                  <  1:4193246/70‑1 (MQ=255)
               aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                                                  >  1:3893752/1‑70 (MQ=255)
               aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg                                                 >  1:5458564/1‑71 (MQ=255)
               aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg                                                 >  1:2517386/1‑71 (MQ=255)
                  aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt                                                    >  1:4116376/1‑65 (MQ=255)
                  aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                                               >  1:192398/1‑70 (MQ=255)
                  aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                                               >  1:5154368/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:2491927/58‑1 (MQ=255)
                          aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                                   <  1:5226499/58‑1 (MQ=255)
                          aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                                   <  1:4275211/58‑1 (MQ=255)
                               tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt                                                    <  1:55517/52‑1 (MQ=255)
                               tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg                                   <  1:3800996/69‑1 (MQ=255)
                               tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg                                   <  1:1632874/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)
                                               gaatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                 >  1:4161099/1‑55 (MQ=255)
                                                aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                                      <  1:4740567/49‑1 (MQ=255)
                                                aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                                      <  1:2961155/49‑1 (MQ=255)
                                                aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGAtt                 >  1:4042398/1‑70 (MQ=255)
                                                aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                <  1:5089868/71‑1 (MQ=255)
                                                 atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                    <  1:1837188/66‑1 (MQ=255)
                                                 atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                <  1:1370517/70‑1 (MQ=255)
                                                  tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGaa                                        >  1:1316831/1‑45 (MQ=255)
                                                  tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                 >  1:3701032/1‑52 (MQ=255)
                                                  tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                 >  1:4582472/1‑52 (MQ=255)
                                                  tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                 >  1:724454/1‑52 (MQ=255)
                                                  tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                <  1:955518/69‑1 (MQ=255)
                                                   gaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAacac                            <  1:4405530/56‑1 (MQ=255)
                                                    aatCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAg             >  1:2506493/1‑70 (MQ=9)
                                                    aatCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAg             >  1:1530604/1‑70 (MQ=9)
                                                     atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                          >  1:3718940/1‑56 (MQ=255)
                                                     atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                          >  1:672959/1‑56 (MQ=255)
                                                     atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                          >  1:5572557/1‑56 (MQ=255)
                                                       cccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa           >  1:1279910/1‑69 (MQ=255)
                                                        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg         <  1:3958831/70‑1 (MQ=255)
                                                        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg         <  1:121317/70‑1 (MQ=255)
                                                        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg         >  1:594413/1‑70 (MQ=255)
                                                        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg         >  1:3322455/1‑70 (MQ=255)
                                                        ccGATGAGCTTACTGTAGTAAGTGATTCAGGGGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg         <  1:2059265/70‑1 (MQ=255)
                                                        ccGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg         >  1:3565310/1‑70 (MQ=9)
                                                        ccGATGAGCTCACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg         <  1:4671174/70‑1 (MQ=255)
                                                         cGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa       <  1:1773347/71‑1 (MQ=255)
                                                           aTGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa       >  1:5325034/1‑69 (MQ=255)
                                                           aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg         <  1:1367101/67‑1 (MQ=9)
                                                           aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg         <  1:128231/67‑1 (MQ=9)
                                                           aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg         <  1:1081083/67‑1 (MQ=9)
                                                            tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag    <  1:3058619/71‑1 (MQ=255)
                                                            tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag    <  1:564370/71‑1 (MQ=255)
                                                            tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag    <  1:5655139/71‑1 (MQ=255)
                                                             gAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGCGACAGATTGAAGGATGAagag    >  1:2168033/1‑70 (MQ=255)
                                                             gAGCTTACTCAGGTAAGTGATTCGTGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATGAagtg    >  1:3677193/1‑68 (MQ=255)
                                                              aGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTa  >  1:3899367/1‑71 (MQ=255)
                                                              aGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATGAAgtgt   <  1:2636176/70‑1 (MQ=9)
                                                                cTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                         <  1:3382326/46‑1 (MQ=255)
                                                                cTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                         <  1:567898/46‑1 (MQ=255)
                                                                 ttACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                         <  1:1212478/45‑1 (MQ=255)
                                                                 ttACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                          <  1:5697365/44‑1 (MQ=255)
                                                                  tACTGTAGTAAGTGATTCAGGTGAAAGAATGAAGCCaa                               >  1:2602198/1‑38 (MQ=255)
                                                                   aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCa                                <  1:4020884/36‑1 (MQ=255)
                                                                   aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga                       >  1:3698710/1‑45 (MQ=255)
                                                                   aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga                       >  1:4518654/1‑45 (MQ=255)
                                                                   aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTa  >  1:5067942/1‑66 (MQ=255)
                                                                      gtagtaAGTGATTCACGTGAAAGAATGCAGCCAACACag                          >  1:3747104/1‑39 (MQ=255)
                                                                      |                                                              
AAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTA  >  Exported/2434252‑2434384

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: