Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,523,557 C→G 5.1% S53T (AGT→ACT)  dinD ← DNA‑damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,523,5570CG 5.1% 428.6 / 6.7 136S53T (AGT→ACT) dinDDNA‑damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base C (58/71);  new base G (5/2);  total (63/73)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

TCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAACCCAGTAG  >  minE/2523487‑2523627
                                                                      |                                                                      
 ctTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACt                                                                       >  1:2485837/1‑71 (MQ=255)
   tttGGCATTTGAGACAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGt                                                                     >  1:4030250/1‑71 (MQ=255)
       gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACTACTGTtct                                                                  <  1:2587336/70‑1 (MQ=255)
       gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                                  <  1:2798215/70‑1 (MQ=255)
       gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                                  <  1:1329977/70‑1 (MQ=255)
       gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                                  <  1:2127272/70‑1 (MQ=255)
        cATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                                  <  1:2966129/69‑1 (MQ=255)
           ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTTTTCTCACAg                                                             >  1:163541/1‑71 (MQ=255)
           ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTca                                                                >  1:163260/1‑68 (MQ=255)
           ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTca                                                                >  1:1555599/1‑68 (MQ=255)
           ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAg                                                             >  1:4233716/1‑71 (MQ=255)
            tGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAg                                                             >  1:3254531/1‑70 (MQ=255)
                ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                        <  1:3990991/71‑1 (MQ=255)
                ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                        <  1:4021030/71‑1 (MQ=255)
                ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                        <  1:4214359/71‑1 (MQ=255)
                ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                        <  1:847126/71‑1 (MQ=255)
                     tttGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGct                                                          >  1:3066302/1‑64 (MQ=255)
                      ttGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCCCAg                                                             >  1:3274290/1‑60 (MQ=255)
                            atatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGt                                            >  1:3358136/1‑71 (MQ=255)
                            atatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGt                                            <  1:3603310/71‑1 (MQ=255)
                            atatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGt                                            >  1:2915006/1‑71 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                                  <  1:1166583/48‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTca                                                                <  1:1680192/50‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTca                                                                <  1:1974258/50‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTca                                                                <  1:2722946/50‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                           <  1:3665626/71‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                           <  1:871809/71‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTg                                                                      <  1:3663983/44‑1 (MQ=255)
                             tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTg                                                                      <  1:384921/44‑1 (MQ=255)
                              atCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                           >  1:2084363/1‑70 (MQ=255)
                              atCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                           >  1:902989/1‑70 (MQ=255)
                               tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         >  1:3720820/1‑71 (MQ=255)
                               tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         <  1:3805777/71‑1 (MQ=255)
                               tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTa                                          <  1:3165929/70‑1 (MQ=255)
                               tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTa                                          <  1:3791605/70‑1 (MQ=255)
                               tCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTggt                                                                    >  1:2912782/1‑42 (MQ=25)
                               tCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTggt                                                                    >  1:2297072/1‑42 (MQ=25)
                                 gAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATa                                       >  1:1299488/1‑71 (MQ=255)
                                  aGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                      <  1:550264/71‑1 (MQ=255)
                                  aGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                      <  1:419819/71‑1 (MQ=255)
                                  aGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                      <  1:1193181/71‑1 (MQ=255)
                                   ggATCTCCTCGAAATTGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATa                                       <  1:2257815/69‑1 (MQ=255)
                                   ggATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATa                                       <  1:2162150/69‑1 (MQ=255)
                                    gATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                           <  1:1346808/64‑1 (MQ=255)
                                    gATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         <  1:3181341/66‑1 (MQ=255)
                                    gATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         <  1:12230/66‑1 (MQ=255)
                                     aTCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACa                                     >  1:1780398/1‑69 (MQ=255)
                                     aTCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACa                                     >  1:3113724/1‑69 (MQ=255)
                                     aTCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                                   >  1:3558895/1‑71 (MQ=255)
                                     aTCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAg                                    >  1:1021338/1‑70 (MQ=255)
                                       ctcctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttctt                                                      >  1:2319931/1‑50 (MQ=255)
                                       ctcctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         >  1:329188/1‑63 (MQ=255)
                                       ctcctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         >  1:253590/1‑63 (MQ=255)
                                         cctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttct                                                       <  1:1680596/47‑1 (MQ=255)
                                          ctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGtata                               >  1:3374354/1‑70 (MQ=255)
                                           tcGAAATTGTCATCAATTGTGTGACCACTGTTCTCACAGGcttctt                                                      <  1:1217846/46‑1 (MQ=255)
                                           tcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttctt                                                      <  1:1222873/46‑1 (MQ=255)
                                           tcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttctt                                                      <  1:3465829/46‑1 (MQ=255)
                                            cGAAATGGTCATCAATTGTGTGACCACTGTTCTcaca                                                              >  1:26830/1‑37 (MQ=255)
                                            cGAAATGGTCATCAATTGTGTGACCACTGTTCTcaca                                                              >  1:4189684/1‑37 (MQ=255)
                                             gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:1398288/71‑1 (MQ=255)
                                             gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           >  1:3941033/1‑71 (MQ=255)
                                             gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           >  1:3987722/1‑71 (MQ=255)
                                             gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           >  1:2417207/1‑71 (MQ=255)
                                              aaaTGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:1867601/70‑1 (MQ=255)
                                                aTGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGAc                        <  1:2441943/71‑1 (MQ=255)
                                                 tGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGAc                        <  1:1933135/70‑1 (MQ=255)
                                                 tGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGAc                        <  1:2985134/70‑1 (MQ=255)
                                                    tcatcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATa                    <  1:792880/71‑1 (MQ=255)
                                                     catcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                                   >  1:944369/1‑55 (MQ=255)
                                                      atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         >  1:2247523/1‑48 (MQ=255)
                                                      atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                      <  1:1620359/51‑1 (MQ=255)
                                                      atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                      <  1:2555963/51‑1 (MQ=255)
                                                      atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:1841637/62‑1 (MQ=255)
                                                      atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATa                    <  1:3803575/69‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                                    <  1:2214192/36‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                                    <  1:2504958/36‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                                    <  1:1297285/36‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                                    <  1:2547762/36‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                                    <  1:624801/36‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                                    <  1:1747065/36‑1 (MQ=255)
                                                       tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGtat                                <  1:455188/56‑1 (MQ=255)
                                                        caATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGa                         >  1:2153156/1‑62 (MQ=255)
                                                         aaTTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATa                    >  1:1710747/1‑66 (MQ=255)
                                                         aaTTGTGTCACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATa                    <  1:666143/66‑1 (MQ=255)
                                                          aTTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:1358470/58‑1 (MQ=255)
                                                            tgtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATaaa                             >  1:753328/1‑54 (MQ=255)
                                                             gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                         >  1:976541/1‑41 (MQ=255)
                                                             gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGa              >  1:3957409/1‑68 (MQ=255)
                                                             gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAAt            >  1:3975420/1‑70 (MQ=255)
                                                             gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAAt            >  1:1400988/1‑70 (MQ=255)
                                                             gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAAt            >  1:860899/1‑70 (MQ=255)
                                                              tgtgACCACTGTCCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATa           <  1:303231/70‑1 (MQ=255)
                                                                tgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAAcc        >  1:832968/1‑71 (MQ=255)
                                                                  aCCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                                   >  1:4186259/1‑42 (MQ=255)
                                                                    cACTGTTCTCACAGTCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:1600449/48‑1 (MQ=255)
                                                                    cACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:2411612/48‑1 (MQ=255)
                                                                    cACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                           <  1:26091/48‑1 (MQ=255)
                                                                     aCTGTTCTCAGAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATa                    <  1:2459227/54‑1 (MQ=255)
                                                                      cTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAACCCagtag  >  1:3663455/1‑71 (MQ=255)
                                                                      |                                                                      
TCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAACCCAGTAG  >  minE/2523487‑2523627

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: