Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 5,628 A→G 23.5% intergenic (+98/+55) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE5,6280AG23.5% 176.6 / 54.5 102intergenic (+98/+55)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/33);  new base G (11/13);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 103 total aligned reads displayed.

TCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGT  >  minE/5560‑5696
                                                                    |                                                                    
tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCAt                                                                     >  1:1566045/1‑70 (MQ=255)
  aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt                                                                  >  1:626255/1‑71 (MQ=255)
         gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg                                                           >  1:2140600/1‑71 (MQ=255)
          aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                          >  1:2157516/1‑71 (MQ=255)
           tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                          >  1:3454739/1‑70 (MQ=255)
             cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                       <  1:5461251/71‑1 (MQ=255)
             cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                       >  1:4258947/1‑71 (MQ=255)
              gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                                                       >  1:3075199/1‑70 (MQ=255)
                cGCTGGCGCGTCTTATCAGTCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                                                    <  1:2522751/71‑1 (MQ=255)
                     gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                                                                      >  1:4682794/1‑48 (MQ=255)
                     gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg                                                <  1:2819361/70‑1 (MQ=255)
                        cgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                            <  1:2253435/71‑1 (MQ=255)
                          tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa                                                   <  1:1069708/62‑1 (MQ=255)
                            ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                              >  1:467451/1‑49 (MQ=255)
                             tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGt                                                             >  1:5188918/1‑49 (MQ=255)
                             tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc                                       >  1:3893917/1‑71 (MQ=255)
                                 aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc                                    <  1:2229800/70‑1 (MQ=255)
                                 aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc                                    <  1:2775898/70‑1 (MQ=255)
                                 aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCa                                   >  1:4362773/1‑71 (MQ=255)
                                   gCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                          <  1:3846967/62‑1 (MQ=255)
                                     cTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATcc                                <  1:3393367/70‑1 (MQ=255)
                                      tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa                                                   <  1:5540013/50‑1 (MQ=255)
                                      tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                              >  1:3953771‑M1/1‑71 (MQ=255)
                                      tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                              >  1:5297533‑M1/1‑71 (MQ=255)
                                       aCGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                              >  1:3297847‑M1/1‑70 (MQ=255)
                                        cGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                            >  1:3059320‑M1/1‑71 (MQ=255)
                                        cGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                            >  1:2783007‑M1/1‑71 (MQ=255)
                                        cGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                            >  1:153691‑M1/1‑71 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                              >  1:2783382/1‑36 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                              >  1:3526172/1‑36 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                              >  1:1939379/1‑36 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTg                                                              <  1:3792198/35‑1 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:1288116/71‑1 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          >  1:3747856/1‑71 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:3994012/71‑1 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          >  1:1654681/1‑71 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          >  1:2323434/1‑71 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          >  1:2146278/1‑71 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:4113015/71‑1 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:1825340/71‑1 (MQ=255)
                                           tAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          >  1:1791558/1‑70 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTa                                                            >  1:2031531/1‑35 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTa                                                            >  1:154/1‑35 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAgg                                                      >  1:344469/1‑41 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAgg                                                      >  1:3636197/1‑41 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTg                         >  1:5380938/1‑70 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGa                        >  1:5387744/1‑71 (MQ=255)
                                             aTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGC                             <  1:1566591‑M1/65‑1 (MQ=255)
                                              ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTc                                           >  1:4311763/1‑50 (MQ=255)
                                              ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgcc                                      <  1:79045/55‑1 (MQ=255)
                                              ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgcc                                      <  1:2981545/55‑1 (MQ=255)
                                              ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                                 >  1:427967/1‑60 (MQ=255)
                                                cTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                                  >  1:3338352/1‑41 (MQ=255)
                                                cTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                                  >  1:691059/1‑41 (MQ=255)
                                                cTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                                  >  1:5620756/1‑41 (MQ=255)
                                                 tGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataa                     <  1:5640649‑M1/69‑7 (MQ=255)
                                                 tGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataa                     <  1:162167‑M1/69‑7 (MQ=255)
                                                 tGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaat                       <  1:1857813‑M1/67‑5 (MQ=255)
                                                  gCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                                                    <  1:388539/37‑1 (MQ=255)
                                                   cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcg                 <  1:3860884‑M1/71‑11 (MQ=255)
                                                   cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcg                 >  1:5612008‑M1/1‑61 (MQ=255)
                                                   cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaa                    <  1:5534453‑M1/68‑8 (MQ=255)
                                                   cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaa                    <  1:2088042‑M1/68‑8 (MQ=255)
                                                    aaTATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgc                <  1:4347374‑M1/71‑12 (MQ=255)
                                                     atatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgca               <  1:4381488‑M1/71‑13 (MQ=255)
                                                     atatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaat                       <  1:3884668‑M1/63‑5 (MQ=255)
                                                     atatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                            <  1:1848515/42‑1 (MQ=255)
                                                     atatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                            <  1:3812902/42‑1 (MQ=255)
                                                      tatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATga                         >  1:5604673‑M1/1‑58 (MQ=255)
                                                      tatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATga                         >  1:2398733‑M1/1‑58 (MQ=255)
                                                      tatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa                   <  1:110977/66‑1 (MQ=255)
                                                      tatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa                   <  1:4012164/66‑1 (MQ=255)
                                                      tatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACt               >  1:284536/1‑70 (MQ=255)
                                                       atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                            <  1:1638927/40‑1 (MQ=255)
                                                       atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                            <  1:2689296/40‑1 (MQ=255)
                                                       atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                          >  1:3930461/1‑42 (MQ=255)
                                                       atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                          >  1:1252321/1‑42 (MQ=255)
                                                        tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactc            <  1:2571599‑M1/71‑16 (MQ=255)
                                                        tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact             >  1:637586‑M1/1‑56 (MQ=255)
                                                        tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact             >  1:1579168‑M1/1‑56 (MQ=255)
                                                        tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact             >  1:1537235‑M1/1‑56 (MQ=255)
                                                        tatTGAATCTGCATGCTTTTGTAGGAAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTaa            >  1:5696364/1‑71 (MQ=255)
                                                          ttGAATTTGCGTGCTTTTGTAGGCAGGTTAAGGCGTTCACGCCGCATCCGGCATgaataaagcgca               >  1:579422‑M1/1‑54 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAAc           >  1:1645072/1‑70 (MQ=255)
                                                           tGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaa          >  1:526587‑M1/1‑53 (MQ=255)
                                                           tGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaa          >  1:1531146‑M1/1‑53 (MQ=255)
                                                             aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaacaa       >  1:4185550‑M1/1‑51 (MQ=255)
                                                             aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaata                      <  1:4928258‑M1/56‑6 (MQ=255)
                                                             aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaata                      <  1:1233799‑M1/56‑6 (MQ=255)
                                                             aaTCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctg       >  1:3503473/1‑71 (MQ=255)
                                                                 tGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaacaatct    >  1:1513021‑M1/1‑47 (MQ=255)
                                                                 tGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaacaatct    <  1:5064246‑M1/70‑24 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:3818629/48‑1 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:420100/48‑1 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                          <  1:3312909/48‑1 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACtt              <  1:2577814/60‑1 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACtt              <  1:2492083/60‑1 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCg   >  1:5346091/1‑71 (MQ=255)
                                                                  gCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATcc                                <  1:5001401/41‑1 (MQ=255)
                                                                   cATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGt  <  1:1482905/70‑1 (MQ=255)
                                                                    |                                                                    
TCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGT  >  minE/5560‑5696

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: