Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,623,632 G→A 48.0% E229E (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,6320GA48.0% 10.0 / 89.2 102E229E (GAG→GAAftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base G (31/22);  new base A (28/21);  total (59/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

AGAAAATCTCGGTTCCGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACC  >  minE/2623566‑2623695
                                                                  |                                                               
aGAAAATCTCGGTTCCGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTgagga                                                               <  1:2294736/69‑1 (MQ=255)
              ccGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAag                                                       <  1:837343/63‑1 (MQ=255)
              ccGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAag                                                       <  1:355461/63‑1 (MQ=255)
              ccGGATTCATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAag                                                       <  1:3119025/63‑1 (MQ=255)
               cGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAag                                                       <  1:78789/62‑1 (MQ=255)
               cGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAag                                                       <  1:1561907/62‑1 (MQ=255)
                   tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                            <  1:2100042/69‑1 (MQ=255)
                   tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                            <  1:3033761/69‑1 (MQ=255)
                   tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                            <  1:2894371/69‑1 (MQ=255)
                   tttATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                            <  1:2293831/69‑1 (MQ=255)
                   tttATCAGCCTGTGCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTg                                            <  1:1671254/69‑1 (MQ=255)
                     tATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAag                                                       <  1:420540/56‑1 (MQ=255)
                         aGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                      <  1:121917/69‑1 (MQ=255)
                         aGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                      <  1:463969/69‑1 (MQ=255)
                         aGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGcc                                      <  1:156624/69‑1 (MQ=255)
                            cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                   >  1:2002535/1‑69 (MQ=255)
                            cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                   >  1:2788213/1‑69 (MQ=255)
                            cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                   >  1:1438208/1‑69 (MQ=255)
                            cTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                   >  1:590117/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:437462/1‑69 (MQ=255)
                                tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTgg                               >  1:1595835/1‑69 (MQ=255)
                                tCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTgg                               >  1:2928785/1‑69 (MQ=255)
                                        aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                       >  1:2845015/1‑53 (MQ=255)
                                        aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                       >  1:1968045/1‑53 (MQ=255)
                                        aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                       >  1:225417/1‑53 (MQ=255)
                                        aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                       >  1:1214477/1‑53 (MQ=255)
                                        aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                       >  1:842562/1‑53 (MQ=255)
                                        aaaaaaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGc                                       >  1:2462622/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:1973201/69‑1 (MQ=255)
                                            aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                   <  1:236337/69‑1 (MQ=255)
                                            aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                   <  1:147726/69‑1 (MQ=255)
                                            aaaTCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                   <  1:957601/69‑1 (MQ=255)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGAt                                   >  1:2475052/1‑43 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2845571/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:286605/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2937253/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2993965/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2999580/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2678697/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:3223331/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:524887/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:530891/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:988592/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1258448/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1249620/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1642928/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1249970/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2141025/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1155659/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1569402/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2000469/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:2072459/1‑48 (MQ=37)
                                                      tgatCTGTTTGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgggcg                            >  1:1186333/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:1681585/1‑65 (MQ=255)
                                                       gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa            >  1:2748434/1‑65 (MQ=255)
                                                       gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa            >  1:357903/1‑65 (MQ=255)
                                                       gatCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTaa            >  1:2340394/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)
                                                             tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACca                   <  1:2912659/52‑1 (MQ=255)
                                                             tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCAcc  >  1:1635859/1‑69 (MQ=255)
                                                             tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCAcc  >  1:1742845/1‑69 (MQ=255)
                                                             tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCAcc  >  1:1158165/1‑69 (MQ=255)
                                                             tttGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCAcc  >  1:2626072/1‑69 (MQ=255)
                                                             tttGAAGAGCTGGAAGAACATCTGTTGATCGCCGATGTgg                               <  1:5108/40‑1 (MQ=21)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2342942/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:3413363/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:3405038/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:125092/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:471829/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:1242001/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:1439706/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2366974/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2441016/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2726344/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2808830/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:203652/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2864660/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:1706627/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:1657870/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:1656814/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:2995933/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:1518068/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:3122100/40‑1 (MQ=25)
                                                             tttGAAGAGCTGGAAGAACAGCTGTTGATCGCCGATGTgg                               <  1:3404661/40‑1 (MQ=25)
                                                               tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc     >  1:425164/1‑64 (MQ=255)
                                                               tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc     >  1:1541322/1‑64 (MQ=255)
                                                               tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc     >  1:235628/1‑64 (MQ=255)
                                                               tGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATc     >  1:1001996/1‑64 (MQ=255)
                                                                  |                                                               
AGAAAATCTCGGTTCCGGATTTATCAGCCTGTTCCGCGGTAAAAAAATCGACGATGATCTGTTTGAGGAGCTGGAAGAGCAGCTTTTGATCGCCGATGTGGGTGTGGAAACCACACGTAAAATTATCACC  >  minE/2623566‑2623695

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: