Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,409,233 G→T 14.3% G290C (GGC→TGC)  PP_2106 → ammonium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,409,2330GT14.3% 80.6 / 3.5 28G290C (GGC→TGC) PP_2106ammonium transporter
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base T (2/2);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

ATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCG  >  NC_002947/2409098‑2409362
                                                                                                                                       |                                                                                                                                 
aTGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCg                                                                                                                                   >  2:179241/1‑136 (MQ=255)
aTGGTGGGTGGCACCCTGTCAGCGTTGCTGGACGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCg                                                                                                                                   >  2:588281/1‑136 (MQ=255)
        gggCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCtggtg                                                                                                                           <  2:355025/135‑1 (MQ=255)
           cACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCg                                                                                                                        >  1:648542/1‑136 (MQ=255)
                   cAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGgcgcg                                                                                                                >  2:45781/1‑135 (MQ=255)
                     gCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGt                                                                                                              >  2:575630/1‑136 (MQ=255)
                                               gggCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCaa                                                                                    <  1:414602/136‑1 (MQ=255)
                                               gggCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCaa                                                                                    <  1:516127/136‑1 (MQ=255)
                                                        gCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCTGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTgg                                                                           >  1:122767/1‑136 (MQ=255)
                                                        gCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCTGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTgg                                                                           >  1:501672/1‑136 (MQ=255)
                                                         cATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCACTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGa                                                                          >  2:476919/1‑136 (MQ=255)
                                                              cGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATc                                                                     <  1:179241/136‑1 (MQ=255)
                                                              cGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATc                                                                     <  1:588281/136‑1 (MQ=255)
                                                              cGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATc                                                                     <  1:440354/136‑1 (MQ=255)
                                                                      tGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGt                                                             <  1:575630/136‑1 (MQ=255)
                                                                      gggCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGt                                                             <  1:247820/135‑1 (MQ=255)
                                                                              cTGGTGGCGATTTGCGCTGGTTCCGATCTGATGCACCCTGTGGGGGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCg                                                     <  1:220835/136‑1 (MQ=255)
                                                                                  tggCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACGGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATg                                                 <  2:644147/136‑1 (MQ=255)
                                                                                     cGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCCGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGcc                                              >  1:433645/1‑136 (MQ=255)
                                                                                       ttttGGGCAGGTTCCTATCTGATGCACCCGGTGGGTGGGCTCGCCTCCGGCCTGGTGGGGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGGTGGAAAATCGACGACGTCCTGGGCGTATGGCCGt                                            <  1:313982/135‑1 (MQ=255)
                                                                                        tttGCGCAGGTTCCTATCTGATGCACCCGGGGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGtt                                           <  2:533131/136‑1 (MQ=255)
                                                                                                     cGATCTGATGCACCGGGTGGGTGGGCTCGCCACGGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCtgtgt                              <  2:648542/136‑1 (MQ=255)
                                                                                                                 ccgggtgggtgCGCTCGCCACCGGCCGGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGt                  <  1:563867/133‑1 (MQ=255)
                                                                                                                 cccGGTGGGTGCGCTCGCCACCTGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGTGGt                  <  2:501672/136‑1 (MQ=255)
                                                                                                                 cccGGTGGGTGCGCTCGCCACCTGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGTGGt                  <  2:122767/136‑1 (MQ=255)
                                                                                                                 cccGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGt                  <  2:51452/136‑1 (MQ=255)
                                                                                                                               tCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGACATCGCCTGTgg    >  2:49083/1‑136 (MQ=255)
                                                                                                                                 gCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCg  >  2:372499/1‑136 (MQ=255)
                                                                                                                                 gCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCg  >  1:491816/1‑136 (MQ=255)
                                                                                                                                 gCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCg  >  1:510618/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                 
ATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCG  >  NC_002947/2409098‑2409362

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: ‑