Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,268,310 T→A 14.9% I99F (ATC→TTC)  pgsA ← phosphatidylglycerophosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,268,3100TA14.9% 112.9 / 12.6 47I99F (ATC→TTC) pgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base T (15/25);  new base A (2/5);  total (17/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CACGCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCACCATGGCGA  >  minE/1268242‑1268375
                                                                    |                                                                 
cACGCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGatca                                                                 >  1:3503007/1‑71 (MQ=255)
 aCGCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGatcat                                                                >  1:280320/1‑71 (MQ=255)
 aCGCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCAAGGGCGatcat                                                                >  1:3182565/1‑71 (MQ=255)
  cGCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGataga                                                                >  1:853181/1‑67 (MQ=255)
   gCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCg                                                              <  1:3197070/71‑1 (MQ=255)
   gCTACTGCGTTTACCCAACTCCCCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCg                                                              <  1:601460/71‑1 (MQ=255)
     tACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGtt                                                            >  1:1932873/1‑71 (MQ=255)
       cTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGatc                                                                  >  1:1245592/1‑63 (MQ=255)
       cTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGtt                                                            <  1:3380479/69‑1 (MQ=255)
       cTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTgc                                                          >  1:1612862/1‑71 (MQ=255)
        tGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGatca                                                                 <  1:2682499/63‑1 (MQ=255)
          cGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTgccgc                                                       <  1:2284828/71‑1 (MQ=255)
          cGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTgccg                                                        <  1:317330/70‑1 (MQ=255)
            tttACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCg                                                              <  1:2962703/62‑1 (MQ=255)
            tttACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTgcc                                                         >  1:1943686/1‑67 (MQ=255)
            tttACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGccgccg                                                     >  1:1725269/1‑71 (MQ=255)
            tttACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGAACATCGtt                                                            <  1:3015597/64‑1 (MQ=255)
             ttACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCgg                                                    >  1:2419404/1‑71 (MQ=255)
              tACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCgg                                                    >  1:1211962/1‑70 (MQ=255)
               aCCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGt                                                   <  1:3127376/70‑1 (MQ=255)
               aCCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGt                                                   <  1:1755971/70‑1 (MQ=255)
                cccAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTa                                                  <  1:1993213/70‑1 (MQ=255)
                  cAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGTCGATCATCGTTGccgccg                                                     <  1:526647/65‑1 (MQ=255)
                  cAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGccgccg                                                     <  1:531381/65‑1 (MQ=255)
                   aaCTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGtt                                                            >  1:2352412/1‑57 (MQ=255)
                    aCTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTc                                             <  1:2333344/71‑1 (MQ=255)
                    aCTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGAACATCGTTGCCGCCGGt                                                   <  1:74767/65‑1 (MQ=255)
                     cTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTca                                            <  1:3660674/71‑1 (MQ=255)
                       ccGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGaacatc                                                               >  1:1729028/1‑50 (MQ=255)
                        cGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTcacc                                          >  1:3240817/1‑70 (MQ=255)
                          ccatccatTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTcacccac                                       >  1:3792106/1‑71 (MQ=255)
                           catccatTCGCGTAGCGCAGAAATAATAATTTCACTGGCGATCATCg                                                              <  1:3577692/47‑1 (MQ=255)
                               catTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGc                                   <  1:487229/70‑1 (MQ=255)
                               catTCGCGTAGCGCAGAAATAATAATTTCACGGGCGAACATCGTTGCCGCCGGTAATGTCACCCACCAGCt                                  <  1:1437829/71‑1 (MQ=255)
                                 ttCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGt                                                             <  1:279861/42‑1 (MQ=255)
                                    gcgTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGa                              <  1:412158/70‑1 (MQ=255)
                                     cgTAGCGCAGAAATAATAATTTCACTGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATa                            <  1:873458/71‑1 (MQ=255)
                                      gTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATaa                           >  1:449736/1‑71 (MQ=255)
                                        aGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCtg                                 >  1:2368594/1‑63 (MQ=255)
                                          cgcAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCt                       <  1:2810833/71‑1 (MQ=255)
                                          cgcAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCt                       <  1:1888999/71‑1 (MQ=255)
                                               aaaTAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGtt                  <  1:1338519/71‑1 (MQ=255)
                                                 ataataatTTCACGGGCGATCATCGTTGCCGCCGGTAAt                                                >  1:3727602/1‑39 (MQ=255)
                                                      aatTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCt                                  <  1:3659134/48‑1 (MQ=255)
                                                      aatTTCACGGGCGAACATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCg                     >  1:1056175/1‑61 (MQ=255)
                                                       atTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCtgt                                <  1:1968193/49‑1 (MQ=255)
                                                        tttCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCAc          >  1:1832784/1‑70 (MQ=255)
                                                          tCACGGGCGAACATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCACCa        <  1:1349/70‑1 (MQ=255)
                                                           cACGGGCGAACATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCACCa        <  1:3829173/69‑1 (MQ=255)
                                                            aCGGGCGATCATCGTTGCCGCCGGTAATGTcacccacc                                      <  1:3298170/38‑1 (MQ=255)
                                                             cGGGCGAACATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCACCATGGc    >  1:689738/1‑71 (MQ=255)
                                                              gggCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGtt                  <  1:3816935/56‑1 (MQ=255)
                                                                gCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTAc                >  1:3345442/1‑56 (MQ=255)
                                                                 cGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCACCATGGCGa  >  1:3945438/1‑69 (MQ=255)
                                                                   atcatcGTTGCCGCCGGTAATGTCACCCACCAGCTGTGAt                             <  1:2941984/40‑1 (MQ=255)
                                                                    |                                                                 
CACGCTACTGCGTTTACCCAACTCCGCCATCCATTCGCGTAGCGCAGAAATAATAATTTCACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACCAGCTGTGATAATGCTCGGTTACCAGCACCATGGCGA  >  minE/1268242‑1268375

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: