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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,664,5290GA36.4% 18.2 / 18.8 22I231I (ATC→ATTguaAGMP synthetase
Reads supporting (aligned to +/- strand):  ref base G (9/5);  new base A (5/3);  total (14/8)
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

TTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTC  >  minE/1664462‑1664598
                                                                   |                                                                     
ttttACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGAt                                                                     >  1:131753/1‑70 (MQ=255)
     ccGATAGCGCGGTTCAGCAGCATTGCTGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACtt                                                                <  1:1472523/70‑1 (MQ=255)
        aTAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTAtc                                                            >  1:663465/1‑71 (MQ=255)
            cgcgGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTAtcgtcg                                                        >  1:670766/1‑71 (MQ=255)
               ggTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCt                                                     <  1:1787550/71‑1 (MQ=255)
                     gcagcaTTGAGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTgg                                               <  1:1010212/71‑2 (MQ=255)
                          aTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTcgcg                                          <  1:2087091/71‑1 (MQ=255)
                           ttGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTc                                             >  1:179265/1‑67 (MQ=255)
                           ttGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACtt                                                                >  1:607141/1‑48 (MQ=255)
                           ttGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACtt                                                                >  1:1761090/1‑48 (MQ=255)
                                gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCg                                    <  1:992635/71‑1 (MQ=255)
                                 ttACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGa                                   <  1:2034712/71‑1 (MQ=255)
                                  tACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAg                                  >  1:1292789/1‑71 (MQ=255)
                                   aCGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGcc                                                      >  1:607164/1‑50 (MQ=255)
                                    cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGc                                 <  1:1627935/70‑1 (MQ=255)
                                       aggaATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTcgcg                                          <  1:1094980/58‑1 (MQ=255)
                                           aTCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCa                          >  1:2102618/1‑70 (MQ=255)
                                              cacaCCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTAc                                                   >  1:2103405/1‑42 (MQ=255)
                                                  ccaccaGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAg                            >  1:2181140/1‑61 (MQ=255)
                                                               cGAGGATGACTTTATCGTCGCCTACCTGCTCGCgg                                         >  1:1324813/1‑35 (MQ=255)
                                                               cGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGc     >  1:1280348/1‑71 (MQ=255)
                                                               cGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAAtttt            >  1:1581772/1‑64 (MQ=255)
                                                                  ggATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTc  <  1:1340555/71‑1 (MQ=255)
                                                                   |                                                                     
TTTTACCGATAGCGCGGTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGCTGGCGTC  >  minE/1664462‑1664598

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