Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,664,529 G→A I231I (ATC→ATT guaA ← GMP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,664,5290GA97.4% 113.4 / ‑4.0 39I231I (ATC→ATTguaAGMP synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/22);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGC  >  minE/1664462‑1664595
                                                                   |                                                                  
ttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAAt                                                                  >  1:1179355/1‑70 (MQ=255)
ttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAAt                                                                  >  1:857232/1‑70 (MQ=255)
ttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAAt                                                                  >  1:258075/1‑70 (MQ=255)
     ccGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACttt                                                            <  1:1923872/71‑1 (MQ=255)
        aTAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTAtc                                                         >  1:2412393/1‑71 (MQ=255)
         tAGCGCGGTGCAGCAGCATTGGGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTAtcg                                                        <  1:348016/71‑1 (MQ=255)
              cgGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGcc                                                   <  1:2430146/71‑1 (MQ=255)
               ggTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGcc                                                   >  1:437908/1‑70 (MQ=255)
                    agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                             >  1:1469050/1‑71 (MQ=255)
                     gcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGc                                            <  1:1064735/71‑1 (MQ=255)
                         caTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTcgc                                        <  1:2720418/71‑1 (MQ=255)
                           ttGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCgg                                      <  1:150049/71‑1 (MQ=255)
                              cGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATg                                   <  1:1410041/71‑1 (MQ=255)
                               ggTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                                  <  1:2948953/71‑1 (MQ=255)
                               ggTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                                  <  1:2028608/71‑1 (MQ=255)
                                 ttACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACttt                                                            <  1:2215658/43‑1 (MQ=255)
                                 ttACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCgg                                      <  1:1071038/65‑1 (MQ=255)
                                   aCGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTAt                                                          >  1:1203542/1‑43 (MQ=255)
                                    cGGAGGAATCCACACCACCATAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGc                              <  1:2863151/70‑1 (MQ=255)
                                    cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTAtcgt                                                       <  1:147125/45‑1 (MQ=255)
                                    cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGc                              <  1:2459475/70‑1 (MQ=255)
                                       aggaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTa                            >  1:2830392/1‑69 (MQ=255)
                                       aggaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACa                          <  1:284008/71‑1 (MQ=255)
                                        ggaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAg                         >  1:1711045/1‑71 (MQ=255)
                                         gaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAg                         >  1:2944469/1‑70 (MQ=255)
                                           aTCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGc                                                    <  1:883671/41‑1 (MQ=255)
                                           aTCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCt                                                  <  1:1800125/43‑1 (MQ=255)
                                             ccACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTAc                                                >  1:855375/1‑43 (MQ=255)
                                                caccacCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtcgtcg                 <  1:2348687/71‑1 (MQ=255)
                                                caccacCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCAtcgag                  >  1:2460834/1‑68 (MQ=255)
                                                  ccaccaGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGa                >  1:1503923/1‑70 (MQ=255)
                                                   caccaGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGa                >  1:188995/1‑69 (MQ=255)
                                                   caccaGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATa              <  1:349521/71‑1 (MQ=255)
                                                       agagagGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAAttt          <  1:826882/71‑1 (MQ=255)
                                                          gagGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGAt                                    >  1:2193929/1‑42 (MQ=255)
                                                            ggCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGGCGATAATTTTCGCt     <  1:816818/71‑1 (MQ=255)
                                                             gCCGAGAATGACTTTATCGTCGCCTACCTGCTcgct                                       <  1:1944171/36‑2 (MQ=39)
                                                               cGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGa                >  1:558919/1‑57 (MQ=255)
                                                               cGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGc  <  1:1926354/71‑1 (MQ=255)
                                                                gagaATGACTTTATCGTCGCCTACCTGCTCGCGGATg                                   <  1:443309/37‑1 (MQ=255)
                                                                   |                                                                  
TTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGC  >  minE/1664462‑1664595

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: