Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,063,055 T→C 5.8% R412R (CGT→CGC)  dnaG → DNA primase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,063,0550TC 5.8% 276.2 / 6.4 97R412R (CGT→CGC) dnaGDNA primase
Reads supporting (aligned to +/- strand):  ref base T (41/49);  new base C (3/2);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TTTGAGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCAAAAGCGGCA  >  minE/2062989‑2063124
                                                                  |                                                                     
tttGAGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAgg                                                                   <  1:1211564/71‑1 (MQ=255)
 ttGAGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAgg                                                                   >  1:1830733/1‑70 (MQ=255)
    aGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAt                                                                <  1:2401029/70‑1 (MQ=255)
    aGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAt                                                                <  1:2959252/70‑1 (MQ=255)
    aGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAt                                                                <  1:3245104/70‑1 (MQ=255)
     gTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTa                                                              <  1:93049/71‑1 (MQ=255)
     gTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTa                                                              <  1:610248/71‑1 (MQ=255)
     gTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGCCAg                                                                    >  1:1061491/1‑65 (MQ=255)
     gTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGCCAg                                                                    >  1:2208078/1‑65 (MQ=255)
             gCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaaca                                                       >  1:1290932/1‑70 (MQ=255)
               aCTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGCCAGGAATTAGGCAACAAAt                                                    >  1:1070862/1‑71 (MQ=255)
                cTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAgg                                                                   <  1:2599084/55‑1 (MQ=255)
                cTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAgg                                                                   <  1:1306394/55‑1 (MQ=255)
                 tACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTa                                                              <  1:2890539/59‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAt                                                                >  1:2611699/1‑55 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAt                                                                <  1:2665726/55‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAt                                                                <  1:732150/55‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACATATTAgg                                                <  1:303/71‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                <  1:248253/71‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                <  1:287740/71‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                >  1:31503/1‑71 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                >  1:1205025/1‑71 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                >  1:1650490/1‑71 (MQ=255)
                   ccATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGCCAGGAATTAGGCAACAAATTAg                                                 <  1:1678163/70‑1 (MQ=255)
                   ccATTGATATCGCAAGTGCC‑GGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                <  1:2251506/70‑1 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                >  1:381156/1‑69 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              <  1:598685/71‑1 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              <  1:2498030/71‑1 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              >  1:3037418/1‑71 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              >  1:3041332/1‑71 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              <  1:1822705/71‑1 (MQ=255)
                     aTTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGCCAGGAAt                                                                <  1:267087/53‑1 (MQ=255)
                      ttGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              <  1:577777/70‑1 (MQ=255)
                        gATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATAc                                           >  1:2240192/1‑71 (MQ=255)
                         atatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAg                                                 <  1:504129/64‑1 (MQ=255)
                         atatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAg                                                 <  1:1497176/64‑1 (MQ=255)
                         atatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATAc                                           <  1:582762/70‑1 (MQ=255)
                         atatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGACAGGAATTAGGCAACAAATTAGGCATACt                                          >  1:2670870/1‑71 (MQ=255)
                         atatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGACAGGAATTAGGCAACAAATTAGGCATACt                                          >  1:3296410/1‑71 (MQ=255)
                          tatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAgg                                                            >  1:2203800/1‑52 (MQ=255)
                          tatCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAgg                                                            >  1:3154424/1‑52 (MQ=255)
                            tCGCAAGTGCCGTGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACt                                          >  1:983148/1‑68 (MQ=255)
                            tCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACt                                          >  1:2271559/1‑68 (MQ=255)
                            tCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtg                                        <  1:2829664/70‑1 (MQ=255)
                              gCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAtt                                                               >  1:1210736/1‑45 (MQ=255)
                              gCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtgat                                      <  1:1456929/70‑1 (MQ=255)
                              gCAAGTGCCGGGCGAAACGCTGAGAATATATCTTCGTCAGGAAt                                                                >  1:274621/1‑44 (MQ=255)
                                aaGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtgatga                                    >  1:3131337/1‑70 (MQ=255)
                                aaGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtgatga                                    >  1:586435/1‑70 (MQ=255)
                                aaGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtgat                                      >  1:2722427/1‑68 (MQ=255)
                                aaGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtga                                       >  1:2541059/1‑67 (MQ=255)
                                 aGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTa                                                  >  1:2786421/1‑55 (MQ=255)
                                 aGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                                                <  1:1819284/57‑1 (MQ=255)
                                 aGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGAca                                  <  1:2913922/71‑1 (MQ=255)
                                 aGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGAca                                  <  1:2180424/71‑1 (MQ=255)
                                  gTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACtt                                         <  1:473701/63‑1 (MQ=255)
                                   tGCCGGGCGAAACTCTGCGAATATATCTTCGTCAGGAAt                                                                >  1:2748794/1‑39 (MQ=255)
                                   tGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTa                                                              >  1:514558/1‑41 (MQ=255)
                                   tGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCa                                              >  1:356688/1‑57 (MQ=255)
                                    gCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGaa                                                                 >  1:966085/1‑37 (MQ=255)
                                    gCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGAcagc                                >  1:1936749/1‑70 (MQ=255)
                                          cGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTg                         <  1:2344526/71‑1 (MQ=255)
                                          cGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTg                         <  1:3172945/71‑1 (MQ=255)
                                           gAAACGCTTCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGa                        <  1:2800985/71‑1 (MQ=255)
                                           gAAACGCTGCGAATATATCTTCGTCATGAATTAGGCAACAAATTAGGCATACTtg                                        <  1:484695/55‑1 (MQ=255)
                                           gAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtg                                        <  1:1639839/55‑1 (MQ=255)
                                            aaaCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGaa                       >  1:2465205/1‑71 (MQ=255)
                                            aaaCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGaa                       >  1:2269957/1‑71 (MQ=255)
                                             aaCGCTGCGAATATATCTTCGTCAGGAATTAGGcaacaa                                                      <  1:2232271/39‑1 (MQ=255)
                                             aaCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATAc                                           >  1:1806788/1‑50 (MQ=255)
                                             aaCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtg                                        >  1:3316922/1‑53 (MQ=255)
                                              aCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAAc                      >  1:931666/1‑70 (MQ=255)
                                              aCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACg                     <  1:424892/71‑1 (MQ=255)
                                              aCGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACg                     <  1:457671/71‑1 (MQ=255)
                                               cGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAAt                                                    >  1:2946678/1‑39 (MQ=255)
                                               cGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAAt                                                    >  1:3046404/1‑39 (MQ=255)
                                                gCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAt                   <  1:2864799/71‑1 (MQ=255)
                                                 cTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:279105/71‑1 (MQ=255)
                                                 cTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:848337/71‑1 (MQ=255)
                                                 cTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:2256754/71‑1 (MQ=255)
                                                 cTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:3151912/71‑1 (MQ=255)
                                                 cTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAAGCATACTTGATGACAGCCAGCTTGAACGAt                   >  1:3014939/1‑70 (MQ=255)
                                                   gCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTaa                >  1:2197349/1‑71 (MQ=255)
                                                   gCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTaa                >  1:2207697/1‑71 (MQ=255)
                                                   gCGAATATATCATCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTaa                >  1:1505044/1‑71 (MQ=255)
                                                    cGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCAt                                             <  1:54278/41‑1 (MQ=255)
                                                        tatatCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGa                    >  1:3173979/1‑62 (MQ=255)
                                                        tatatCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAAt               <  1:498862/67‑1 (MQ=255)
                                                         atatCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTaa                <  1:2785203/65‑1 (MQ=255)
                                                          tatCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtga                                       >  1:1640585/1‑41 (MQ=255)
                                                          tatCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCaaa         <  1:3081810/71‑1 (MQ=255)
                                                             cTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:2291586/59‑1 (MQ=255)
                                                             cTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:689605/59‑1 (MQ=255)
                                                             cTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGAtt                  <  1:1351861/59‑1 (MQ=255)
                                                             cTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCAAAAGc      >  1:124553/1‑71 (MQ=255)
                                                             cTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCAAAAGc      >  1:3207478/1‑71 (MQ=255)
                                                                 gTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCAAAAGCGGCa  <  1:1635927/71‑1 (MQ=255)
                                                                  tCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCAAAAGCGGCa  <  1:3319299/70‑1 (MQ=255)
                                                                  |                                                                     
TTTGAGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACAGCCAGCTTGAACGATTAATGCCAAAAGCGGCA  >  minE/2062989‑2063124

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: ‑