Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 558,552 G→A 25.9% intergenic (+102/‑202) iaaA → / → cmr L‑asparaginase/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE558,5520GA25.9% 102.7 / 36.7 58intergenic (+102/‑202)iaaA/cmrL‑asparaginase/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (9/6);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAA  >  minE/558483‑558620
                                                                     |                                                                    
aaaGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCgg                                                                     <  1:80741/71‑1 (MQ=255)
 aaGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGc                                                                    <  1:2587882/71‑1 (MQ=255)
 aaGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGc                                                                    <  1:352359/71‑1 (MQ=255)
  aGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAgcagca                                                                   >  1:2054398/1‑71 (MQ=255)
    tGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCAt                                                                  <  1:2312009/70‑1 (MQ=255)
    tGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATc                                                                 <  1:945915/71‑1 (MQ=255)
     gATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCg                                                                >  1:1188140/1‑71 (MQ=255)
      aTTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGTGTAAAGCAGCATCGt                                                               <  1:1830799/71‑1 (MQ=255)
      aTTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGt                                                               <  1:107000/71‑1 (MQ=255)
         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCt                                                             <  1:3091545/70‑1 (MQ=255)
          aaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCt                                                             <  1:1806155/69‑1 (MQ=255)
                    aaTACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTc                                                   >  1:1893242/1‑69 (MQ=255)
                      tACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCg                                                                <  1:1514336/54‑1 (MQ=255)
                      tACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGc                                               >  1:2137045/1‑71 (MQ=255)
                         aGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGc                                            <  1:2309779/71‑1 (MQ=255)
                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGc                                                                    >  1:3566350/1‑45 (MQ=255)
                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGc                                                                    >  1:3339826/1‑45 (MQ=255)
                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTc                                                       <  1:536712/58‑1 (MQ=255)
                            cAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCg                                                                      <  1:2478992/42‑1 (MQ=255)
                            cAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTg                                         >  1:1334649/1‑71 (MQ=255)
                             aaGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCaa                                                 >  1:993844/1‑62 (MQ=255)
                              aGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGtt                                       <  1:2813555/71‑1 (MQ=255)
                              aGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGtt                                       <  1:2431906/71‑1 (MQ=255)
                                tCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTa                                     <  1:1320811/71‑1 (MQ=255)
                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTAttt                                                        <  1:1518486/50‑1 (MQ=255)
                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTAt                                                          >  1:1417727/1‑48 (MQ=255)
                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTAt                                                          >  1:2823331/1‑48 (MQ=255)
                                  gTTGAGAATAAAGTGCACGTTAAACTGGGTAAAGCGGCATCGTCTTAttt                                                        <  1:1227181/50‑1 (MQ=255)
                                   ttGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCgg                                                                     <  1:410514/36‑1 (MQ=255)
                                   ttGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTAttt                                                        <  1:433165/49‑1 (MQ=255)
                                    tGAGAATAATGTGCAGGTTAAACTGGGTAAAGCGGCAt                                                                  <  1:449511/38‑1 (MQ=255)
                                    tGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTc                                                              >  1:739790/1‑42 (MQ=255)
                                    tGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTc                                                              >  1:2017776/1‑42 (MQ=255)
                                     gagaATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggt                                 <  1:1892861/70‑1 (MQ=255)
                                      agaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgg                               <  1:1452881/71‑1 (MQ=255)
                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTccc                                                     >  1:2042055/1‑48 (MQ=255)
                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgg                               <  1:3374784/70‑1 (MQ=255)
                                         aTAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgggt                             <  1:4110457/70‑1 (MQ=255)
                                          tAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGttt                                      >  1:3740527/1‑60 (MQ=255)
                                          tAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTa                                     <  1:3237670/61‑1 (MQ=255)
                                          tAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGa                           >  1:1360687/1‑71 (MQ=255)
                                          tAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGa                           >  1:3522957/1‑71 (MQ=255)
                                             aGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTg                        >  1:2599459/1‑71 (MQ=255)
                                              gTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGt                       <  1:2942778/71‑1 (MQ=255)
                                               tGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTa                      <  1:2845345/71‑1 (MQ=255)
                                                gCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTaa                     >  1:3841079/1‑71 (MQ=255)
                                                 cAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCt                                          >  1:3365632/1‑49 (MQ=255)
                                                 cAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACg                                   >  1:2107775/1‑56 (MQ=255)
                                                    gTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGtt                                       <  1:1796858/49‑1 (MQ=255)
                                                    gTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACgg                                  >  1:1743776/1‑54 (MQ=255)
                                                    gTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACgg                                  >  1:2845926/1‑54 (MQ=255)
                                                      tAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTa                                     >  1:3067251/1‑49 (MQ=255)
                                                      tAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                              >  1:2242253/1‑56 (MQ=255)
                                                      tAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa               >  1:3118572/1‑71 (MQ=255)
                                                      tAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa               >  1:1921087/1‑71 (MQ=255)
                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa               >  1:3467174/1‑71 (MQ=255)
                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa               >  1:935934/1‑71 (MQ=255)
                                                       aaaCTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACgg                                  >  1:2046715/1‑51 (MQ=255)
                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACgg                                  >  1:700536/1‑51 (MQ=255)
                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACgg                                  >  1:837813/1‑51 (MQ=255)
                                                              gTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGttt                                      <  1:2300535/40‑1 (MQ=255)
                                                                    cGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTtaaataaa  <  1:387553/70‑1 (MQ=255)
                                                                     |                                                                    
AAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAA  >  minE/558483‑558620

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: