Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,281,038 G→A Q95* (CAG→TAG)  fabR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,281,0380GA97.7% 133.6 / ‑3.8 43Q95* (CAG→TAG) fabRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/19);  total (24/19)
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.90e-01

CTCCATAAATGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACG  >  minE/2280973‑2281102
                                                                 |                                                                
cTCCATAAATGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCa                                                             >  1:98329/1‑71 (MQ=255)
     tAAATGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGt                                                        <  1:951996/71‑1 (MQ=255)
      aaaTGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGtt                                                       >  1:2151179/1‑71 (MQ=255)
             gagaCCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAgcg                                                               <  1:76665/56‑1 (MQ=255)
             gagaCCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAg                                                         <  1:1812043/62‑1 (MQ=255)
                aCCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCAt                                             >  1:1982879/1‑71 (MQ=255)
                  cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCAtt                                            >  1:124953/1‑70 (MQ=255)
                  cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTa                                           <  1:1241694/71‑1 (MQ=255)
                   gAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCAt                                                            >  1:1042288/1‑53 (MQ=255)
                        gCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCg                                      >  1:2178749/1‑70 (MQ=255)
                         cGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTg                                                      >  1:1464362/1‑53 (MQ=255)
                         cGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTg                                                      >  1:1816666/1‑53 (MQ=255)
                         cGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGc                                     >  1:626037/1‑70 (MQ=255)
                           gATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTa                                           <  1:71444/62‑1 (MQ=255)
                           gATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTa                                           >  1:1439564/1‑62 (MQ=255)
                           gATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGctct                                  <  1:889707/71‑1 (MQ=255)
                            aTCACACTCCCGCCTTTAGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGctct                                  >  1:1853386/1‑70 (MQ=255)
                              cacaCTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTg                                                      >  1:2323714/1‑48 (MQ=255)
                               acacTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAg                                                         <  1:1979901/44‑1 (MQ=255)
                                  cTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGt                                                        <  1:51390/42‑1 (MQ=255)
                                  cTCCCGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAAc                                        <  1:323953/58‑1 (MQ=255)
                                     ccGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCtcatc                              <  1:703382/65‑1 (MQ=255)
                                     ccGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAAcc                          <  1:1230826/69‑1 (MQ=255)
                                     ccGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCAt                        >  1:90231/1‑71 (MQ=255)
                                      cGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATg                       >  1:468005/1‑71 (MQ=255)
                                      cGCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATg                       <  1:1985264/71‑1 (MQ=255)
                                       gCCTTTGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATgg                      <  1:1638883/71‑1 (MQ=255)
                                            tGGCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCAt                                             >  1:1476063/1‑43 (MQ=255)
                                             ggCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTaa                                          >  1:2170655/1‑45 (MQ=255)
                                             ggCGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTaa                                          >  1:2086893/1‑45 (MQ=255)
                                               cGATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGa                 <  1:159923/68‑1 (MQ=255)
                                                gATACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                   <  1:287479/65‑1 (MQ=255)
                                                  tACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGt           >  1:2204731/1‑71 (MQ=255)
                                                  tACGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGt           <  1:1097674/71‑1 (MQ=255)
                                                   aCGCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGt           >  1:2147401/1‑70 (MQ=255)
                                                     gCTGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCg        <  1:932504/71‑1 (MQ=255)
                                                       tGACGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                   >  1:1536221/1‑58 (MQ=255)
                                                         aCGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTc      <  1:568547/69‑1 (MQ=255)
                                                         aCGCGCCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa    >  1:102560/1‑71 (MQ=255)
                                                           gcgcCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTc      <  1:649243/67‑1 (MQ=255)
                                                           gcgcCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACg  >  1:1131583/1‑71 (MQ=255)
                                                             gcCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCtcat                               >  1:1615222/1‑40 (MQ=255)
                                                             gcCTAGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCt     >  1:2100232/1‑66 (MQ=255)
                                                                 |                                                                
CTCCATAAATGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACG  >  minE/2280973‑2281102

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: