Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,434,283 T→C 5.5% intergenic (+35/‑101) 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,2830TC 5.5% 207.7 / 2.0 73intergenic (+35/‑101)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGC  >  minE/2434214‑2434352
                                                                     |                                                                     
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:3970646/71‑1 (MQ=255)
    cgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                                  <  1:331616/71‑1 (MQ=255)
    cgcTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                                  <  1:5626912/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:4905170/1‑68 (MQ=255)
               aTCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCaa                                                          >  1:2752797/1‑68 (MQ=255)
                tcgtcgtcgTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                                  <  1:883947/59‑1 (MQ=255)
                tcgtcgtcgTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                                  <  1:349124/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: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)
                                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:823176/39‑1 (MQ=255)
                                    taaaAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                                  <  1:4559240/39‑1 (MQ=255)
                                    taaaAAATGTCATACTCAAAATCATTCAAACTGTATCaa                                                                  <  1:3104090/39‑1 (MQ=255)
                                      aaaaaTGTCATACTCAAAATCATTCAAACTGTATCa                                                                   <  1:2357833/36‑1 (MQ=255)
                                      aaaaaTGTCATACTCAAAATCATTCAAACTGTATCa                                                                   <  1:3731587/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: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: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:5226499/58‑1 (MQ=255)
                                                                aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                   <  1:4275211/58‑1 (MQ=255)
                                                                aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                   <  1:2491927/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)
                                                                     |                                                                     
CTGGCGCTCCGCGTAATCGTCGTCGTTCAGGTTAAATAAAAAATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGC  >  minE/2434214‑2434352

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: