Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,068,677 A→G 5.3% N371N (AAT→AAC aer ← fused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,068,6770AG 5.3% 221.2 / 2.7 75N371N (AAT→AACaerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base A (45/26);  new base G (2/2);  total (47/28)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

ACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAAT  >  minE/2068608‑2068739
                                                                     |                                                              
aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCa                                                                <  1:60279/70‑1 (MQ=255)
  gCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTc                                                             <  1:5172089/71‑1 (MQ=255)
  gCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTc                                                             <  1:1810667/71‑1 (MQ=255)
          ccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAg                                                     <  1:2222441/71‑1 (MQ=255)
          ccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAg                                                     <  1:2172191/71‑1 (MQ=255)
              ccaccaCTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGaa                                                   <  1:5181832/69‑1 (MQ=255)
                     tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGaa                                                   <  1:835635/62‑1 (MQ=255)
                     tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGaa                                                   <  1:143719/62‑1 (MQ=255)
                     tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGaa                                                   <  1:2484975/62‑1 (MQ=255)
                      gCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                                                 >  1:3676150/1‑63 (MQ=255)
                      gCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                                                 >  1:4451688/1‑63 (MQ=255)
                      gCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                                                 >  1:2576390/1‑63 (MQ=255)
                       cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                                                 >  1:1703148/1‑62 (MQ=255)
                       cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTg                                         <  1:2549873/70‑1 (MQ=255)
                              tttGCCCTGTTCGCCGGCACGCGCCGCTTCCACGGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCaa                                 >  1:273898/1‑71 (MQ=255)
                              tttGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCa                                  >  1:403341/1‑70 (MQ=255)
                               ttGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCaa                                 >  1:1781715/1‑70 (MQ=255)
                                     tGTTCGCCGGCGCGCGCCGCTTCCACCGCCGCGTTCAgggg                                                        >  1:5038773/1‑40 (MQ=37)
                                     tGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGt                           >  1:3682060/1‑70 (MQ=255)
                                       ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTaa                        >  1:289922/1‑71 (MQ=255)
                                       ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTaa                        >  1:3922155/1‑71 (MQ=255)
                                       ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTaa                        >  1:2614643/1‑71 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:191333/1‑46 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:719737/1‑46 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:3076514/1‑46 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:299026/1‑46 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:2596079/1‑46 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:5187847/1‑46 (MQ=255)
                                        tCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                                                >  1:512322/1‑46 (MQ=255)
                                         cGCCCGCGCGCGCCGCTTCCACGGCGGCGTTCAGGGCCAGGATATTGGTCTGGAAGGCGatgcc                             >  1:1993573/1‑61 (MQ=255)
                                          gCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAatat                              <  1:2639550/62‑1 (MQ=255)
                                          gCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTCGTCTGAAACGCAATATCGTTAATCa                     <  1:4862538/71‑1 (MQ=255)
                                           ccGGCGCGCGCCGCTTCCACCGCCGCGTTCAGCGCCAGGATATTGGTCTGGAAGGc                                   <  1:3308771/56‑1 (MQ=25)
                                           ccGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    <  1:3550580/71‑1 (MQ=255)
                                              gCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg                  >  1:3789691/1‑70 (MQ=255)
                                              gCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg                  >  1:5128810/1‑70 (MQ=255)
                                               cACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGaa                <  1:1715739/71‑1 (MQ=255)
                                               cACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGaa                <  1:1834508/71‑1 (MQ=255)
                                               cACGCGCCGCTTCCAACGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGaa                <  1:4237223/71‑1 (MQ=255)
                                                  gcgcCGCTTCGACCGCCGCGTTCAGCGCCAGGATATTGGTCTGGAAGGCGATGCCGTCGATCACCGagat              <  1:4356027/70‑4 (MQ=255)
                                                  gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTa             <  1:2718333/71‑1 (MQ=255)
                                                  gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTa             <  1:5057359/71‑1 (MQ=255)
                                                   cgccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTaa            <  1:4932045/71‑1 (MQ=255)
                                                    gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTaa            <  1:697639/70‑1 (MQ=255)
                                                    gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTaa            >  1:4989854/1‑70 (MQ=255)
                                                    gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAAt           >  1:885271/1‑71 (MQ=255)
                                                    gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAAt           >  1:3446729/1‑71 (MQ=255)
                                                    gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAAt           >  1:2278701/1‑71 (MQ=255)
                                                    gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAAt           >  1:2466648/1‑71 (MQ=255)
                                                     ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg          <  1:1426218/71‑1 (MQ=255)
                                                     ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg          <  1:601432/71‑1 (MQ=255)
                                                     ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg          <  1:4930800/71‑1 (MQ=255)
                                                     ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg          <  1:2987709/71‑1 (MQ=255)
                                                     ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg          <  1:1844499/71‑1 (MQ=255)
                                                     ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg          <  1:5015128/71‑1 (MQ=255)
                                                           ccACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGcc     <  1:79837/70‑1 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAAt                       >  1:3719988/1‑50 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATc                      >  1:1051065/1‑51 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATc                      >  1:15331/1‑51 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATc                      >  1:2244833/1‑51 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATc                      >  1:5083953/1‑51 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATc                      >  1:4885081/1‑51 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAAt  >  1:2473425/1‑71 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAAt  >  1:5133052/1‑71 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAAt  >  1:713830/1‑71 (MQ=255)
                                                             aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAAt  >  1:4538841/1‑71 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCg                            >  1:4408153/1‑42 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:2503798/1‑50 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:2867355/1‑50 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:5197193/1‑50 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:1553461/1‑50 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:4664456/1‑50 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:797935/1‑50 (MQ=255)
                                                                gccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg                    >  1:1131523/1‑50 (MQ=255)
                                                                    cATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCaa   >  1:1696594/1‑63 (MQ=255)
                                                                    cATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCaa   >  1:3384665/1‑63 (MQ=255)
                                                                     |                                                              
ACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAAT  >  minE/2068608‑2068739

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: