Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,434,301 A→G 9.8% intergenic (+53/‑83) 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
*minE2,434,3010AG 9.8% 135.4 / 5.3 51intergenic (+53/‑83)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (30/16);  new base G (2/3);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTG  >  minE/2434233‑2434370
                                                                    |                                                                     
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: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:5102755/1‑35 (MQ=255)
                                  aaaaTCATTCAAACTGTATCAAGATGGCAAGAgaa                                                                       >  1:3690673/1‑35 (MQ=255)
                                  aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                                    <  1:4193246/70‑1 (MQ=255)
                                  aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                                    >  1:3893752/1‑70 (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: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:5226499/58‑1 (MQ=255)
                                             aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                     <  1:2491927/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:4304180/71‑1 (MQ=255)
                                                               agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg      <  1:2893665/71‑1 (MQ=255)
                                                               agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg      <  1:5298159/71‑1 (MQ=255)
                                                                gagaATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa                 >  1:279177/1‑59 (MQ=255)
                                                                  gaatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                   >  1:4161099/1‑55 (MQ=255)
                                                                   aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                        <  1:2961155/49‑1 (MQ=255)
                                                                   aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                        <  1:4740567/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)
                                                                    |                                                                     
TCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTG  >  minE/2434233‑2434370

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: