Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,029,463 A→C 31.3% D290A (GAC→GCC)  mcpU_3 → Methyl‑accepting chemotaxis protein McpU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,029,4630AC31.3% 41.1 / 2.7 43D290A (GAC→GCC) mcpU_3Methyl‑accepting chemotaxis protein McpU
Reads supporting (aligned to +/- strand):  ref base A (14/14);  new base C (9/5);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGTGCTGCTG  >  ACB122/4029331‑4029596
                                                                                                                                    |                                                                                                                                     
aGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGCTGCCCCCCAGCTGGCCCAGCGCGTCGGCCCCGGAGAccc                                                                                                                                    >  7:85193/1‑136 (MQ=255)
        aGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATTCCAGCAAGCTGGCCCCGCGCTTCGCCCAGGTAGCCCCGAACAAAg                                                                                                                            >  4:110844/1‑136 (MQ=255)
        aGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCAGGCCCCGCGCTTCGACCAGGTAGACCCGAACAAAg                                                                                                                            >  2:20110/1‑136 (MQ=255)
          ccAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCATCGACCAGGTAGACCCGGCCAAACgc                                                                                                                          >  1:63460/1‑136 (MQ=255)
          ccAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACAACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCACCCACCTTGCCCAGCGCTTAGACCAGGTAGCCCCGCCCAACGgc                                                                                                                          >  7:110765/1‑136 (MQ=255)
                 gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCCAGCAGGCCCCCCGCGTCGACCAGGTAGACCCGACCAAACGCCCCgcgc                                                                                                                   >  4:125057/1‑136 (MQ=255)
                     gCCGCAGCCTCTACGCAGGCCGCACCACGGTCGGCATTCTCACCCCGGGCGGGCGGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGt                                                                                                               <  7:182451/136‑1 (MQ=255)
                     gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCTGGCCCCGCGCTTCGCCCAGGGAGCCCCGCCCAAAGGCGCCGAGCGGGt                                                                                                               >  2:41189/1‑136 (MQ=255)
                     gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCTGGCCCAGCGCGTCGACCAGGGAGCCCCGCCCAAAGGCGCCGAGCTGGt                                                                                                               >  7:64003/1‑136 (MQ=255)
                     gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTAACCCCGGTCGGTCTGTTGGCTGGCTACATCGCCGATTGCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAACGCGCCGCggtggt                                                                                                               >  2:32557/1‑136 (MQ=255)
                          aGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCa                                                                                                          <  7:55561/136‑1 (MQ=255)
                                         cGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCTGGCCCAGCGCTTCGACCAGGTAGCCCCGCCCAAAGGCGCCGCGCTGGTGCGCAAGCTGGCCGATGGc                                                                                            >  6:200358/1‑135 (MQ=255)
                                          gCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGCTGCCAGCAACCTGGCCCCCCGGATGGCCCAGGGAGCCCCGCCCAAAGGCCCCGAGCGGGTGCGCAACCTGGCCGGTGGc                                                                                            >  3:184244/1‑134 (MQ=255)
                                          gCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCaa                                                                                          >  3:120962/1‑136 (MQ=255)
                                           caccacGGTCGGCATTCTCCCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATTCCAGCACGCTGGCCCCGCGCTTCGCCCCGGTAGGCCCGCCCAAAGGCGCCGAGCGGGTGCGCCAGCTGGCCGATGGCacg                                                                                         >  6:35000/1‑134 (MQ=255)
                                           caccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAg                                                                                         >  4:103453/1‑136 (MQ=255)
                                           caccacGGTCGGCATTCTCACCCCGGTCAGTCTGTTGGCTGGCTACAGCGCCGATTCCAGCAAGCTGGCCCAGCGCATCGACCAGGTAGCCCCGACCAAACGCGCCCCGCAGGTGCGCAAGCTGGCCGATGGCAAg                                                                                         >  4:91723/1‑136 (MQ=255)
                                               acGGGCGGCATACTCACCCCGGTCGGGCGGTTGGCGGGCTACAGCGCCGATGCCTGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGa                                                                                     <  5:170285/136‑1 (MQ=255)
                                                      gCATTCTCCCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGCCCAGGGAGCCCCGACCCAAGGCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGCccccccg                                                                              >  3:19549/1‑135 (MQ=255)
                                                              accggggtctgtcTGTTGGCTGGATACAGCGCGTATGCCAGAAAGCTGGCCAAGCGCTTCGACGCTGTATACCGTACAAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGAcc                                                                      <  8:110765/131‑1 (MQ=255)
                                                              aCCCCGGGCGGGCTGGTGGGTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGAcc                                                                      <  2:63460/136‑1 (MQ=255)
                                                                                  tggctACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGCCCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCCCCCTCCACGCCCGGCAGCGCGTGAAGGTGCTg                                                  >  3:5780/1‑136 (MQ=255)
                                                                                      tACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTggcgg                                              >  8:152413/1‑136 (MQ=255)
                                                                                       accgcgcCGAGGCCAGCAAGCGGGCCCAGCGCTTCGACCAGGAAGACCCGACCAAAGGCGCCGAGCGGGTGCGCAAGCGGGCCGATGGCAAGCGGACCACCCACCACGACCGGCAGCGCTTGAAGGTGCTggcggc                                             <  5:81190/133‑1 (MQ=255)
                                                                                       aCAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGAAGCCCCGCCCAAAGGCGCCGAGCTGGTGCGCAAGCGGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTggcggc                                             <  5:158094/136‑1 (MQ=255)
                                                                                        ccgcgcCGATGCCAGCACGCTGGCCCAGCGCTTCGCCCAGGTAGACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGAGGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCa                                            <  1:56359/134‑1 (MQ=255)
                                                                                            gcCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTcc                                        <  8:64003/136‑1 (MQ=255)
                                                                                            gcCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTcc                                        <  7:167294/136‑1 (MQ=255)
                                                                                              cGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGCCCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCgg                                      >  8:105170/1‑136 (MQ=255)
                                                                                                          gcgggCCCCGCGCTGCGCCCAGTTACACCCGACCAACGGCGCCGAGCTGGGGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          <  5:91920/133‑1 (MQ=255)
                                                                                                          gcgggCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCGGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          <  8:151113/133‑1 (MQ=255)
                                                                                                          gCTGGCCCAGCGCTTCGACCAGGTAGCCCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          <  7:105170/136‑1 (MQ=255)
                                                                                                          gCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          >  4:82757/1‑136 (MQ=255)
                                                                                                          gCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          <  4:120962/136‑1 (MQ=255)
                                                                                                          gCTGGCCCAGCGCTGCCACCCGGTACACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          <  1:41189/136‑1 (MQ=255)
                                                                                                           cTGCCCCCGCGCTTCGACCAGGTCGCCCCCCCCAACGGCGCCGAGCTGGGGCGCAAGCTGGCCGATGGAAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                          <  6:119705/135‑1 (MQ=255)
                                                                                                            ggccccAGCGCTGCTACCTGGTAGCCCCGACCAAGGGCGCCGAGCTGGTGCGCAAGCTGGCCGTTGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGc                        <  4:184244/133‑1 (MQ=255)
                                                                                                               cccAGCGGTGCGACCAGCTAGCCCCGCCCAACGGCGCCGAGCGGGGGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCa                     <  8:85193/136‑1 (MQ=255)
                                                                                                                ccAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAg                    >  6:104716/1‑136 (MQ=255)
                                                                                                                        tCGACCTGGTAGACCCGACCAACGGCGCCGAGCTGGTGCGAAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTggg             <  3:91723/135‑1 (MQ=255)
                                                                                                                           aCCAGCTAGACCCGACCAAAGGCGCCGAGCTTGTGCGCAAGCTGGCCGATGGCAAGCTGACCACACTCCACGACCGGCAGCCCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCAGGGGCGt         >  1:125575/1‑136 (MQ=255)
                                                                                                                                  aGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGTGctgctg  >  1:13374/1‑136 (MQ=255)
                                                                                                                                  aGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGCGCGTGctgctg  >  1:86921/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                     
AGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGTGCTGCTG  >  ACB122/4029331‑4029596

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