Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,065,799 T→G 67.8% F150V (TTT→GTT)  rcsC_6 → Sensor histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,065,7990TG67.8% 69.8 / 34.1 52F150V (TTT→GTT) rcsC_6Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base G (21/14);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATGCCGCCGAAGCCGCCAACCGCA  >  ACB122/3065669‑3065934
                                                                                                                                  |                                                                                                                                       
tACCGCGGGGCGCGCCCGTGAGCACTGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGGGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGt                                                                                                                                     <  1:173925/135‑1 (MQ=255)
  ccGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGc                                                                                                                                  >  4:101914/1‑136 (MQ=255)
        tGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACAt                                                                                                                            >  4:384167/1‑136 (MQ=255)
             gCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATGTCGATGTGCTGATGACCCTGCTGCTCAAACACGGCCATGAAGGGCTGGTGGAGCGGTTTGTCGCCGACATTTCCg                                                                                                                       >  4:83540/1‑136 (MQ=255)
             gCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCg                                                                                                                       >  7:152663/1‑136 (MQ=255)
                   gAGCACAGCCTGCACGGCTACGAAACCCGTCGGCGGCGCAAGGACGGCAGCCATATCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCa                                                                                                                 <  4:113252/136‑1 (MQ=255)
                    aGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCa                                                                                                                 <  3:361825/135‑1 (MQ=255)
                    aGCACAGCCTGCACGGCTACGAAACCCGGCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAACCCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCa                                                                                                                 <  2:98729/135‑1 (MQ=255)
                     gCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCa                                                                                                                 <  2:375702/134‑1 (MQ=255)
                     gCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCAGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCaa                                                                                                                >  5:333967/1‑135 (MQ=255)
                           ccTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGcc                                                                                                         <  4:209344/136‑1 (MQ=255)
                            cTGCACGGCTACTAAACCCGGCTGCGGCGCAAGGACGGCAGCCAGCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGcc                                                                                                         <  2:217291/135‑1 (MQ=255)
                            cTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGc                                                                                                          >  5:155626/1‑134 (MQ=255)
                            cTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGCCATTACCGAACGCAAGCTGGc                                                                                                          >  2:230122/1‑134 (MQ=255)
                            cTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCAGCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCAGGAAGGGCTGCTGGAGCGGTTTGTCGCCGCCAGTACCGAACGCAAGCTGGca                                                                                                         >  3:263634/1‑134 (MQ=255)
                            cTGCACGGCTACGAAACCCGTCTGCGGCGCAAGCACGGCAGCCATCTCGATGTGCTGATGAACCTGCAGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTCGccc                                                                                                        >  7:361464/1‑136 (MQ=255)
                              gCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGCCCTGGAGGAGGGGGTTGTCTCCGACATTACCGAACTCCAGCTCGCCca                                                                                                       >  2:177639/1‑135 (MQ=255)
                               cACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATCAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCca                                                                                                       <  4:339224/134‑1 (MQ=255)
                                actgCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTTCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCca                                                                                                       <  1:15841/130‑1 (MQ=255)
                                aCGGCTACGAAACCCGGCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCca                                                                                                       <  1:290603/133‑1 (MQ=255)
                                aCGGCTACCAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGACCTGATCAACCTGCAGCTCAAACCCGGCCATGACGGGCTGGGGGAGGGGGATGTCGCCGACATTACCGAACGCAAGCTGGCCca                                                                                                       <  5:64108/133‑1 (MQ=255)
                                 cGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCcagcag                                                                                                   >  8:340907/1‑136 (MQ=255)
                                      aCGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGcc                                                                                               >  6:84916/1‑135 (MQ=255)
                                      aCGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCt                                                                                              >  3:321552/1‑136 (MQ=255)
                                      aCGAAACCAGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCt                                                                                              >  4:59433/1‑136 (MQ=255)
                                       cGACACCCGGCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTg                                                                                             <  8:94816/136‑1 (MQ=255)
                                        gAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTg                                                                                             <  1:214083/135‑1 (MQ=255)
                                        caaaCCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGcc                                                                                               <  1:394300/132‑1 (MQ=255)
                                         caaCCCGACTGCGGCGCACGGCCGGCAGCCAGCTCGATGTGCTGATGAACCTGCTGCGCACACCCGGCCAGGACGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGc                                                                                            <  1:198603/134‑1 (MQ=255)
                                          aaCCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGGCTGca                                                                                           >  1:200920/1‑135 (MQ=255)
                                          aaCCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGcaa                                                                                          >  8:61618/1‑136 (MQ=255)
                                                 cTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAACCCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCca                                                                                                       <  7:296402/116‑1 (MQ=255)
                                                 cTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACCACTc                                                                                    >  3:74228/1‑135 (MQ=255)
                                                 cTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCca                                                                                                       >  8:296402/1‑116 (MQ=255)
                                                     ggcgCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATg                                                                                >  5:236220/1‑135 (MQ=255)
                                                       cgcAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGAc                                                                              <  6:139130/135‑1 (MQ=255)
                                                            ggACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGCCGAACTg                                                                        >  8:169125/1‑136 (MQ=255)
                                                              aCGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGa                                                                      <  4:130002/136‑1 (MQ=255)
                                                                   aGCCATCTCGATGTGATGATGAACCTGCTGCTCAAACCCGGCCAAGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGAAAGCTGGCCCAGCAGCGCCTGCAACACCTCAATGACGAACTGCAGCAAc                                                                 >  7:200195/1‑136 (MQ=255)
                                                                                    gatgaACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTAcc                                                >  7:68379/1‑136 (MQ=255)
                                                                                     atgaACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTAc                                                 >  6:241469/1‑134 (MQ=255)
                                                                                     atgaACCTGCTGCTCAAACCCAGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTAc                                                 >  5:160553/1‑134 (MQ=255)
                                                                                      tgaACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGa                                              <  8:344479/136‑1 (MQ=255)
                                                                                       gaACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGAt                                             >  8:257143/1‑136 (MQ=255)
                                                                                               tgctCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCCGCAGCGCCTGCACCAACTCAAAGACGAACAGGAGCAACGGGATGCTCCGCGTAACGATGAATTGc                                      >  5:355884/1‑135 (MQ=255)
                                                                                                gctCAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCt                                     >  1:193654/1‑135 (MQ=255)
                                                                                                 cgcAAACCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTg                                    <  3:101914/133‑1 (MQ=255)
                                                                                                      aCCCGGCCATGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGgcg                              >  8:72810/1‑136 (MQ=255)
                                                                                                              aTGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATgc                      >  3:57081/1‑136 (MQ=255)
                                                                                                              aTGAAGGGCTGGTGGAGGGGGTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATg                       >  1:150282/1‑135 (MQ=255)
                                                                                                                        ggtggAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATGCCGCCGAAgcc            <  3:59433/136‑1 (MQ=255)
                                                                                                                         gtggAGGGGGTTGTCGCAGACATTACCGAACTCAAGCTCGCCCAGCAGCGCCTGCAAGAACTCAATGACGAAGTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATGCCGCCGAAggc            <  1:177639/135‑3 (MQ=255)
                                                                                                                                  tttGTCGCCGACATTACCGAACGCAAGCTGGCCCATCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATGCCGCCGAAGCCGCCAACCgca  <  7:178820/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
TACCGCGGTGCTCGCCCGTGAGCACAGCCTGCACGGCTACGAAACCCGTCTGCGGCGCAAGGACGGCAGCCATCTCGATGTGCTGATGAACCTGCTGCTCAAACCCGGCCATGAAGGGCTGGTGGAGGGGTTTGTCGCCGACATTACCGAACGCAAGCTGGCCCAGCAGCGCCTGCAACAACTCAATGACGAACTGGAGCAACGGGTTGCTGCGCGTACCGATGAATTGCTGGAGGCGCGTGATGCCGCCGAAGCCGCCAACCGCA  >  ACB122/3065669‑3065934

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