Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,867,615 G→A 37.3% T429I (ACC→ATC)  ygcU ← predicted FAD containing dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,867,6150GA37.3% 64.1 / 70.6 67T429I (ACC→ATC) ygcUpredicted FAD containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (14/11);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTACAGTC  >  minE/1867552‑1867682
                                                               |                                                                   
cTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTttcttc                                                              <  1:491824/71‑1 (MQ=255)
cTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttctt                                                               <  1:957227/70‑1 (MQ=255)
   cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTca                                                             >  1:1118458/1‑69 (MQ=255)
     gTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTttctt                                                               >  1:1918189/1‑65 (MQ=255)
     gTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTttctt                                                               >  1:516821/1‑65 (MQ=255)
     gTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAg                                                          <  1:1428153/70‑1 (MQ=255)
     gTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAg                                                          <  1:1549281/70‑1 (MQ=255)
     gTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCCAATGGTttctt                                                               >  1:1231791/1‑65 (MQ=255)
      tGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTca                                                             <  1:1200003/66‑1 (MQ=255)
       gAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAgat                                                        >  1:955910/1‑70 (MQ=255)
        aaCGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAgatga                                                      >  1:369690/1‑71 (MQ=255)
        aaCGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAgatga                                                      >  1:445121/1‑71 (MQ=255)
            cgATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCtt                                                  <  1:697156/71‑1 (MQ=255)
            cgATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCtt                                                  <  1:948599/71‑1 (MQ=255)
               tGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCttgtt                                               >  1:1521349/1‑71 (MQ=255)
               tGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgtt                                               >  1:779445/1‑71 (MQ=255)
               tGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgtt                                               >  1:431716/1‑71 (MQ=255)
                gTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTttc                                                                 >  1:1915810/1‑52 (MQ=255)
                gTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTttc                                                                 >  1:1670283/1‑52 (MQ=255)
                   tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAga                                                         >  1:736781/1‑57 (MQ=255)
                   tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGc                                           >  1:77034/1‑71 (MQ=255)
                   tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTATTCACAga                                                         >  1:579975/1‑57 (MQ=255)
                   tACCGAGACCGTGGTGGTGCACCATCTAACCGCCGAGGCGAATGGTTTCTTCACAga                                                         >  1:2073/1‑57 (MQ=255)
                         tACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATtg                                     <  1:1084745/71‑1 (MQ=255)
                         tACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATtg                                     <  1:751241/71‑1 (MQ=255)
                              tggtggtgCACCATCGAACCGCCGAGGCGAATGAtt                                                                   <  1:1897291/36‑1 (MQ=255)
                                 tggtgCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCt                                                   <  1:318249/49‑1 (MQ=255)
                                    tgCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTc                          <  1:1535203/71‑1 (MQ=255)
                                    tgCACCATCCAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTc                          <  1:831215/71‑1 (MQ=255)
                                        ccATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAtt                      >  1:1819259/1‑71 (MQ=255)
                                        ccATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAtt                      >  1:567540/1‑71 (MQ=255)
                                        ccATCGAACCGCCGAGGCGAAAGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAtt                      <  1:277575/71‑1 (MQ=255)
                                         cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAttt                     <  1:1248316/71‑1 (MQ=255)
                                         cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAttt                     <  1:2057033/71‑1 (MQ=255)
                                         cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAtt                      <  1:1853991/70‑1 (MQ=255)
                                          aTCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCgg                                         <  1:199399/50‑1 (MQ=255)
                                          aTCGAACCGCCGAGGCGAATGATTTCTTCAAAgatga                                                      >  1:300472/1‑37 (MQ=38)
                                           tCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTca                   <  1:1293160/71‑2 (MQ=255)
                                            cGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCct                  >  1:2153056/1‑71 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttg                                                 >  1:323393/1‑39 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttg                                                 >  1:20253/1‑39 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttg                                                 >  1:525936/1‑39 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAt                                       >  1:297684/1‑49 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAt                                       >  1:1179024/1‑49 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAtt                      >  1:2000133/1‑66 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCctc                 >  1:202500/1‑71 (MQ=255)
                                             gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCctc                 >  1:892461/1‑71 (MQ=255)
                                              aaCCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCctct                <  1:755960/71‑1 (MQ=255)
                                              aaCCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCctct                <  1:474154/71‑1 (MQ=255)
                                                     gAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCt         <  1:655070/71‑1 (MQ=255)
                                                      aGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCtt        <  1:1832866/71‑1 (MQ=255)
                                                      aGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCtt        <  1:194005/71‑1 (MQ=255)
                                                      aGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAtt                      <  1:502917/57‑1 (MQ=255)
                                                       ggCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCtt        <  1:2128998/70‑1 (MQ=255)
                                                       ggCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTa       <  1:1684332/71‑1 (MQ=255)
                                                        gCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTa       >  1:1752670/1‑70 (MQ=255)
                                                         cGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATtgtg                                   <  1:2203937/41‑1 (MQ=255)
                                                         cGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTc                    >  1:985683/1‑56 (MQ=255)
                                                           aaTGGTTTCTTCACAGATGATCTTGTTGAGCGGATtgt                                    <  1:628518/38‑1 (MQ=255)
                                                           aaTGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTg                            <  1:56447/46‑1 (MQ=255)
                                                           aaTGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTACAGt   <  1:1880458/71‑1 (MQ=255)
                                                           aaTGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTACAGt   <  1:1722814/71‑1 (MQ=255)
                                                           aaTGATTTCTTCACAGATGATCTTGTTGAGCGGATtgt                                    <  1:1950874/38‑1 (MQ=255)
                                                            aTGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTACAGt   >  1:735569/1‑70 (MQ=255)
                                                            aTGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTACAGTc  <  1:1069273/71‑1 (MQ=255)
                                                             tGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTAc                               >  1:1868810/1‑41 (MQ=255)
                                                             tGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTAGTAc                               >  1:153030/1‑41 (MQ=255)
                                                               |                                                                   
CTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTACAGTC  >  minE/1867552‑1867682

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