Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,357 A→G 13.3% intergenic (+109/‑27) 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,3570AG13.3% 196.8 / 17.4 75intergenic (+109/‑27)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base G (4/6);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

AAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTG  >  Exported/2434287‑2434417
                                                                      |                                                            
aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca                                                              >  1:2247963/1‑71 (MQ=255)
 aGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca                                                              >  1:4882327/1‑70 (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:5298159/71‑1 (MQ=255)
         agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                                                     <  1:2893665/71‑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)
               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:672959/1‑56 (MQ=255)
                  atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                                           >  1:5572557/1‑56 (MQ=255)
                  atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                                           >  1:3718940/1‑56 (MQ=255)
                    cccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                                            >  1:1279910/1‑69 (MQ=255)
                     ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                          >  1:3322455/1‑70 (MQ=255)
                     ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                          >  1:594413/1‑70 (MQ=255)
                     ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                          <  1:121317/70‑1 (MQ=255)
                     ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                          <  1:3958831/70‑1 (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:564370/71‑1 (MQ=255)
                         tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag                                     <  1:3058619/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)
                                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)
                                    tagtaAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCa                          >  1:2404347/1‑71 (MQ=255)
                                       taAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                                                          <  1:618693/36‑1 (MQ=255)
                                         aGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                                            >  1:1790166/1‑48 (MQ=255)
                                           tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc                   >  1:2311966/1‑71 (MQ=255)
                                           tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc                   >  1:2891341/1‑71 (MQ=255)
                                           tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc                   >  1:1550582/1‑71 (MQ=255)
                                           tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATCAAGAGTATGGGTTCCACAtcgtc                   >  1:5458551/1‑71 (MQ=255)
                                            gATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtcg                  >  1:4667734/1‑71 (MQ=255)
                                              ttCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                        <  1:1675654/47‑1 (MQ=255)
                                              ttCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                        <  1:3243562/47‑1 (MQ=255)
                                              ttCAGGGGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                        <  1:3009514/47‑1 (MQ=255)
                                                cAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTg               <  1:5350950/70‑1 (MQ=255)
                                                 aGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGt              <  1:2491598/70‑1 (MQ=255)
                                                    tGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                        >  1:3101460/1‑41 (MQ=255)
                                                    tGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGc          <  1:3294857/71‑1 (MQ=255)
                                                     gAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:1985509/1‑49 (MQ=255)
                                                     gAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:3942289/1‑49 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:962305/1‑48 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:5610772/1‑48 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:922984/1‑48 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:4383222/1‑48 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                               >  1:1439269/1‑48 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATgg                                >  1:4693506/1‑47 (MQ=255)
                                                      aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATgg                                >  1:3802152/1‑47 (MQ=255)
                                                         gAATGCAGCCAACACAGAGAGAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCa     >  1:4839376/1‑71 (MQ=255)
                                                         gAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCa     >  1:1944370/1‑71 (MQ=255)
                                                         gAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCa     >  1:1017520/1‑71 (MQ=255)
                                                          aaTGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCt                       >  1:2739377/1‑52 (MQ=255)
                                                          aaTGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCt                       >  1:2340350/1‑52 (MQ=255)
                                                          aaTGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCt                       >  1:2115110/1‑52 (MQ=255)
                                                            tGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTg  <  1:2741606/71‑1 (MQ=255)
                                                             gCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTg  <  1:4640623/70‑1 (MQ=255)
                                                                      |                                                            
AAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTG  >  Exported/2434287‑2434417

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: