Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,523,551 T→A 5.6% H55L (CAC→CTC)  dinD ← DNA‑damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,523,5510TA 5.6% 359.2 / 10.7 125H55L (CAC→CTC) dinDDNA‑damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base T (51/67);  new base A (5/2);  total (56/69)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 140 total aligned reads displayed.

ATGCTCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAACC  >  minE/2523483‑2523621
                                                                    |                                                                      
aTGCTCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATtgtgtg                                                                       <  1:121678/70‑1 (MQ=255)
 tGCTCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATtgtgtg                                                                       <  1:3270182/69‑1 (MQ=255)
 tGCTCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATtgtgtg                                                                       <  1:3941791/69‑1 (MQ=255)
   ctctTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCa                                                                   <  1:1372787/71‑1 (MQ=255)
   ctctTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCa                                                                   <  1:3760894/71‑1 (MQ=255)
    tctTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCAc                                                                  >  1:2175737/1‑71 (MQ=255)
    tctTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCAc                                                                  >  1:1627583/1‑71 (MQ=255)
     ctTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACt                                                                 >  1:2485837/1‑71 (MQ=255)
       tttGGCATTTGAGACAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGt                                                               >  1:4030250/1‑71 (MQ=255)
        ttGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGtt                                                              >  1:148163/1‑71 (MQ=255)
        ttGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGtt                                                              >  1:2509807/1‑71 (MQ=255)
           gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACTACTGTtct                                                            <  1:2587336/70‑1 (MQ=255)
           gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                            <  1:1329977/70‑1 (MQ=255)
           gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                            <  1:1395121/70‑1 (MQ=255)
           gCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                            <  1:2127272/70‑1 (MQ=255)
            cATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtct                                                            <  1:2966129/69‑1 (MQ=255)
               ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTtctc                                                           >  1:2270685/1‑67 (MQ=255)
               ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTca                                                          >  1:163260/1‑68 (MQ=255)
               ttGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAg                                                       >  1:4233716/1‑71 (MQ=255)
               ttGAGCCAATTTTGACCATATCGAGGATATCCTCGAAATGGTCATCAATTGTGTGACCACTGtt                                                              >  1:1942751/1‑64 (MQ=255)
                    ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                  <  1:847126/71‑1 (MQ=255)
                    ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                  <  1:4021030/71‑1 (MQ=255)
                    ccAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttc                                                  <  1:4214359/71‑1 (MQ=255)
                         tttGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGct                                                    >  1:3066302/1‑64 (MQ=255)
                          ttGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGa                                                                      >  1:515070/1‑45 (MQ=255)
                          ttGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCCCAg                                                       >  1:3274290/1‑60 (MQ=255)
                              ccATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGgcgcgc                                        <  1:3197965/71‑1 (MQ=255)
                                atatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGt                                      >  1:2915006/1‑71 (MQ=255)
                                atatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGt                                      >  1:3358136/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:2128866/71‑1 (MQ=255)
                                 tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                     <  1:871809/71‑1 (MQ=255)
                                 tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                     <  1:3665626/71‑1 (MQ=255)
                                 tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTg                                                                <  1:3663983/44‑1 (MQ=255)
                                 tatCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTg                                                                <  1:384921/44‑1 (MQ=255)
                                  atCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGtt                                     >  1:902989/1‑70 (MQ=255)
                                   tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                   <  1:3805777/71‑1 (MQ=255)
                                   tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                   >  1:3720820/1‑71 (MQ=255)
                                   tCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTa                                    <  1:3165929/70‑1 (MQ=255)
                                   tCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTggt                                                              >  1:3521858/1‑42 (MQ=25)
                                   tCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTggt                                                              >  1:2912782/1‑42 (MQ=25)
                                   tCGAGGATCTCCTCGAAATGGTCATCAATTGTGAGAACAGTgg                                                               >  1:1841657/1‑42 (MQ=37)
                                     gAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAAt                                  <  1:29999/70‑1 (MQ=255)
                                      aGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                <  1:1193181/71‑1 (MQ=255)
                                      aGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                <  1:550264/71‑1 (MQ=255)
                                      aGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                <  1:419819/71‑1 (MQ=255)
                                       ggATCTCCTCGAAATTGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATa                                 <  1:2257815/69‑1 (MQ=255)
                                       ggATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATa                                 <  1:2162150/69‑1 (MQ=255)
                                        gATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                   <  1:3181341/66‑1 (MQ=255)
                                        gATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                   <  1:12230/66‑1 (MQ=255)
                                         aTCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACa                               >  1:3113724/1‑69 (MQ=255)
                                         aTCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                             >  1:2395387/1‑71 (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:1295933/1‑63 (MQ=255)
                                           ctcctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                   >  1:253590/1‑63 (MQ=255)
                                              ctcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGtata                         >  1:3374354/1‑70 (MQ=255)
                                               tcGAAATTGTCATCAATTGTGTGACCACTGTTCTCACAGGcttctt                                                <  1:1217846/46‑1 (MQ=255)
                                               tcGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGcttctt                                                <  1:660133/46‑1 (MQ=255)
                                                cGAAATGGTCATCAATTGTGTGACCACTGTTCTcaca                                                        >  1:26830/1‑37 (MQ=255)
                                                cGAAATGGTCATCAATTGTGTGACCACTGTTCTcaca                                                        >  1:4189684/1‑37 (MQ=255)
                                                 gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                     >  1:3987722/1‑71 (MQ=255)
                                                 gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                     <  1:1321428/71‑1 (MQ=255)
                                                 gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                     >  1:3941033/1‑71 (MQ=255)
                                                 gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                     >  1:2417207/1‑71 (MQ=255)
                                                 gAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                     >  1:3798550/1‑71 (MQ=255)
                                                     tGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGAc                  <  1:2985134/70‑1 (MQ=255)
                                                     tGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGAc                  <  1:1933135/70‑1 (MQ=255)
                                                         catcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                             >  1:1745200/1‑55 (MQ=255)
                                                         catcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                             >  1:944369/1‑55 (MQ=255)
                                                         catcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAgg                             >  1:2726451/1‑55 (MQ=255)
                                                          atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTaa                                   >  1:2450713/1‑48 (MQ=255)
                                                          atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                <  1:2555963/51‑1 (MQ=255)
                                                          atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                <  1:1620359/51‑1 (MQ=255)
                                                          atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATAc                                <  1:4004280/51‑1 (MQ=255)
                                                          atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGt                     <  1:1841637/62‑1 (MQ=255)
                                                          atcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATa              <  1:3803575/69‑1 (MQ=255)
                                                           tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                              <  1:361276/36‑1 (MQ=255)
                                                           tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                              <  1:1297285/36‑1 (MQ=255)
                                                           tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                              <  1:2214192/36‑1 (MQ=255)
                                                           tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                              <  1:2547762/36‑1 (MQ=255)
                                                           tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTg                                              <  1:1747065/36‑1 (MQ=255)
                                                           tcaATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGtat                          <  1:455188/56‑1 (MQ=255)
                                                            caATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGa                   >  1:3870066/1‑62 (MQ=255)
                                                            caATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATAct            >  1:1642996/1‑69 (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:860899/1‑70 (MQ=255)
                                                                 gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAAt      >  1:3350386/1‑70 (MQ=255)
                                                                 gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAAt      >  1:1400988/1‑70 (MQ=255)
                                                                 gtgtgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAAt      >  1:3975420/1‑70 (MQ=255)
                                                                    tgACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAAcc  >  1:832968/1‑71 (MQ=255)
                                                                    |                                                                      
ATGCTCTTTTGGCATTTGAGCCAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCATCAATTGTGTGACCACTGTTCTCACAGGCTTCTTTGGCGCGCGTTAATACAGGTATAAAGTGACGATACTCTGAATAACC  >  minE/2523483‑2523621

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: