Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,651,936 C→T 15.1% A57V (GCG→GTG)  glgC → glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,651,9360CT15.1% 140.7 / 14.3 53A57V (GCG→GTG) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/15);  new base T (3/5);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTC  >  minE/2651866‑2652004
                                                                      |                                                                    
aaGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgc                                                                      <  1:22006/71‑1 (MQ=255)
 aGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgcg                                                                     >  1:971541/1‑71 (MQ=255)
       tAACCAATAAGCGAGCAAAACCAGCCGTACACTTTGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTc                                                                >  1:1271472/1‑70 (MQ=255)
        aaCCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCt                                                               >  1:87949/1‑70 (MQ=255)
           cAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgcg                                                                     <  1:442735/61‑1 (MQ=255)
           cAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgcg                                                                     <  1:965176/61‑1 (MQ=255)
           cAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACt                                                           >  1:241223/1‑71 (MQ=255)
           cAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACt                                                           >  1:618241/1‑71 (MQ=255)
                aGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCAt                                                       >  1:1873987/1‑70 (MQ=255)
                           ccGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTaa                                                             >  1:99543/1‑53 (MQ=255)
                            cGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATc                                          <  1:462849/71‑1 (MQ=255)
                                 gTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCTCTGTCTAACTGCATCAACTCCGGGATCcgtcg                                     <  1:1496148/71‑1 (MQ=255)
                                 gTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATc                                                      <  1:100056/54‑1 (MQ=255)
                                 gTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATc                                                      <  1:1698225/54‑1 (MQ=255)
                                 gTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcg                                     <  1:969690/71‑1 (MQ=255)
                                  tACACTTCGGCGGTAAGTTCCGCATTATCGACTTTgc                                                                      <  1:1496901/37‑1 (MQ=255)
                                   acacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCaa                                                    >  1:642572/1‑54 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGTGCTGTCTAAc                                                            >  1:1917790/1‑45 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACt                                                           >  1:342725/1‑46 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACt                                                           >  1:1271905/1‑46 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                                   <  1:1967721/70‑1 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                                   <  1:1353666/70‑1 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                                  <  1:827145/71‑1 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                                  <  1:1725065/71‑1 (MQ=255)
                                    cacTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                                  <  1:1967747/71‑1 (MQ=255)
                                          ggcggTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGt                            >  1:555663/1‑71 (MQ=255)
                                           gcggTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTg                           <  1:159919/71‑1 (MQ=255)
                                            cggTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATg                                 >  1:565539/1‑64 (MQ=255)
                                            cggTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATg                                 >  1:1765545/1‑64 (MQ=255)
                                              gTAAGTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTAt                                  <  1:817757/61‑1 (MQ=255)
                                               tAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCa                       >  1:1291389/1‑71 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGa                          >  1:347578/1‑67 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCAc                      >  1:1865985/1‑71 (MQ=255)
                                                aaGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCAc                      >  1:123722/1‑71 (MQ=255)
                                                 aGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCAcc                     <  1:618668/71‑1 (MQ=255)
                                                  gTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcgtcgt                                    <  1:1684846/55‑1 (MQ=255)
                                                  gTTCCGCATTATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCcgtcgt                                    <  1:445937/55‑1 (MQ=255)
                                                  gTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtc                                      <  1:2170352/53‑1 (MQ=255)
                                                        cATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGt                            <  1:1064715/57‑1 (MQ=255)
                                                         aTTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcat                                    >  1:666027/1‑46 (MQ=255)
                                                         aTTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcat                                    >  1:663513/1‑46 (MQ=255)
                                                          ttATCGACTTTGTGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGt                 <  1:1567905/66‑1 (MQ=255)
                                                           tATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTa                                   >  1:1557223/1‑47 (MQ=255)
                                                            aTCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc          >  1:1725831/1‑71 (MQ=255)
                                                            aTCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc          >  1:664178/1‑71 (MQ=255)
                                                            aTCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCAACCAGTACCAGTc          >  1:1688373/1‑71 (MQ=255)
                                                              cGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtc                                      >  1:1101894/1‑41 (MQ=255)
                                                              cGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtc                                      >  1:1497430/1‑41 (MQ=255)
                                                              cGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCcgtc                                      >  1:1921538/1‑41 (MQ=255)
                                                              cGACTTTGCGCTGTCTAACTGCATCAACTCAGGGATCcgtc                                      >  1:1844160/1‑41 (MQ=255)
                                                                aCTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc          >  1:1891387/1‑67 (MQ=255)
                                                                aCTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc          >  1:156671/1‑67 (MQ=255)
                                                                aCTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTc          >  1:1893915/1‑67 (MQ=255)
                                                                  tttGCGCTGTCTAACTGCATCAACTCCGGGATCcgtcgt                                    <  1:1862464/39‑1 (MQ=255)
                                                                    tGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCa                       >  1:779260/1‑50 (MQ=255)
                                                                    tGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCa                       >  1:1640009/1‑50 (MQ=255)
                                                                    tGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCacac    <  1:389411/69‑1 (MQ=255)
                                                                    tGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACActc  >  1:855974/1‑71 (MQ=255)
                                                                      |                                                                    
AAGGATTTAACCAATAAGCGAGCAAAACCGGCCGTACACTTCGGCGGTAAGTTCCGCATTATCGACTTTGCGCTGTCTAACTGCATCAACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTC  >  minE/2651866‑2652004

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: