Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 573,576 A→G 20.7% G135G (GGA→GGG ybjD → conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE573,5760AG20.7% 226.8 / 57.7 106G135G (GGA→GGGybjDconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (38/46);  new base G (5/17);  total (43/63)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGG  >  minE/573511‑573637
                                                                 |                                                             
cTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATcc                                                          >  1:3662240/1‑71 (MQ=255)
   gAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATcc                                                          <  1:1811338/68‑1 (MQ=255)
   gAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGAt                                                       >  1:2407953/1‑71 (MQ=255)
     aGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTg                                                     <  1:3484633/71‑1 (MQ=255)
          agagTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTc                                                >  1:2877887/1‑71 (MQ=255)
          agagTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTc                                                >  1:1352184/1‑71 (MQ=255)
           gagTGCGGAAGACGGCAGGGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGCTTGATGTCg                                               <  1:221604/71‑1 (MQ=255)
           gagTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCg                                               <  1:1346331/71‑1 (MQ=255)
           gagTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCg                                               <  1:63869/71‑1 (MQ=255)
           gagTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCg                                               <  1:2338866/71‑1 (MQ=255)
           gagTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCg                                               <  1:3849279/71‑1 (MQ=255)
            agTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGCCGa                                              <  1:3734750/71‑1 (MQ=255)
             gTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAg                                             <  1:722185/71‑1 (MQ=255)
              tGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGAt                                                       <  1:2149560/60‑1 (MQ=255)
              tGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGAt                                                       <  1:2386634/60‑1 (MQ=255)
                cGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGa                                           <  1:3715114/70‑1 (MQ=255)
                   aaGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGatat                                        >  1:2104044/1‑70 (MQ=255)
                   aaGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGatat                                        >  1:378705/1‑70 (MQ=255)
                   aaGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGatat                                        >  1:3400898/1‑70 (MQ=255)
                    aGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTa                                      <  1:1734768/71‑1 (MQ=255)
                    aGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTa                                      <  1:840138/71‑1 (MQ=255)
                      aCGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAAc                                    <  1:1347896/71‑1 (MQ=255)
                       cGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGATATTAACg                                   <  1:2883912/71‑1 (MQ=255)
                       cGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGATATTAACg                                   <  1:2776619/71‑1 (MQ=255)
                       cGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAAc                                    >  1:2637270/1‑70 (MQ=255)
                        ggCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGa                                  >  1:3398891/1‑71 (MQ=255)
                         gCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGAt                                 >  1:210652/1‑71 (MQ=255)
                         gCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGAt                                 >  1:3223848/1‑71 (MQ=255)
                         gCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGAt                                 <  1:1364521/71‑1 (MQ=255)
                         gCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGAt                                 <  1:2527020/71‑1 (MQ=255)
                           aGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGat                                          <  1:1463232/60‑1 (MQ=255)
                           aGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGat                                          <  1:490906/60‑1 (MQ=255)
                           aGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCa                               <  1:2769241/71‑1 (MQ=255)
                           aGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCa                               <  1:3304066/71‑1 (MQ=255)
                            gCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGat                                          <  1:2971255/59‑1 (MQ=255)
                            gCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAg                              <  1:3108405/71‑1 (MQ=255)
                            gCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAg                              >  1:2933774/1‑71 (MQ=255)
                             cGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCAATTGATGTCGAGGATATTaa                                     <  1:3912786/62‑1 (MQ=255)
                              gTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAgg                             >  1:1112228/1‑70 (MQ=255)
                               tgatgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:1211445/44‑1 (MQ=255)
                               tgatgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:1881657/44‑1 (MQ=255)
                               tgatgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:4121505/44‑1 (MQ=255)
                               tgatgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:487206/44‑1 (MQ=255)
                                gatgaCACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGa                                           >  1:3472624/1‑54 (MQ=255)
                                gatgaCACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCAc                          >  1:1867063/1‑71 (MQ=255)
                                 atgaCACTGCGCAGTTTTCTCGATAAAGACGGACATCCGa                                                        >  1:937613/1‑40 (MQ=255)
                                 atgaAACTGCGCAGTTTTCTCGATAAAGACGGACATCCGa                                                        >  1:1052164/1‑40 (MQ=255)
                                  tgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:3100915/41‑1 (MQ=255)
                                  tgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:3065666/41‑1 (MQ=255)
                                  tgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:2971298/41‑1 (MQ=255)
                                  tgaCACTGCGCAGTTTTCTCGATAAAGACGGGCATCCGAtt                                                      <  1:543631/41‑1 (MQ=255)
                                  tgaCACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGc                        >  1:2018564/1‑71 (MQ=255)
                                    acacTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGcc                       <  1:1496171/70‑1 (MQ=255)
                                     cacTGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCAt                     <  1:521971/71‑1 (MQ=255)
                                     cacTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGata                                         >  1:2905704/1‑51 (MQ=255)
                                     cacTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGata                                         >  1:885867/1‑51 (MQ=255)
                                     cacTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGata                                         >  1:3504220/1‑51 (MQ=255)
                                     cacTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCAt                     >  1:1082237/1‑71 (MQ=255)
                                      acTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGAtct                                                 >  1:586721/1‑40 (MQ=255)
                                        tGCGCAGTTTTCTCGATAAAGACGGGCATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATc                    <  1:3620141/69‑1 (MQ=255)
                                        tGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTg                  >  1:3814083/1‑71 (MQ=255)
                                         gcgcAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTgg                 >  1:1602490/1‑71 (MQ=255)
                                         gcgcAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTgg                 >  1:2870535/1‑71 (MQ=255)
                                          cgcAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATc                                >  1:1125548/1‑55 (MQ=255)
                                          cgcAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTg                  <  1:2673406/69‑1 (MQ=255)
                                            cAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGc              <  1:3158800/71‑1 (MQ=255)
                                               ttttCTCGATAAAGACGGACATCCGATTGATGTCGAGGa                                           <  1:1836661/39‑1 (MQ=255)
                                               ttttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAgg                             <  1:880918/53‑1 (MQ=255)
                                               ttttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGc                            <  1:75646/54‑1 (MQ=255)
                                               ttttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGc                            <  1:3173645/54‑1 (MQ=255)
                                               ttttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGtt           >  1:3207297/1‑71 (MQ=255)
                                                tttCTCGATAAAGACGGGCATCCGATTGATGTCGAGGATATTAACGAt                                 >  1:1877092/1‑48 (MQ=255)
                                                tttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGtt           >  1:785292/1‑70 (MQ=255)
                                                tttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGtt           >  1:123023/1‑70 (MQ=255)
                                                tttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGtt           >  1:1307848/1‑70 (MQ=255)
                                                 ttCTCGATAAAGACGGGCATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGttt          <  1:3296513/70‑1 (MQ=255)
                                                 ttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATc                    <  1:745234/60‑1 (MQ=255)
                                                 ttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATc                    <  1:711930/60‑1 (MQ=255)
                                                 ttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGttt          <  1:2189917/70‑1 (MQ=255)
                                                 ttCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGttt          >  1:10607/1‑70 (MQ=255)
                                                   ctcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGc                            <  1:2162895/50‑1 (MQ=255)
                                                   ctcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGc                            <  1:1356063/50‑1 (MQ=255)
                                                    tcGATAAAGACGGACATCCGATTGATGTCGAGGATAtt                                       <  1:2394780/38‑1 (MQ=255)
                                                    tcGATAAAGACGGACATCCGATTGATGTCGAGGATAtt                                       <  1:2437491/38‑1 (MQ=255)
                                                    tcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGttt          <  1:222621/67‑1 (MQ=255)
                                                    tcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAAt       >  1:2291869/1‑70 (MQ=255)
                                                    tcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAAt       >  1:1081895/1‑70 (MQ=255)
                                                    tcGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATg      >  1:2921238/1‑71 (MQ=255)
                                                      gATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGcc    <  1:3879224/71‑1 (MQ=255)
                                                       aTAAAGACGGGCATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCg   >  1:1728188/1‑71 (MQ=255)
                                                       aTAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAgg                             <  1:1752018/45‑1 (MQ=255)
                                                       aTAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAgg                             <  1:3584624/45‑1 (MQ=255)
                                                       aTAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTg                  <  1:594980/56‑1 (MQ=255)
                                                        tAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCgg  <  1:1029845/71‑1 (MQ=255)
                                                        tAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCgg  <  1:3228505/71‑1 (MQ=255)
                                                         aaaGACGGACATCCGATTGATGTCGAGGATATTAAGGATCAGGCACGcc                       >  1:3249167/1‑49 (MQ=255)
                                                         aaaGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGcc                       >  1:2122945/1‑49 (MQ=255)
                                                         aaaGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGt                >  1:3677792/1‑56 (MQ=255)
                                                           aGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTgg                 >  1:3067642/1‑53 (MQ=255)
                                                              cGGACATCCGATTGATGTCGAGGATATTAAGGATCAGGCACGCCATCTGGTg               >  1:2128629/1‑52 (MQ=255)
                                                                 |                                                             
CTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAGGCACGCCATCTGGTGCGTTTAATGCCGG  >  minE/573511‑573637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: