Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,880,619 C→A 9.3% G271V (GGT→GTT)  rcsC_9 ← Sensor histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,880,6190CA 9.3% 134.7 / 3.5 43G271V (GGT→GTT) rcsC_9Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base A (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGAC  >  ACB143/3880487‑3880735
                                                                                                                                    |                                                                                                                    
aTACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCAt                                                                                                                   >  5:69416/1‑136 (MQ=255)
 tACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATc                                                                                                                  <  8:107789/136‑1 (MQ=255)
         cTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTc                                                                                                          >  7:37353/1‑136 (MQ=255)
           gttgtGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCg                                                                                                         >  5:6786/1‑135 (MQ=255)
                       cgcCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAtctc                                                                                             >  5:106595/1‑135 (MQ=255)
                       cgcCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGAtctc                                                                                             >  3:213983/1‑135 (MQ=255)
                         ccTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAATCCAGAGGTAGAGACCATCGCGGTGTCGCAGTCCGATCTCCat                                                                                          >  4:251693/1‑135 (MQ=255)
                         ccTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGAACATCGCGGTGTCGCAGGCGGATCTCCAg                                                                                          >  4:182093/1‑136 (MQ=255)
                         ccTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGAACATCGCGGTGTCGCAGGCGGATCTCCAg                                                                                          >  8:171135/1‑136 (MQ=255)
                            ctGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAgctg                                                                                       >  1:137732/1‑136 (MQ=255)
                                   aCGCGCAGCCTCTTCACGGTCTGCGGTGATGTNCCGGCCAACCGCCGTCAGGAGGCNGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCcgcg                                                                                <  2:256364/136‑1 (MQ=255)
                                    cgcgCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                               <  8:104011/136‑1 (MQ=255)
                                    cgcgCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                               <  8:276516/136‑1 (MQ=255)
                                    cgcgCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                               <  5:22429/136‑1 (MQ=255)
                                    cgcgCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                               <  8:156234/136‑1 (MQ=255)
                                    cgcgCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGAACATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                               <  3:182093/136‑1 (MQ=255)
                                    cgcgCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGAACATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGt                                                                               <  7:171135/136‑1 (MQ=255)
                                               ttCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAtt                                                                      <  4:146486/134‑1 (MQ=255)
                                               ttCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCAtt                                                                      <  8:182691/134‑1 (MQ=255)
                                                 cACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTgg                                                                    <  4:275448/134‑1 (MQ=255)
                                                  aCGGTCTGCGGTGATGTCCAGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGAGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGTCTGCATTGGAAc                                                                 <  8:54911/136‑1 (MQ=255)
                                                        tGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCg                                                           >  4:111005/1‑136 (MQ=255)
                                                        tGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCg                                                           >  7:34075/1‑136 (MQ=255)
                                                         gCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCg                                                           >  7:25140/1‑135 (MQ=255)
                                                             tGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGacac                                                      <  8:127743/136‑1 (MQ=255)
                                                                  tCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGCCCTTCGAGGTGTGGCAGGGTGATCTCCAGCTNCTCCGNGTGGCNGCATGGGTACAGCGCGGACGCTGc                                                   >  2:231639/1‑134 (MQ=255)
                                                                     cGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATc                                              <  5:98534/136‑1 (MQ=255)
                                                                     cGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATc                                              >  3:297336/1‑136 (MQ=255)
                                                                       gCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCa                                                                                                       <  8:18064/77‑1 (MQ=255)
                                                                       gCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCa                                                                                                       >  7:18064/1‑77 (MQ=255)
                                                                          aaCCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGCCACTGCCAATCGTACt                                         >  1:121556/1‑136 (MQ=255)
                                                                              gccgTCAGAAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATAGTACTTctg                                     >  8:247186/1‑134 (MQ=255)
                                                                               ccgTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCgg                                     >  6:257228/1‑135 (MQ=255)
                                                                               ccgGCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGt                                    <  1:100111/136‑1 (MQ=255)
                                                                                            gCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCggg                       >  7:146053/1‑136 (MQ=255)
                                                                                               gTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATc                                              <  6:57385/110‑1 (MQ=255)
                                                                                               gTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATc                                              >  5:57385/1‑110 (MQ=255)
                                                                                                   aaCCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCTTACTTCGGTACGCTCTGTCGGGTGCACCa                >  3:286453/1‑136 (MQ=255)
                                                                                                   aaCCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCa                <  8:36209/136‑1 (MQ=255)
                                                                                                               cAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCg    <  7:110011/136‑1 (MQ=255)
                                                                                                               cAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCg    <  2:199762/136‑1 (MQ=255)
                                                                                                               cAACGAAACCAGTGGGAGTGACCATCGCGGGGTCGCAGGCGGATCTCGAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCAcc     <  6:25008/135‑1 (MQ=255)
                                                                                                                  cGAAACCAGTGGTAGTGACCATCGCGGAGTCGCAGGCGGATCTCCAGCTGCTCCTCGTGGCAGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGAc  >  7:203203/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                    
ATACGCTCACTGTTGTGCATCAGCGCCTCTGCCTGACGCGCAGCCTCTTCACGGTCTGCGGTGATGTCCCGGCCAACCGCCGTCAGGAGGCTGCCGTCAAACCTCGCACTCCAACGAAACCAGTGGTAGTGACCATCGCGGTGTCGCAGGCGGATCTCCAGCTGCTCCGCGTGGCTGCATTGGAACAGCCCGGACACTGCCAATCGTACTTCGGTACGCTCTGTCGGGTGCACCAGCTCCATCACCGAC  >  ACB143/3880487‑3880735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑