Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,731,574 T→C 10.9% T780A (ACG→GCG) ‡ PP_1526 ← beta‑(1‑3)‑glucosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,731,5740TC10.9% 118.4 / 3.7 44T780A (ACG→GCG) ‡PP_1526beta‑(1‑3)‑glucosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (2/3);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAGA  >  NC_002947/1731428‑1731711
                                                                                                                                                  |                                                                                                                                         
cACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCg                                                                                                                                        >  2:964937/1‑150 (MQ=255)
cACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCg                                                                                                                                        >  2:181955/1‑150 (MQ=255)
          cgcATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCCACCCATGGCTGTCAGCGTTCTTCCGCACGCCGAAGAACgct                                                                                                                              <  1:1168727/150‑3 (MQ=255)
                tCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCAtgctg                                                                                                                        >  2:511780/1‑150 (MQ=255)
                  cGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATgctgct                                                                                                                      >  1:190702/1‑150 (MQ=255)
                    aaCTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTgg                                                                                                                    <  1:964937/150‑1 (MQ=255)
                      cTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCgaagaa                                                                                                                                  >  2:1173621/1‑134 (MQ=255)
                      cTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCgaagaa                                                                                                                                  <  1:1173621/134‑1 (MQ=255)
                             gCCCCTGCACCAGGTAGATACCCGACGCGGCTCGCCACGTCAGCACCATGCTGCACAGCTCTTCAGCTGCATCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAAc                                                                                                           <  1:1095938/150‑1 (MQ=255)
                             gCCCCTGCACCAGGCACCAACCCAACGCCGAACACAACCGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAAc                                                                                                           <  2:1142845/150‑1 (MQ=255)
                               cccTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCa                                                                                                         <  2:190702/150‑1 (MQ=255)
                                        aGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCCGCACGCCgaagaa                                                                                                                                  <  2:1401562/116‑1 (MQ=255)
                                             gATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTga                                                                                                                 >  1:1367208/1‑128 (MQ=255)
                                             gATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTga                                                                                                                 <  2:1367208/128‑1 (MQ=255)
                                               tACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTgg                                                                                         <  1:40961/150‑1 (MQ=255)
                                                   cAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGAt                                                                                     >  2:607686/1‑150 (MQ=255)
                                                    aaCGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATg                                                                                    <  1:807709/150‑1 (MQ=255)
                                                    aaCGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATg                                                                                    >  2:218976/1‑150 (MQ=255)
                                                    aaCGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATg                                                                                    >  2:1162149/1‑150 (MQ=255)
                                                      cgccgcGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGAAACTCTCGGTTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGgt                                                                                  <  1:449616/150‑1 (MQ=255)
                                                              cccACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTCGCGATGGTGTGCGACa                                                                          >  2:1382143/1‑150 (MQ=255)
                                                                      agcaCCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCCTGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGc                                                                  <  1:218976/150‑1 (MQ=255)
                                                                        caCCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGccc                                                                <  1:280728/150‑1 (MQ=255)
                                                                        caCCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGccc                                                                >  2:1393077/1‑150 (MQ=255)
                                                                             tgatgaACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCa                                                           <  1:511780/150‑1 (MQ=255)
                                                                              gatgaACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAg                                                          >  2:1416278/1‑150 (MQ=255)
                                                                                 gaACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAgcgc                                                       >  1:1373999/1‑150 (MQ=255)
                                                                                     aGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCt                                                   >  2:870797/1‑150 (MQ=255)
                                                                                     aGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCt                                                   <  1:181955/150‑1 (MQ=255)
                                                                                     aGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCt                                                   >  1:498336/1‑150 (MQ=255)
                                                                                       ctctTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCt                                                   >  2:516261/1‑148 (MQ=255)
                                                                                       ctctTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCt                                                   <  1:516261/148‑1 (MQ=255)
                                                                                       ctctTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGc                                                 >  2:928761/1‑150 (MQ=255)
                                                                                                     cGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGa                                   <  1:1162149/150‑1 (MQ=255)
                                                                                                     cGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGa                                   <  1:1393077/150‑1 (MQ=255)
                                                                                                              ccaccaGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAccc                          <  1:502109/150‑1 (MQ=255)
                                                                                                                accaGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAccccc                        >  1:1309219/1‑150 (MQ=255)
                                                                                                                  caGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAc                            >  1:124889/1‑144 (MQ=255)
                                                                                                                  caGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAc                            <  2:124889/144‑1 (MQ=255)
                                                                                                                         ccATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAc                            >  1:577715/1‑137 (MQ=255)
                                                                                                                         ccATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCAc                            <  2:577715/137‑1 (MQ=255)
                                                                                                                           aTGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTcggcg             <  1:928761/150‑1 (MQ=255)
                                                                                                                                    aGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGaa    >  2:1151193/1‑150 (MQ=255)
                                                                                                                                      cGTTCTTCCGCACGCCGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAga  >  2:131822/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
CACCCAGAAGCGCATGTCCGAACTCGGCAGCCCCTGCACCAGGTAGATACCCAACGCCGCGCCCCACAGCAGCACCATGATGAACAGCTCTTCACGGGCTTCGGAAATCGCCACCAGCAACCCATGGCTGTCAGCGTTCTTCGGCGTGCGGAAGAACGGCATGCTGCTGGTGAAGAAACCATACAGTACCGCCTTGGCGATGGTGTGCGACAGTGCCAGCCCGGCCAGCGCAGCTGCGAAGGCATCCTTGAGGTTCACCCCCACCGCTCGGCGGTACAGGAAGA  >  NC_002947/1731428‑1731711

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