Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 842,247 A→G 27.9% P109P (CCA→CCG rssB → response regulator of RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE842,2470AG27.9% 76.2 / 31.0 75P109P (CCA→CCGrssBresponse regulator of RpoS
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base G (2/19);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 6.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAA  >  minE/842185‑842309
                                                              |                                                              
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAg                                                               >  1:2603445/1‑64 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:2764632/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:2681343/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:1300889/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:2434588/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:1373085/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:3280352/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:1943196/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:345369/1‑65 (MQ=255)
gAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGt                                                              >  1:582683/1‑65 (MQ=255)
   aaTATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:981020/67‑1 (MQ=255)
   aaTATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:3056463/67‑1 (MQ=255)
   aaTATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:2605483/67‑1 (MQ=255)
   aaTATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:2132268/67‑1 (MQ=255)
   aaTATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:349066/67‑1 (MQ=255)
   aaTATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:1582528/67‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:538269/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:301774/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:2810511/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:1094162/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:2734064/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:1011407/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:973122/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:3095635/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:2576148/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:3118346/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:1329933/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:2043232/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:3387651/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:3427965/69‑1 (MQ=255)
     tatGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:1759341/69‑1 (MQ=255)
      atGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:3280661/68‑1 (MQ=255)
      atGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:1658326/68‑1 (MQ=255)
      atGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:1576959/68‑1 (MQ=255)
      atGGCAGATATTGCCAAAGCGTTACGTCTGGGTGTTGAAGATGTTTTGCTGAAACAGGTTAAAGATCt                                                     <  1:512019/68‑1 (MQ=255)
              tatTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCtt                                            >  1:1390638/1‑69 (MQ=255)
              tatTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCtt                                            >  1:1451356/1‑69 (MQ=255)
                   ccAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAg                                                         <  1:2859894/51‑1 (MQ=255)
                     aaaGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCga                                      <  1:2185119/68‑1 (MQ=255)
                     aaaGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCga                                      <  1:898052/68‑1 (MQ=255)
                                tCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTg                           <  1:2955552/68‑1 (MQ=255)
                                tCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTg                           <  1:2690678/68‑1 (MQ=255)
                                tCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTg                           <  1:690604/68‑1 (MQ=255)
                                tCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTg                           <  1:710652/68‑1 (MQ=255)
                                tCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTg                           <  1:3217295/68‑1 (MQ=255)
                                             gATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATccc             >  1:1730674/1‑69 (MQ=255)
                                             gATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATccc             >  1:1734533/1‑69 (MQ=255)
                                             gATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATccc             >  1:1853037/1‑69 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGTGCTAGATGGTTTTTGCCTGTCTCTATCCCAg           >  1:2750095/1‑69 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGtctc                   <  1:1044391/61‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGtctc                   <  1:165524/61‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           <  1:2901585/69‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           >  1:966129/1‑69 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           >  1:1254708/1‑69 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           >  1:827847/1‑69 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           <  1:1955829/69‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           <  1:1833094/69‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           <  1:391266/69‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAg           <  1:1693874/69‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCACAg           <  1:1485652/69‑1 (MQ=255)
                                               tGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCACGAGATGGTTTTTGCCTGTCTCTATCCCAg           <  1:3432003/69‑1 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGtt                               <  1:225736/48‑1 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGtt                               <  1:131603/48‑1 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:463796/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:453733/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:2881308/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:2589226/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:701274/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:2207180/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:1186764/1‑69 (MQ=255)
                                                gTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGc          >  1:3134436/1‑69 (MQ=255)
                                                   ttGCTGAAACCGGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATg       >  1:2564578/1‑69 (MQ=255)
                                                   ttGCTGAAACCGGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATg       >  1:2065715/1‑69 (MQ=255)
                                                   ttGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATg       >  1:304426/1‑69 (MQ=255)
                                                        gAAACCAGTTAAAGATCCGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTaa  <  1:3363705/69‑1 (MQ=255)
                                                              |                                                              
GAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAA  >  minE/842185‑842309

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: