Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 204,556 A→G 100% intergenic (+156/‑28) alaV → / → rrlH tRNA‑Ala/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE204,5560AG87.1% 133.4 / 7.5 76intergenic (+156/‑28)alaV/rrlHtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (10/0);  new base G (29/37);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCT  >  minE/204503‑204617
                                                     |                                                             
cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACTCGATGATGGATCGCAAGAAAc                                                   >  1:720765/1‑66 (MQ=11)
cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGa                                                      <  1:266285/63‑1 (MQ=255)
cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc                                                >  1:2962856/1‑69 (MQ=255)
cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc                                                >  1:2256100/1‑69 (MQ=255)
          aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                        >  1:496312/1‑67 (MQ=25)
          aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                        >  1:388263/1‑67 (MQ=25)
          aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                        >  1:2397509/1‑67 (MQ=25)
          aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                        >  1:541844/1‑67 (MQ=25)
          aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                        >  1:1215931/1‑67 (MQ=25)
          aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                        >  1:263186/1‑67 (MQ=25)
            aGTTGTTTGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg                                     <  1:1857957/68‑1 (MQ=255)
            aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg                                     <  1:3009114/68‑1 (MQ=25)
            aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg                                     <  1:2310571/68‑1 (MQ=25)
            aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg                                     <  1:415897/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:2736356/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:737290/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:3063870/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:3253540/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:2324132/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:3374924/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:1558341/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:2435209/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:1123118/68‑1 (MQ=25)
                     tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                            <  1:1346909/68‑1 (MQ=25)
                           tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg                     >  1:1059686/1‑69 (MQ=25)
                           tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg                     >  1:1909679/1‑69 (MQ=25)
                           tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg                     >  1:3354438/1‑69 (MQ=25)
                           tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg                     >  1:1038151/1‑69 (MQ=25)
                           tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg                     >  1:3229825/1‑69 (MQ=25)
                           tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg                     >  1:2521328/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc                    >  1:2817025/1‑68 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:3266250/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:3334054/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:2023376/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:1483747/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:605410/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:2716096/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:1055969/1‑69 (MQ=25)
                             tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   >  1:2837402/1‑69 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:2901164/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:388577/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:506467/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:597652/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:1012858/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:2779144/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:1206128/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:1508595/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:1803927/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:1027351/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:198503/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:2029654/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:2092389/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:2109007/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:2323388/67‑1 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                 >  1:3203029/1‑69 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                 >  1:2128705/1‑69 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                 >  1:2955887/1‑69 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                 >  1:1877849/1‑69 (MQ=25)
                               aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                 >  1:1917337/1‑69 (MQ=25)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:3178112/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:114611/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:458026/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:2958314/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:1010527/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:1997439/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:2130271/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:3321832/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:938742/1‑69 (MQ=35)
                                    ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg            >  1:2083933/1‑69 (MQ=35)
                                       gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:947607/68‑1 (MQ=25)
                                       gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:844293/68‑1 (MQ=25)
                                       gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:981/68‑1 (MQ=25)
                                       gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:449054/68‑1 (MQ=25)
                                       gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:2137363/68‑1 (MQ=25)
                                        cAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:2066169/67‑1 (MQ=25)
                                          acacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg          <  1:2843214/65‑1 (MQ=25)
                                                tgatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                   <  1:3290557/50‑1 (MQ=255)
                                                tgatgaATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt  >  1:1373339/1‑67 (MQ=35)
                                                tgatgaATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt  >  1:357248/1‑67 (MQ=35)
                                                     |                                                             
CTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCT  >  minE/204503‑204617

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: