Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 5,581 G→A 9.8% intergenic (+51/+102) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE5,5810GA 9.8% 165.5 / 6.9 62intergenic (+51/+102)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/23);  new base A (3/3);  total (36/26)
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.86e-01

CAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGG  >  minE/5511‑5650
                                                                      |                                                                     
cAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTgg                                                                       >  1:4560764/1‑71 (MQ=255)
            aCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCa                                                           <  1:82975/71‑1 (MQ=255)
                  aaaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                     >  1:4630571/1‑71 (MQ=255)
                      gACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGt                                                  >  1:3674851/1‑70 (MQ=255)
                       aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGtt                                                 <  1:3212405/70‑1 (MQ=255)
                       aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTa                                                <  1:3452523/71‑1 (MQ=255)
                        cAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                               <  1:4412898/71‑1 (MQ=255)
                        cAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                               >  1:4630314/1‑71 (MQ=255)
                             gCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                                          >  1:1866947/1‑71 (MQ=255)
                                 ggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                                          <  1:910727/67‑1 (MQ=255)
                                 ggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                                          <  1:53032/67‑1 (MQ=255)
                                  gTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGTCGCGTCTTAt                                                             >  1:2487018/1‑47 (MQ=255)
                                  gTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTAt                                                             >  1:3831613/1‑47 (MQ=255)
                                   tAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGt                                                                  >  1:5568036/1‑41 (MQ=255)
                                     cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg                                  >  1:5540228‑M1/4‑53 (MQ=255)
                                     aCAATCAGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatatat                                  >  1:740861/1‑71 (MQ=255)
                                      cAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATAtt                                 <  1:2741661/71‑1 (MQ=255)
                                        aTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTg                                >  1:5468280/1‑70 (MQ=255)
                                        aTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGa                               <  1:4792640/71‑1 (MQ=255)
                                           cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAgg                                                         >  1:340860/1‑42 (MQ=255)
                                           cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATc                            >  1:4198018/1‑71 (MQ=255)
                                           cGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggat                            <  1:3323539‑M1/71‑25 (MQ=255)
                                            ggCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggat                            >  1:3154301‑M1/1‑46 (MQ=255)
                                             gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                          >  1:1033395/1‑71 (MQ=255)
                                              cATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                          >  1:3429588/1‑70 (MQ=255)
                                              cATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataag                         <  1:577759‑M1/71‑28 (MQ=255)
                                                ttCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt                           <  1:4489936/67‑1 (MQ=255)
                                                ttCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagc                        >  1:264260‑M1/1‑42 (MQ=255)
                                                ttCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccgg                              >  1:5316084‑M1/1‑42 (MQ=255)
                                                 tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                      <  1:5617861/39‑1 (MQ=255)
                                                 tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCAt                       >  1:1566045/1‑70 (MQ=255)
                                                   aacgcctGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaag                    >  1:2528465‑M1/3‑39 (MQ=255)
                                                   aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt                    >  1:626255/1‑71 (MQ=255)
                                                    acgccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttg                                      <  1:1843152‑M1/51‑15 (MQ=255)
                                                    acgccTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaag                    <  1:1266928‑M1/69‑33 (MQ=255)
                                                      gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgc                 <  1:270646‑M1/71‑36 (MQ=255)
                                                      gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcg                  >  1:5269286‑M1/1‑36 (MQ=255)
                                                      gcCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgc                 <  1:4307580‑M1/71‑36 (MQ=255)
                                                       ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgca                <  1:1028253‑M1/71‑37 (MQ=255)
                                                        cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgca                <  1:3871861‑M1/70‑37 (MQ=255)
                                                         tGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatc              >  1:4757524‑M1/1‑33 (MQ=255)
                                                         tGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccgg                              <  1:2817629‑M1/55‑23 (MQ=255)
                                                          gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg             >  1:2140600/1‑71 (MQ=255)
                                                          gATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgca                >  1:3802728‑M1/1‑32 (MQ=255)
                                                           aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg            >  1:2157516/1‑71 (MQ=255)
                                                            tGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtag                                   >  1:1294376‑M1/1‑30 (MQ=255)
                                                            tGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtag                                   >  1:1451296‑M1/1‑30 (MQ=255)
                                                            tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg            >  1:3454739/1‑70 (MQ=255)
                                                             gCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatccgg           >  1:1124834‑M1/1‑29 (MQ=255)
                                                              cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg         >  1:4258947/1‑71 (MQ=255)
                                                              cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg         >  1:4075143/1‑71 (MQ=255)
                                                              cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg         <  1:5461251/71‑1 (MQ=255)
                                                              cGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggca         >  1:5385333‑M1/1‑28 (MQ=255)
                                                               gACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatccggcat        <  1:3128084‑M1/71‑45 (MQ=255)
                                                               gACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatccggcat        <  1:3170691‑M1/71‑45 (MQ=255)
                                                               gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTg                                <  1:697496/47‑1 (MQ=255)
                                                               gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg         >  1:3075199/1‑70 (MQ=255)
                                                                 cGCTGGCGCGTCTTATCAGTCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt      <  1:2522751/71‑1 (MQ=255)
                                                                  gCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacg       >  1:781071‑M1/1‑24 (MQ=255)
                                                                   cTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacg       <  1:4337157‑M1/68‑46 (MQ=255)
                                                                   cTGGCGCGTCTTATCAGGCCTACacgctgcgattttgta                                    >  1:297465‑M1/1‑23 (MQ=255)
                                                                   cTGGCGCGTCTTATCAGGCCTACacgctgcgattttgta                                    >  1:3077778‑M1/1‑23 (MQ=255)
                                                                    tGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                          <  1:4933170/48‑1 (MQ=255)
                                                                      gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                        >  1:4682794/1‑48 (MQ=255)
                                                                      gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg  <  1:2819361/70‑1 (MQ=255)
                                                                      |                                                                     
CAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGG  >  minE/5511‑5650

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: