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,5560AG100.0% 185.9 / NA 62intergenic (+156/‑28)alaV/rrlHtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/27);  total (35/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGT  >  minE/204491‑204623
                                                                 |                                                                   
aaaaTTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGc                                                                <  1:208490/71‑1 (MQ=9)
  aaTTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCaa                                                              <  1:2568704/71‑1 (MQ=255)
    ttGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCaa                                                              >  1:3688905/1‑69 (MQ=255)
     tGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGa                                                            <  1:1392035/70‑1 (MQ=255)
      gAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                                           <  1:2323940/70‑1 (MQ=255)
       aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaaa                                                          >  1:2536449/1‑70 (MQ=255)
           acTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCt                                                     <  1:195645/71‑1 (MQ=255)
            cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCtt                                                    >  1:3817665/1‑71 (MQ=255)
               aacaatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCgg                                                 >  1:332373/1‑71 (MQ=255)
               aacaatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCgg                                                 >  1:2611497/1‑71 (MQ=255)
               aacaatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCgg                                                 <  1:2138633/71‑1 (MQ=255)
                     gAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgt                                            >  1:3130730/1‑70 (MQ=25)
                     gAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgt                                            >  1:1858693/1‑70 (MQ=25)
                      aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCaa                                                              >  1:993210/1‑51 (MQ=255)
                      aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCg                                                  <  1:2466782/63‑1 (MQ=25)
                      aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgt                                            >  1:1132867/1‑69 (MQ=25)
                        aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                                        >  1:1861921/1‑71 (MQ=21)
                            gttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGtt                                      >  1:2434665/1‑69 (MQ=25)
                              tCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtg                                             >  1:327633/1‑60 (MQ=255)
                               cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                         >  1:2512618/1‑47 (MQ=255)
                               cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                         >  1:3097743/1‑47 (MQ=255)
                               cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                         >  1:2614614/1‑47 (MQ=255)
                                gTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCg                                 <  1:3572672/70‑1 (MQ=25)
                                gTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCg                                 <  1:3474191/70‑1 (MQ=25)
                                gTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCg                                 <  1:689009/70‑1 (MQ=25)
                                 tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCtt                                                    <  1:2188858/50‑1 (MQ=255)
                                 tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCtt                                                    <  1:1649193/50‑1 (MQ=255)
                                 tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGAc                               <  1:3175320/71‑1 (MQ=21)
                                   aGTCTCTCAAATTTTCGCAACACGATGATGGATCGCGAGAAACATCTTc                                                   >  1:1143203/1‑49 (MQ=11)
                                   aGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCg                                                  >  1:2578747/1‑50 (MQ=255)
                                   aGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTa                             <  1:2148629/71‑1 (MQ=21)
                                    gTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                              >  1:1534779/1‑69 (MQ=25)
                                    gTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                              >  1:1249111/1‑69 (MQ=25)
                                      ctctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGtt                                      >  1:1117239/1‑59 (MQ=255)
                                         tcAAATTTTCGCAACACGATGATGGATCGCAAGaa                                                           <  1:2434833/35‑1 (MQ=255)
                                         tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                              <  1:2364204/64‑1 (MQ=25)
                                         tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                         <  1:3919776/69‑1 (MQ=25)
                                           aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTac                      <  1:2901577/70‑1 (MQ=25)
                                            aaTTTTCGCAACACGATGATGGATCGCAAGAAACATc                                                      <  1:1749396/37‑1 (MQ=255)
                                            aaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGa                                          >  1:3369659/1‑49 (MQ=255)
                                            aaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                         <  1:3566370/66‑1 (MQ=25)
                                              ttttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg                  <  1:1375511/71‑1 (MQ=21)
                                              ttttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg                  >  1:3575004/1‑71 (MQ=21)
                                                ttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg                >  1:2266522/1‑71 (MQ=21)
                                                 tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGt                 <  1:3155149/69‑1 (MQ=25)
                                                     aaCACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGc           <  1:482091/71‑1 (MQ=21)
                                                     aaCACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGc           <  1:1601206/71‑1 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc                          >  1:1558225/1‑54 (MQ=255)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:1091849/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:1672472/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:1177620/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:3219746/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:230552/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:610220/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:614141/1‑71 (MQ=21)
                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc         >  1:2772308/1‑71 (MQ=21)
                                                        acGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAg                                         <  1:1645201/38‑1 (MQ=255)
                                                        acGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt        >  1:278047/1‑71 (MQ=21)
                                                         cGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTg       >  1:2845123/1‑71 (MQ=21)
                                                          gatgatGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                              >  1:107968/1‑47 (MQ=255)
                                                          gatgatGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTgg      <  1:3925653/71‑1 (MQ=21)
                                                           atgatgGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGc     <  1:1564299/71‑1 (MQ=21)
                                                            tgatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGt                 <  1:3535631/58‑1 (MQ=255)
                                                            tgatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTgg      >  1:330614/1‑69 (MQ=25)
                                                             gatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTaca                     >  1:3151695/1‑53 (MQ=255)
                                                              atggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAg   <  1:2941969/70‑1 (MQ=25)
                                                              atggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGt  <  1:2783215/71‑1 (MQ=21)
                                                                 |                                                                   
AAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGT  >  minE/204491‑204623

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: