Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,623,644 G→A 42.9% E233E (GAG→GAA ftsY → fused Signal Recognition Particle (SRP) receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,623,6440GA42.9% 41.6 / 85.0 114E233E (GAG→GAAftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base G (44/21);  new base A (28/21);  total (72/42)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGGCGC  >  minE/2623585‑2623712
                                                           |                                                                    
tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                                             <  1:2894371/69‑1 (MQ=255)
tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                                             <  1:2293831/69‑1 (MQ=255)
tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                                             <  1:3033761/69‑1 (MQ=255)
tttATCAGCCTGTGCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                                             <  1:1671254/69‑1 (MQ=255)
      aGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                                       <  1:463969/69‑1 (MQ=255)
      aGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                                       <  1:156624/69‑1 (MQ=255)
      aGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                                       <  1:121917/69‑1 (MQ=255)
         cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                                    >  1:1438208/1‑69 (MQ=255)
         cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                                    >  1:2002535/1‑69 (MQ=255)
         cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                                    >  1:590117/1‑69 (MQ=255)
         cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                                    >  1:2788213/1‑69 (MQ=255)
         cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                                    >  1:1789734/1‑69 (MQ=255)
             tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGAtgtg                                                 >  1:3341415/1‑68 (MQ=255)
             tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTgg                                                >  1:3061184/1‑69 (MQ=255)
             tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTgg                                                >  1:1595835/1‑69 (MQ=255)
             tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTgg                                                >  1:2928785/1‑69 (MQ=255)
             tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTgg                                                >  1:437462/1‑69 (MQ=255)
                     aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                                        >  1:225417/1‑53 (MQ=255)
                     aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                                        >  1:1214477/1‑53 (MQ=255)
                     aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                                        >  1:842562/1‑53 (MQ=255)
                      aaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                                        >  1:2830007/1‑52 (MQ=255)
                      aaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGaa                                        >  1:2977405/1‑68 (MQ=255)
                      aaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATAGCCGATGTGGGTGTGGa                                         >  1:1640572/1‑67 (MQ=255)
                         aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                                    <  1:147726/69‑1 (MQ=255)
                         aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                                    <  1:957601/69‑1 (MQ=255)
                         aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                                    <  1:1973201/69‑1 (MQ=255)
                         aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                                    <  1:236337/69‑1 (MQ=255)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                                    >  1:2475052/1‑43 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2678697/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:286605/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2000469/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2845571/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2887329/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2141025/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2072459/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:1642928/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:530891/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:988592/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:1155659/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:1186333/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:1249970/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2937253/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:1258448/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:524887/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:1569402/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                                             >  1:2993965/1‑48 (MQ=37)
                                   tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCATGTgggc                                              >  1:1565781/1‑47 (MQ=21)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTTGGTGTGGAAACCACACGTaa                             >  1:1082915/1‑65 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa                             >  1:1304813/1‑65 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa                             >  1:357903/1‑65 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa                             >  1:2340394/1‑65 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa                             >  1:2748434/1‑65 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAAtt                         <  1:2183212/69‑1 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAAtt                         <  1:110895/69‑1 (MQ=255)
                                    gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAAt                          >  1:2735477/1‑68 (MQ=255)
                                     atCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                                       <  1:1003905/38‑1 (MQ=255)
                                          tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCAcc                   >  1:2626072/1‑69 (MQ=255)
                                          tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCAcc                   >  1:1742845/1‑69 (MQ=255)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:3122100/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:1518068/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:3404661/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:1656814/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:2995933/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:1657870/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:1706627/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:3405038/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:203652/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:1439706/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:2864660/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:471829/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:2366974/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:2808830/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:125092/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:2441016/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:1242001/40‑1 (MQ=25)
                                          tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                                                <  1:2726344/40‑1 (MQ=25)
                                            tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc                      >  1:425164/1‑64 (MQ=255)
                                            tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc                      >  1:1001996/1‑64 (MQ=255)
                                            tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc                      >  1:235628/1‑64 (MQ=255)
                                            tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc                      >  1:1769935/1‑64 (MQ=255)
                                                gagCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACTTAAAATTATCACCAATCt              <  1:307574/68‑1 (MQ=255)
                                                gagCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCt              <  1:3141327/68‑1 (MQ=255)
                                                gagCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCt              <  1:2604227/68‑1 (MQ=255)
                                                gagCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCt              <  1:2081303/68‑1 (MQ=255)
                                                  gCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGa            >  1:3263775/1‑68 (MQ=255)
                                                  gCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGa            >  1:1210186/1‑68 (MQ=255)
                                                  gCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGa            >  1:818364/1‑68 (MQ=255)
                                                  gCTGGAACAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGa            >  1:1756777/1‑68 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:1481676/1‑36 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:2401100/1‑36 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:1579593/1‑36 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:714075/1‑36 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:3264254/1‑36 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:1020204/1‑36 (MQ=255)
                                                     ggAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGa                                         >  1:2634925/1‑36 (MQ=255)
                                                          agcagcTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGgcg   <  1:709683/69‑1 (MQ=255)
                                                           gcagcTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGgcgc  >  1:1460880/1‑69 (MQ=255)
                                                           gcagcTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGgcgc  >  1:1367679/1‑69 (MQ=255)
                                                           gcagcTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGgcgc  >  1:759094/1‑69 (MQ=255)
                                                           gcagcTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGgcgc  >  1:3327199/1‑69 (MQ=255)
                                                           |                                                                    
TTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACCAATCTGACGGAAGGCGC  >  minE/2623585‑2623712

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: