Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,020,631 Δ1 bp 100% coding (222/927 nt) yqhG → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,020,6310T.86.1% 115.8 / 16.1 36coding (222/927 nt)yqhGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base . (0/31);  total (0/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGG  >  minE/2020590‑2020675
                                         |                                            
gcgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGttt                     <  1:1158543/67‑1 (MQ=255)
gcgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGttt                     <  1:414784/67‑1 (MQ=255)
gcgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGttt                     <  1:383912/67‑1 (MQ=255)
gcgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGttt                     <  1:222152/67‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTAGGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:779235/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1441637/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:770007/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:764616/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:711989/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:626396/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:526567/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:522358/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:492773/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:917918/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1078772/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:929725/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:936344/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:2245375/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1391994/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:2167292/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:2081712/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1518341/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1547100/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1658888/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1709031/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1803454/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1851795/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1887502/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:1943354/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:19676/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:970002/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:2002269/63‑1 (MQ=255)
 cgcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:212052/63‑1 (MQ=255)
  gcTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:27237/62‑1 (MQ=255)
       cccTGAACGACAATCTTACGCTGCCGGGCTATAT‑‑‑‑TGGCTGGATTATGATCCGGTTtg                    <  1:83949/57‑1 (MQ=255)
                   aTCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:1762915/67‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:1420944/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:971661/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:899617/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:819391/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:456318/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:118336/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:1992954/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  <  1:1931550/46‑1 (MQ=255)
                                        ttGGCTGGCTGGATTATGATCCGGATTGTGATTGTCAGGATTTTgg  <  1:2305174/46‑1 (MQ=255)
                                         |                                            
GCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGG  >  minE/2020590‑2020675

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: