Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,953,938 Δ1 bp 7.1% coding (585/1977 nt) PP_5606 ← YVTN beta‑propeller repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,953,9380C. 7.1% 181.3 / 7.4 56coding (585/1977 nt)PP_5606YVTN beta‑propeller repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (29/23);  new base . (2/2);  total (31/25)
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.06e-01

ATCACTGTTACCCCACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTG  >  NC_002947/3953798‑3954083
                                                                                                                                             |                                                                                                                                                 
aTCACTGTTACCCCACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGAc                                                                                                                                           <  8:86385/149‑1 (MQ=255)
aTCACTGTTACCCCACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGAc                                                                                                                                           >  3:288798/1‑149 (MQ=255)
           cccACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCt                                                                                                                                >  4:116840/1‑149 (MQ=255)
             cACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCAc                                                                                                                                                 >  6:9035/1‑130 (MQ=255)
             cACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCAc                                                                                                                                                 <  5:9035/130‑1 (MQ=255)
                    ggCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATc                                                                                                                        <  4:327587/148‑1 (MQ=255)
                     gCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGa                                                                                                                      >  2:379219/1‑149 (MQ=255)
                               tgctgTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAAcc                                                                                                            >  8:317182/1‑149 (MQ=255)
                                gctgTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCt                                                                                                           >  4:184941/1‑149 (MQ=255)
                                          ccAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTgtcgtc                                                                                                 >  8:350457/1‑149 (MQ=255)
                                            aaCCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGc                                                                                               >  1:38847/1‑149 (MQ=255)
                                                     ccGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCAc                                                                                                                                                 <  8:88877/90‑1 (MQ=255)
                                                     ccGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCAc                                                                                                                                                 >  7:88877/1‑90 (MQ=255)
                                                       gATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGgccgcc                                                                                    <  4:288798/149‑1 (MQ=255)
                                                         tGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCa                                                                                  <  3:184941/149‑1 (MQ=255)
                                                           cTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGt                                                                                <  2:38847/149‑1 (MQ=255)
                                                            tGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTc                                                                               >  6:75042/1‑149 (MQ=255)
                                                            tGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCAGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTc                                                                               >  4:322588/1‑149 (MQ=255)
                                                             ggTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTc                                                                               >  6:230169/1‑148 (MQ=255)
                                                             ggTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTc                                                                               >  7:117377/1‑148 (MQ=255)
                                                              gTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTAAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGt                                                                                >  5:270758/1‑146 (MQ=255)
                                                              gTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTAAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGt                                                                                <  6:270758/146‑1 (MQ=255)
                                                               tCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCAt                                                                            >  7:42145/1‑149 (MQ=255)
                                                                 cGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATgg                                                                          <  6:234055/149‑1 (MQ=255)
                                                                 cGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATg                                                                           <  3:116840/148‑1 (MQ=255)
                                                                  ggCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATggg                                                                         <  7:350457/149‑1 (MQ=255)
                                                                      gCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCtggtgg                                                                                         >  2:324840/1‑129 (MQ=255)
                                                                      gCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGgccgc                                                                     >  6:184427/1‑149 (MQ=255)
                                                                      gCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGgccgc                                                                     >  7:308427/1‑149 (MQ=255)
                                                                        aGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGccc                                                                   <  5:75042/149‑1 (MQ=255)
                                                                           gCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTGCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTgg                                                                <  5:85638/149‑1 (MQ=255)
                                                                            ccACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTgg                                                                >  4:322258/1‑148 (MQ=255)
                                                                                gcgcAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCt                                                           <  8:42145/149‑1 (MQ=255)
                                                                                 cgcAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTc                                                          <  8:308427/149‑1 (MQ=255)
                                                                                  gcAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTc                                                          >  1:376116/1‑148 (MQ=255)
                                                                                     cgaccTGGACGCCGGCGACCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATg                                                      >  8:106812/5‑149 (MQ=255)
                                                                                      cGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATgc                                                     <  3:311477/149‑1 (MQ=255)
                                                                                       accTGGACGCCGGCGACCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGTGCCGCCCTGGCAGCTCGATgcg                                                    >  3:130783/3‑149 (MQ=255)
                                                                                        gCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATgcgc                                                   <  7:317182/149‑1 (MQ=255)
                                                                                                         aGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCAc                                  >  4:137282/1‑149 (MQ=255)
                                                                                                         aGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCAc                                  >  3:110919/1‑149 (MQ=255)
                                                                                                                ccgagcgGGGCAAGGTCTCGGTGATCGATACCGAAACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGccccc                            >  4:257293/5‑148 (MQ=38)
                                                                                                                ccgagcgGGGCAAGGTCTCGGTGATCGATACCGAAACCTTCACCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGccccc                            >  6:164534/5‑148 (MQ=37)
                                                                                                                          gAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTg                                       >  4:298802/1‑127 (MQ=255)
                                                                                                                          gAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTg                                       <  3:298802/127‑1 (MQ=255)
                                                                                                                              gTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCTCCACCGACACCGTCgg             >  1:46799/1‑149 (MQ=255)
                                                                                                                                 tCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGt           <  4:110919/148‑1 (MQ=255)
                                                                                                                                 tCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGt           <  5:184427/148‑1 (MQ=255)
                                                                                                                                   ggtgaTCGATACCGAAACCTTCACCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGt           <  5:164534/142‑1 (MQ=255)
                                                                                                                                    gtgaTCGATACCGAAACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGAt       <  3:257293/146‑1 (MQ=255)
                                                                                                                                       aTCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGCCACGcccc                             >  7:316521/1‑125 (MQ=255)
                                                                                                                                       aTCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGCCACGcccc                             <  8:316521/125‑1 (MQ=255)
                                                                                                                                         cGATCACCGAGACCTTGCCCCGCTCGGGTAACGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACt   >  1:30136/1‑149 (MQ=255)
                                                                                                                                          gATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTg  <  3:137282/149‑1 (MQ=255)
                                                                                                                                            tCACCGTGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACt   <  3:67657/146‑1 (MQ=255)
                                                                                                                                            tCACCGTGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACt   >  4:67657/1‑146 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
ATCACTGTTACCCCACCTTTGGCCGGGTCGCTGCTGTTGTTGCCAACCCACAGCCGATGCTGGTCCGGCTGCAGGGCCACGCGCAGCGGCTGGTCGCCGGCGTCCAGGTCGGCCACGCGGGTGAAGGTTTCGGTATCGATCACCGAGACCTTGCCCCGCTCGGGTATCGACACGAAAACCTGCTTGTCGTCGCTGGTGGCCGCCCAGTCCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTG  >  NC_002947/3953798‑3954083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: