Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,255,897 A→C 28.4% intergenic (‑152/+72) pnuC ← / ← mcpP Nicotinamide riboside transporter PnuC/Methyl‑accepting chemotaxis protein McpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,255,8970AC28.4% 22.6 / 4.0 36intergenic (‑152/+72)pnuC/mcpPNicotinamide riboside transporter PnuC/Methyl‑accepting chemotaxis protein McpP
Reads supporting (aligned to +/- strand):  ref base A (22/3);  new base C (6/5);  total (28/8)
Fisher's exact test for biased strand distribution p-value = 4.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGGAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCAATGCCACCC  >  ACB122/3255767‑3256030
                                                                                                                                   |                                                                                                                                          
gggAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGggc                                                                                                                                        >  1:50014/1‑135 (MQ=255)
gggAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGggc                                                                                                                                        >  7:7218/1‑135 (MQ=255)
       ggCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGACGGCAGGGCCGGccca                                                                                                                                 >  7:429994/1‑134 (MQ=255)
        gCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCCGCCGGGCCGGCCCtt                                                                                                                               >  1:115657/1‑135 (MQ=255)
         cGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCt                                                                                                                                >  3:206457/1‑133 (MQ=255)
           agaTTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCAGGCCCttt                                                                                                                              >  3:99183/1‑133 (MQ=255)
                gTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGACCCTTTCGCGGGt                                                                                                                     >  3:370362/1‑135 (MQ=255)
                         cGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCCGGGGCGGCCCGTTCGCCGGGTGGAccgct                                                                                                            >  1:414278/1‑132 (MQ=21)
                                             cgcAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTAGCGGGTGCACCCGCTCCTACGACGGCTTGTGCCgc                                                                                      >  8:330240/1‑134 (MQ=255)
                                                  cTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGCACCCGCTCCGACGACGGCTTGTGCCGGacccc                                                                                 >  7:302644/1‑133 (MQ=255)
                                                  cTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCGTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGcaccac                                                                                 >  8:439075/1‑135 (MQ=255)
                                                  cTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCGTTCGCGGGTGAACCCGCGCCGACGCCGGCGTGTGGCGCacccc                                                                                 >  7:28923/1‑133 (MQ=17)
                                                   tCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCGTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGCaccacc                                                                                >  6:226599/1‑135 (MQ=255)
                                                          tttGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGCACCCGCTCCTACGACGGCTTGTGCCGCACCCCCGATGgac                                                                         >  1:39760/1‑135 (MQ=255)
                                                             gTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCGTTCGCGGGTGAACCCGCGCCTACGGCGGCTTGTGGCGCCCCCCCGATGgacgcc                                                                      >  5:454439/1‑133 (MQ=14)
                                                               aGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCCGGGCCGGCCCGTTCGCGGGGGAGCCCGCTCCTACGTCGGCGTGTGCCGCACCACCGAGGgccgcc                                                                      >  7:18697/1‑131 (MQ=14)
                                                               aGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCGTTCGCGGGTGGGCCCGCGCCTACGACGGGGTTTGGCGCACCACCGCTGgccggc                                                                      >  8:310215/1‑131 (MQ=14)
                                                                     gggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTGCTACGACGGCTTGTGCCGAACCCCCGCTGGCCGCCGTAGGc                                                                >  4:128749/1‑133 (MQ=17)
                                                                      ggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCCGGGCCGGCCCTTTCGCGGGTGCACCCGCGGCGTCGCCGGCTTGTGGCGAACCACCGATGGCCGACGTAGGCTCa                                                             >  6:67635/1‑135 (MQ=17)
                                                                       gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGCACCCGCTCCGACGCCGGGTTGTGCCGCACCACCGATGGCCGCCGGAGGGTCAg                                                            >  4:451559/1‑135 (MQ=17)
                                                                       gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGAACCACCGCTGGACGACGTAGGCTCAg                                                            >  8:26425/1‑135 (MQ=255)
                                                                         ttACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCCGGGCCGGCCCGTTCGCGGGTGGACCCGCTGCGACGCCGGCGTGTGCCGGACCCCCGAGGGCCGGCGTAGGGGCAGAc                                                          >  1:55246/1‑135 (MQ=14)
                                                                         ttACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTAGCGGGTGAACCCGGTCCGACGCCGGCTTGTGCCGCACCCCCGATGGACGGCGTAGGGTGAGAc                                                          >  7:461415/1‑135 (MQ=17)
                                                                                    tACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCGTTGGCGGGGGGACCCGCGCCTACGACGGGTTGTGCCGCACCCCCGCGGGGCGCCGTTGGCTCCGCCCCGGCACTGcc                                               >  6:318545/1‑135 (MQ=17)
                                                                                            cggcggCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCCGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGAAACAGCGATGGACGACCTAGGCGCAGACCCCGAACTGCCGcacccac                                       >  8:152810/1‑130 (MQ=21)
                                                                                               cggcAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGcc                                                                                        >  8:464447/1‑83 (MQ=255)
                                                                                               cggcAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGcc                                                                                        <  7:464447/83‑1 (MQ=255)
                                                                                                    aCGGGGCACGGGGGGGGGAAATTGGCCCGCCGGGGCGGGCCGTTCGGCGGGTCACCCGGCCCCAACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCa                                <  7:439075/135‑1 (MQ=255)
                                                                                                    aCGGGGACCGGGGGGGGGAAAATGGCCGGCCGGGCCGGGCCGTTTGGGGGTTAACCCCGCGCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCa                                <  7:330240/135‑1 (MQ=255)
                                                                                                             ggcggGGGGAAATTGGCCGGCCGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGCCGGCTTGTGCCGCACCCCCGAAGGGCGACGTAGGCTCAGGCCCGGAACTGCCGCAACAACTGCTCCAGCGCTTAGca                      >  2:288356/1‑134 (MQ=25)
                                                                                                              gggggggcAATTTGGCCGCCCGGGCGGCCCCGTTCGCGGGTCACCCGGCCCCAACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTc                     <  3:338241/133‑1 (MQ=255)
                                                                                                              ggcgggggAAATTGCGCGGGCCGGGCGGGCCCTTGCGGGGGGGAACCCGCCCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCt                      <  6:454439/132‑1 (MQ=255)
                                                                                                              gcggGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTc                     <  1:288356/135‑1 (MQ=255)
                                                                                                               gggggggAAACTGGCCGGGCGGGCCGGCCCCTTGGGGGGGGCACCCGGCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTc                     <  2:249372/134‑1 (MQ=255)
                                                                                                                     gAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGcc              >  8:395378/1‑135 (MQ=255)
                                                                                                                          gggccggccGGGGCGGCCCGTTGCGGGGGGAACCCCGCCCCAACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCaa          <  8:302644/134‑1 (MQ=255)
                                                                                                                                  cAGGGCCGTCCCTTTCCCGGGTGAACCCGATCCTACGACGGCTTGTGCCGAACCACCTATGGACGACGTAGGCTCAGACCCGGAACTGCAGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCAATGCCAcca  <  8:251244/135‑2 (MQ=255)
                                                                                                                                   |                                                                                                                                          
GGGAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTGGCCGGCAGGGCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCAATGCCACCC  >  ACB122/3255767‑3256030

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: