Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,209,603 G→A 7.0% A262T (GCC→ACC)  PP_2813 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,209,6030GA 7.0% 164.7 / 2.6 57A262T (GCC→ACC) PP_2813BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CTGTGGGCAGGCGCGGTACTGGCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTT  >  NC_002947/3209462‑3209750
                                                                                                                                             |                                                                                                                                                   
cTGTGGGCAGGCGCGGTACTGGCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCAACCGCAAAGCAAGGGCTCTGGTACGTGCGTCGCGGTAGCCGGctct                                                                                                                                             >  2:443266/1‑150 (MQ=255)
 tgtgGGCAGGCGCGGTACTGGCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCtctc                                                                                                                                            <  2:223323/150‑1 (MQ=255)
     ggCAGGCGCGGTACTGGCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCgtgt                                                                                                                                        <  1:450979/150‑1 (MQ=255)
         ggCGCGGTACTGGCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTg                                                                                                                                    <  1:785711/150‑1 (MQ=255)
                    ggCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCt                                                                                                                         <  1:249616/150‑1 (MQ=255)
                    ggCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCt                                                                                                                         >  2:1123781/1‑150 (MQ=255)
                     gCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCTGTCGGGCTg                                                                                                                        >  2:279198/1‑150 (MQ=255)
                                 ctgctgctgGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGc                                                                                                            <  2:98125/150‑1 (MQ=255)
                                  tgctgctgGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCtctc                                                                                                                                            >  2:996642/1‑117 (MQ=255)
                                  tgctgctgGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCtctc                                                                                                                                            <  1:996642/117‑1 (MQ=255)
                                     tgctgGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTg                                                                                                                        <  1:121501/134‑1 (MQ=255)
                                     tgctcGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTATCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTg                                                                                                                        >  2:121501/1‑134 (MQ=255)
                                        tgGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCaga                                                                                                     <  1:1418310/150‑1 (MQ=255)
                                        tgGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCaga                                                                                                     <  1:1159131/150‑1 (MQ=255)
                                          ggcggcCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGACTCCGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGagc                                                                                                   <  1:1168372/150‑1 (MQ=255)
                                          ggcggcCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGACTGCCAACGCAAAGCAAGGTCGCTGTCACGTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGGAGagc                                                                                                   >  1:883580/1‑150 (MQ=255)
                                              gcCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCg                                                                                               >  2:995239/1‑150 (MQ=255)
                                                cAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCgg                                                                                                                            <  1:1333497/119‑1 (MQ=255)
                                                cAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCgg                                                                                                                            >  2:1333497/1‑119 (MQ=255)
                                                             tGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGc                                                                                                                                                    <  1:439878/82‑1 (MQ=255)
                                                             tGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGc                                                                                                                                                    >  2:439878/1‑82 (MQ=255)
                                                                 gCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCa                                                                            >  2:1415845/1‑150 (MQ=255)
                                                                 gCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCa                                                                            <  2:3030/150‑1 (MQ=255)
                                                                            cGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGa                                                                 <  1:1104753/150‑1 (MQ=255)
                                                                             ggACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGc                                                                                                                                                    >  2:364620/1‑66 (MQ=255)
                                                                             ggACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGc                                                                                                                                                    <  1:364620/66‑1 (MQ=255)
                                                                                 ggCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCgt                                                                                                                                          <  1:619321/72‑1 (MQ=255)
                                                                                 ggCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCgt                                                                                                                                          >  2:619321/1‑72 (MQ=255)
                                                                                 ggCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACcgcg                                                            >  1:1396009/1‑150 (MQ=255)
                                                                                    aaGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGAt                                                         <  2:883580/150‑1 (MQ=255)
                                                                                          tttCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCtt                                                   >  2:281369/1‑150 (MQ=255)
                                                                                          tttCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCtt                                                   >  2:253640/1‑150 (MQ=255)
                                                                                                       cAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCTGTCGGGCTGGAGGGGCTTAGCGTGCAGAGCAGCGATGACGGGCACAGTTTCCAGGAGACGCAGATCGGGTATGGCCTTTCACTGCcggcgg                                      >  2:290810/1‑150 (MQ=255)
                                                                                                       cAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCTCAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCa                                                <  1:1154857/140‑1 (MQ=255)
                                                                                                       cAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCa                                                >  2:1154857/1‑140 (MQ=255)
                                                                                                       cAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGAcggcgg                                      >  1:1062591/1‑150 (MQ=255)
                                                                                                        aaCGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGATCGGGCTGGATGGGc                                                                                                                 <  1:1388395/74‑1 (MQ=255)
                                                                                                        aaCGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGATCGGGCTGGATGGGc                                                                                                                 >  2:1388395/1‑74 (MQ=255)
                                                                                                        aaCGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAACCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGcaca                                                                                   >  1:373085/1‑104 (MQ=255)
                                                                                                        aaCGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAACCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGcaca                                                                                   >  1:373222/1‑104 (MQ=255)
                                                                                                        aaCGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAACCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGcaca                                                                                   <  2:373085/104‑1 (MQ=255)
                                                                                                        aaCGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAACCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGcaca                                                                                   <  2:373222/104‑1 (MQ=255)
                                                                                                           gCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCgtt                                  <  1:541089/150‑1 (MQ=255)
                                                                                                             aaaGCAAGGTCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGCGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGCAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCgttgt                                <  2:714662/150‑1 (MQ=255)
                                                                                                             aaaGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCgttgt                                <  1:1415845/150‑1 (MQ=255)
                                                                                                             aaaGCAAGGGCGTTCTCACTTGCGGAGCGGTAGCCGGCTCTCGTGTGCTGGCGGGCGGGCTGGATGGGCTCAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCgttgt                                >  2:54066/1‑150 (MQ=255)
                                                                                                                        tctGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGaa                     >  1:536859/1‑150 (MQ=255)
                                                                                                                        tctGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGaa                     >  1:537168/1‑150 (MQ=255)
                                                                                                                          tGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCCTTCACTGACGGCGGCGTTGTTCAACGCCGAAgg                   >  1:758908/1‑150 (MQ=255)
                                                                                                                            tCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATTGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGAt                                                         <  1:983691/110‑1 (MQ=255)
                                                                                                                            tCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGAt                                                         >  2:983691/1‑110 (MQ=255)
                                                                                                                               cTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTAc              >  1:295864/1‑150 (MQ=255)
                                                                                                                               cTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTAc              >  1:311480/1‑150 (MQ=255)
                                                                                                                               cTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTCAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTAc              <  1:54066/150‑1 (MQ=255)
                                                                                                                                      aGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTAtt       <  1:995239/150‑1 (MQ=255)
                                                                                                                                          gTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGttt   <  1:909939/150‑1 (MQ=255)
                                                                                                                                           tAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGtttt  <  1:1123781/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTGTGGGCAGGCGCGGTACTGGCCGATGGACGTCTGCTGCTGGGCGGCCAGCGCGGTACCCTGTTGCAGGGTTCAGCGGACGGCAAGGGCTTTCACCGAGTGCCAACGCAAAGCAAGGGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTT  >  NC_002947/3209462‑3209750

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