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

GGGAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCAATGCCACCC  >  ACB122/3255767‑3256030
                                                                                                                                   |                                                                                                                                          
gggAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGggc                                                                                                                                        >  1:50014/1‑135 (MQ=255)
gggAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGggc                                                                                                                                        >  7:7218/1‑135 (MQ=255)
       ggCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GACGGCAGGGCCGGccca                                                                                                                                 >  7:429994/1‑134 (MQ=255)
        gCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCCGCCGGGCCGGCCC‑tt                                                                                                                               >  1:115657/1‑135 (MQ=255)
         cGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑t                                                                                                                                >  3:206457/1‑133 (MQ=255)
           agaTTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCAGGCCC‑ttt                                                                                                                              >  3:99183/1‑133 (MQ=255)
                gTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGACCC‑TTTC‑G‑CGGGt                                                                                                                     >  3:370362/1‑135 (MQ=255)
                         cGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCCGGGGCGGCCC‑GTTC‑GCCGGGTGGA‑ccgct                                                                                                            >  1:414278/1‑132 (MQ=21)
                                             cgcAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTA‑G‑CGGGTGCA‑CCC‑GCTCCTACGACGGCTTGTGCCgc                                                                                      >  8:330240/1‑134 (MQ=255)
                                                  cTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGCA‑CCC‑GCTCCGACGACGGCTTGTGCCGGacccc                                                                                 >  7:302644/1‑133 (MQ=255)
                                                  cTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑GTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGcaccac                                                                                 >  8:439075/1‑135 (MQ=255)
                                                  cTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑GTTC‑G‑CGGGTGAA‑CCC‑GCGCCGACGCCGGCGTGTGGCGCacccc                                                                                 >  7:28923/1‑133 (MQ=17)
                                                   tCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑GTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGCaccacc                                                                                >  6:226599/1‑135 (MQ=255)
                                                          tttGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGCA‑CCC‑GCTCCTACGACGGCTTGTGCCGCACCCCCGATGgac                                                                         >  1:39760/1‑135 (MQ=255)
                                                             gTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑GTTC‑G‑CGGGTGAA‑CCC‑GCGCCTACGGCGGCTTGTGGCGCCCCCCCGATGgacgcc                                                                      >  5:454439/1‑133 (MQ=14)
                                                               aGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCCGGGCCGGCCC‑GTTC‑G‑CGGGGGAG‑CCC‑GCTCCTACGTCGGCGTGTGCCGCACCACCGAGGgccgcc                                                                      >  7:18697/1‑131 (MQ=14)
                                                               aGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑GTTC‑G‑CGGGTGGG‑CCC‑GCGCCTACGACGGGGTTTGGCGCACCACCGCTGgccggc                                                                      >  8:310215/1‑131 (MQ=14)
                                                                     gggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTGCTACGACGGCTTGTGCCGAACCCCCGCTGGCCGCCGTAGGc                                                                >  4:128749/1‑133 (MQ=17)
                                                                      ggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCCGGGCCGGCCC‑TTTC‑G‑CGGGTGCA‑CCC‑GCGGCGTCGCCGGCTTGTGGCGAACCACCGATGGCCGACGTAGGCTCa                                                             >  6:67635/1‑135 (MQ=17)
                                                                       gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGCA‑CCC‑GCTCCGACGCCGGGTTGTGCCGCACCACCGATGGCCGCCGGAGGGTCAg                                                            >  4:451559/1‑135 (MQ=17)
                                                                       gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGAACCACCGCTGGACGACGTAGGCTCAg                                                            >  8:26425/1‑135 (MQ=255)
                                                                         ttACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCCGGGCCGGCCC‑GTTC‑G‑CGGGTGGA‑CCC‑GCTGCGACGCCGGCGTGTGCCGGACCCCCGAGGGCCGGCGTAGGGGCAGAc                                                          >  1:55246/1‑135 (MQ=14)
                                                                         ttACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTA‑G‑CGGGTGAA‑CCC‑GGTCCGACGCCGGCTTGTGCCGCACCCCCGATGGACGGCGTAGGGTGAGAc                                                          >  7:461415/1‑135 (MQ=17)
                                                                                    tACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑GTTG‑G‑CGGGGGGA‑CCC‑GCGCCTACGACGGGTTGTGCCGCACCCCCGCGGGGCGCCGTTGGCTCCGCCCCGGCACTGcc                                               >  6:318545/1‑135 (MQ=17)
                                                                                            cggcggCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCCGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGAAACAGCGATGGACGACCTAGGCGCAGACCCCGAACTGCCGcacccac                                       >  8:152810/1‑130 (MQ=21)
                                                                                               cggcAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGcc                                                                                        >  8:464447/1‑83 (MQ=255)
                                                                                               cggcAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGcc                                                                                        <  7:464447/83‑1 (MQ=255)
                                                                                                    aCGGGGCACGGGGGGGGGAAATTG‑GCCCGCCGGGGCGGGCC‑GTTCGG‑CGGGTCAC‑CCG‑GCCCCAACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCa                                <  7:439075/135‑1 (MQ=255)
                                                                                                    aCGGGGACCGGGGGGGGGAAAATG‑GCCGGCCGGGCCGGGCC‑GTTT‑G‑GGGGTTAACCCC‑GCGCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCa                                <  7:330240/135‑1 (MQ=255)
                                                                                                             ggcggGGGGAAATTG‑GCCGGCCGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGCCGGCTTGTGCCGCACCCCCGAAGGGCGACGTAGGCTCAGGCCCGGAACTGCCGCAACAACTGCTCCAGCGCTTAGca                      >  2:288356/1‑134 (MQ=25)
                                                                                                              gggggggcAATTTG‑GCCGCCCGGGCGGCCCC‑GTTC‑G‑CGGGTCAC‑CCG‑GCCCCAACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTc                     <  3:338241/133‑1 (MQ=255)
                                                                                                              ggcgggggAAATTGCGCGGGCCGGGCGGGCCC‑TTGC‑G‑GGGGGGAA‑CCC‑GCCCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCt                      <  6:454439/132‑1 (MQ=255)
                                                                                                              gcggGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTc                     <  1:288356/135‑1 (MQ=255)
                                                                                                               gggggggAAACTG‑GCCGGGCGGGCCGGCCC‑CTTG‑G‑GGGGGGCA‑CCCGGCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTc                     <  2:249372/134‑1 (MQ=255)
                                                                                                                     gAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGcc              >  8:395378/1‑135 (MQ=255)
                                                                                                                          gg‑gccggccGGGGCGGCCCGTTGC‑G‑GGGGGAAC‑CCC‑GCCCCAACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCaa          <  8:302644/134‑1 (MQ=255)
                                                                                                                                  cAGGGCCGTCCC‑TTTC‑C‑CGGGTGAA‑CCC‑GATCCTACGACGGCTTGTGCCGAACCACCTATGGACGACGTAGGCTCAGACCCGGAACTGCAGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCAATGCCAcca  <  8:251244/135‑2 (MQ=255)
                                                                                                                                   |                                                                                                                                          
GGGAAGGGGCGAGATTGTAGCGGGTCGGCCGGTGAAGGTGTAGGGCGCAACTCAACCCTTTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGGGGAAATTG‑GCCGGCAGGGCCGGCCC‑TTTC‑G‑CGGGTGAA‑CCC‑GCTCCTACGACGGCTTGTGCCGAACCACCGATGGACGACGTAGGCTCAGACCCGGAACTGCCGCAACAACTGCTCCAGCTCTTCGCTCAACCGCCCCAATGCCACCC  >  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: ‑