Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,351,511 A→G 36.1% V457A (GTG→GCG)  arcB ← hybrid sensory histidine kinase in two‑component regulatory system with ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,351,5110AG36.1% 29.3 / 34.6 36V457A (GTG→GCG) arcBhybrid sensory histidine kinase in two‑component regulatory system with ArcA
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base G (9/4);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCA  >  W3110S.gb/3351442‑3351577
                                                                     |                                                                   
gCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCAc                                                                    <  1:1409895/71‑1 (MQ=255)
      cgacgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCGCCTGGt                                                               <  1:1085867/70‑1 (MQ=255)
      cgacgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGt                                                               >  1:2207580/1‑70 (MQ=255)
       gacgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTaa                                                             <  1:2337893/71‑1 (MQ=255)
       gacgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTaa                                                             <  1:1479729/71‑1 (MQ=255)
         cgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCGCCTGGTAATa                                                           <  1:1946160/71‑1 (MQ=255)
         cgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATa                                                           <  1:846788/71‑1 (MQ=255)
          gagaAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCGCCTGGTAAt                                                            >  1:880333/1‑69 (MQ=255)
          gagaAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCGCCTGGTAATa                                                           >  1:176401/1‑70 (MQ=255)
                cGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGc                                                     >  1:2472592/1‑70 (MQ=255)
                        ccAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAAtttt                                           <  1:36214/71‑1 (MQ=255)
                        ccAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAAttt                                            >  1:22745/1‑70 (MQ=255)
                              ccggtgccggtgGCAGGTTTACCGCCATGACTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAg                                     >  1:2792102/1‑71 (MQ=255)
                                    ccggtgGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCAt                                <  1:1445508/70‑1 (MQ=255)
                                       gtgGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTAt                                         <  1:78277/58‑1 (MQ=255)
                                       gtgGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTg                            >  1:2039610/1‑71 (MQ=255)
                                          gCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATa                                                           <  1:2791654/38‑1 (MQ=255)
                                            aGGTTTACCGCCATGACTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAGCTCATc                               >  1:2212133/1‑63 (MQ=255)
                                               tttACCGCCATGACTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGc                    >  1:1100091/1‑71 (MQ=255)
                                                ttACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAAATTTTATCCAGCTCATCCTGCGGAATGc                    >  1:395170/1‑71 (MQ=255)
                                                ttACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCa                                 <  1:1205079/57‑1 (MQ=255)
                                                   ccGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGa                 <  1:348412/70‑1 (MQ=255)
                                                         tGACTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAg                                     >  1:1081445/1‑44 (MQ=255)
                                                           aCTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGa                 >  1:2563620/1‑62 (MQ=255)
                                                           aCTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATc                                        >  1:2033148/1‑39 (MQ=255)
                                                           aCTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATc                                        >  1:1934707/1‑39 (MQ=255)
                                                           aCTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCagag         <  1:1461329/70‑1 (MQ=255)
                                                            cTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATAc              <  1:1813981/64‑1 (MQ=255)
                                                            cTGTCTTTCGCCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATAc              <  1:380113/64‑1 (MQ=255)
                                                            cTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTc       <  1:965867/71‑1 (MQ=255)
                                                               tCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTc       >  1:126128/1‑68 (MQ=255)
                                                               tCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTc    >  1:630564/1‑71 (MQ=255)
                                                               tCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTc    >  1:2095434/1‑71 (MQ=255)
                                                                 tttCGCCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCa  >  1:2509759/1‑71 (MQ=255)
                                                                 tttCGCCTGGTAATACACGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCa  >  1:1998191/1‑71 (MQ=255)
                                                                 tttCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGcc                   <  1:1222490/54‑1 (MQ=255)
                                                                     |                                                                   
GCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCA  >  W3110S.gb/3351442‑3351577

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: