Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 826,990 G→A 14.9% N173N (AAC→AAT narL ← DNA‑binding response regulator in two‑component regulatory system with NarX (or NarQ)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE826,9900GA14.9% 178.7 / 20.8 67N173N (AAC→AATnarLDNA‑binding response regulator in two‑component regulatory system with NarX (or NarQ)
Reads supporting (aligned to +/- strand):  ref base G (38/19);  new base A (3/7);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 3.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TCTTCAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACATCGCG  >  minE/826920‑827058
                                                                      |                                                                    
tcttcAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTg                                                                      >  1:3958990/1‑71 (MQ=255)
 cttcAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTg                                                                      <  1:177533/70‑1 (MQ=255)
 cttcAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTg                                                                      <  1:1445105/70‑1 (MQ=255)
    cAGCATGTGCTTTACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGt                                                                     <  1:111606/68‑1 (MQ=255)
    cAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGt                                                                     <  1:1147690/68‑1 (MQ=255)
    cAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGt                                                                     <  1:414506/68‑1 (MQ=255)
    cAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCg                                                                  >  1:3375510/1‑71 (MQ=255)
     aGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCgg                                                                 >  1:1021884/1‑71 (MQ=255)
           tgCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAAcc                                                           >  1:146018/1‑71 (MQ=255)
                  cGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCa                                                    <  1:2031980/71‑1 (MQ=255)
                      cACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAAt                                                  >  1:2545881/1‑69 (MQ=255)
                       aCTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAg                                               >  1:1215267/1‑71 (MQ=255)
                        cTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCaa                                                   <  1:2462129/66‑1 (MQ=255)
                           ttACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCt                                                         >  1:2992064/1‑57 (MQ=255)
                              cTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTgag                                         <  1:3644285/70‑1 (MQ=255)
                               tgtgCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAAccc                                                          <  1:2149272/52‑1 (MQ=255)
                               tgtgCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAAccc                                                          <  1:3029747/52‑1 (MQ=255)
                                gtgCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGc                                                                >  1:1809819/1‑45 (MQ=255)
                                 tgCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:3982976/1‑71 (MQ=255)
                                   cTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATc                                   >  1:2925188/1‑71 (MQ=255)
                                      tCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATcgc                                 <  1:3040549/70‑1 (MQ=255)
                                      tCGGTGATATCCAGGCGGCGGGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGATAATATcgcg                                >  1:139631/1‑71 (MQ=255)
                                       cGGTGATATCCAGGCGGCGGGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTgag                                         >  1:318033/1‑61 (MQ=255)
                                       cGGTGATATCCAGGCGGCGGGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGAGAATATcgcg                                <  1:2128160/70‑1 (MQ=255)
                                        ggTGATATCCAGGCGGCGGGCAATCATCTTGTTCgg                                                                 >  1:3292379/1‑36 (MQ=255)
                                        ggTGATATCCAGGCGGCGGGCAATCATCTTATTCGGCAAACCCTGCGCAATCAGCTTGAGAATATCgcgc                               <  1:3425557/70‑1 (MQ=255)
                                         gTGATATCCAGGCGGCGGGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGc                                              <  1:1777104/54‑1 (MQ=255)
                                             tatCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGaa                                       >  1:818911/1‑57 (MQ=255)
                                             tatCCAGGCGGCGGGCAATCATCTTGTTAGGCAAACCCTGGGCAATCAGCTTGAGaa                                       >  1:172436/1‑57 (MQ=255)
                                               tCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGaa                                       <  1:4113244/55‑1 (MQ=255)
                                               tCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTgg                       <  1:47960/71‑1 (MQ=255)
                                                 cAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTg                                                        >  1:1038910/1‑36 (MQ=255)
                                                 cAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:1283258/1‑55 (MQ=255)
                                                 cAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTgggg                     >  1:4059526/1‑71 (MQ=255)
                                                 cAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTgggg                     >  1:2553406/1‑71 (MQ=255)
                                                    gcggcgGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGaa                                       >  1:2409763/1‑50 (MQ=255)
                                                    gcggcgGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGtt                   >  1:3946863/1‑70 (MQ=255)
                                                    gcggcgGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTa                  <  1:1919354/71‑1 (MQ=255)
                                                    gcggcgGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTa                  <  1:2890090/71‑1 (MQ=255)
                                                    gcggcgGGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTa                  <  1:2134456/71‑1 (MQ=255)
                                                     cggcggGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:2452850/1‑71 (MQ=255)
                                                       gcggGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTgggg                     <  1:2428141/65‑1 (MQ=255)
                                                       gcggGCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTgggg                     <  1:2764155/65‑1 (MQ=255)
                                                         gggCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:1497627/1‑67 (MQ=255)
                                                          ggCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGt            >  1:2982234/1‑71 (MQ=255)
                                                          ggCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGt            <  1:2144187/71‑1 (MQ=255)
                                                          ggCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGt            >  1:1885130/1‑71 (MQ=255)
                                                          ggCAATCATCTTATTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTgg                       <  1:1281114/60‑1 (MQ=255)
                                                            cAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGtt           >  1:1920157/1‑70 (MQ=255)
                                                             aaTCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTaa         <  1:3274092/71‑1 (MQ=255)
                                                              atcatcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAAc        <  1:3632495/71‑1 (MQ=255)
                                                               tcatcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTc                             <  1:2660929/49‑1 (MQ=255)
                                                               tcatcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:1748097/1‑61 (MQ=255)
                                                               tcatcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:1083874/1‑61 (MQ=255)
                                                               tcatcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:3112414/1‑61 (MQ=255)
                                                                catcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACAt      >  1:3893112/1‑71 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGaa                                       >  1:1314338/1‑37 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:2932235/1‑39 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:2450057/1‑39 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:354609/1‑39 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:2146109/1‑59 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTaa                 >  1:1497749/1‑59 (MQ=255)
                                                                 atcTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGGTAACTGGTTAACATc     <  1:3306192/71‑1 (MQ=255)
                                                                   cTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATc                                   <  1:3603459/39‑1 (MQ=255)
                                                                    ttGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:3102834/1‑36 (MQ=255)
                                                                    ttGTTCGGCAAACCCTGGGCAATCAGCTTGAGAata                                     >  1:1462625/1‑36 (MQ=255)
                                                                    ttGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATc                                   >  1:2211201/1‑38 (MQ=255)
                                                                    ttGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACATcgcg  >  1:2626996/1‑71 (MQ=255)
                                                                    ttGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACATcgcg  >  1:1191117/1‑71 (MQ=255)
                                                                    ttGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACATcgcg  >  1:3292245/1‑71 (MQ=255)
                                                                      gTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACa       <  1:2190546/64‑1 (MQ=255)
                                                                      |                                                                    
TCTTCAGCATGTGCTTGACGTGCACTTTTACTGTGCTTTCGGTGATATCCAGGCGGCGGGCAATCATCTTGTTCGGCAAACCCTGGGCAATCAGCTTGAGAATATCGCGCTCGCGTGGGGTTAACTGGTTAACATCGCG  >  minE/826920‑827058

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: