Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,640,493 C→T 8.0% R175* (CGA→TGA)  PP_3207 → cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,640,4930CT 8.0% 138.4 / 2.2 50R175* (CGA→TGA) PP_3207cyclase
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (2/2);  total (25/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.52e-02

GCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGATATCGGCTATGGGCAGATGGAGA  >  NC_002947/3640349‑3640640
                                                                                                                                                 |                                                                                                                                                   
gCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTGGTACCTGCTGGAGCGAgg                                                                                                                                                 >  7:243183/1‑149 (MQ=255)
         gcgcGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAgg                                                                                                                                          <  6:114862/147‑1 (MQ=255)
           gcgTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACAACCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGggtgg                                                                                                                                      >  4:101799/1‑149 (MQ=255)
                   ggCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTAcc                                                                                                                              <  3:50029/149‑1 (MQ=255)
                      aTGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGAt                                                                                                                           <  7:279332/149‑1 (MQ=255)
                      aTGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGAt                                                                                                                           >  5:228013/1‑149 (MQ=255)
                            tGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGggtggt                                                                                                                                     <  5:128666/133‑1 (MQ=255)
                            tGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGggtggt                                                                                                                                     >  6:128666/1‑133 (MQ=255)
                             gCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTg                                                                                                                      >  2:43693/1‑147 (MQ=255)
                                gCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTAcc                                                                                                                              >  1:301936/1‑136 (MQ=255)
                                gCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTAcc                                                                                                                              <  2:301936/136‑1 (MQ=255)
                                   gCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTgg                                                                                                               >  5:131832/1‑148 (MQ=255)
                                      ggACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATg                                                                                                           >  1:303931/1‑149 (MQ=255)
                                       gACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATgc                                                                                                          <  1:130633/149‑1 (MQ=255)
                                        aCATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATgcg                                                                                                         <  3:147894/149‑1 (MQ=255)
                                         cATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATgcgc                                                                                                        >  3:319956/1‑149 (MQ=255)
                                          aTCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATgcgc                                                                                                        >  7:235716/1‑148 (MQ=255)
                                           tCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCa                                                                                                      <  2:303931/149‑1 (MQ=255)
                                              tGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCAtt                                                                                                    >  5:83445/1‑148 (MQ=255)
                                                cTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCaa                                                                                                 <  6:131832/149‑1 (MQ=255)
                                                cTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCaa                                                                                                 <  1:43693/149‑1 (MQ=255)
                                                cTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGGCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCaa                                                                                                 <  8:235716/149‑1 (MQ=255)
                                                 ttGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTTGTCGGTACCGATGCCTGGAGTTGGGATGCGCAATTCaat                                                                                                >  4:96552/1‑148 (MQ=255)
                                                 ttGTGAACACCCGCGAGGGCAGCTTGTTCGGCCAGCCTGGCTATATCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGTATGAGCAATTCAAg                                                                                                >  4:91849/1‑149 (MQ=255)
                                                   gtgAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGtt                                                                                              >  2:2769/1‑149 (MQ=255)
                                                      aaCACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCAc                                                                                           >  8:233305/1‑149 (MQ=255)
                                                        cacCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGTGAGGGGTGAGGGTGGTCGGTACCGATGCCTg                                                                                                                      <  8:192282/120‑1 (MQ=255)
                                                        cacCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGTGAGGGGTGAGGGTGGTCGGTACCGATGCCTg                                                                                                                      <  8:187352/120‑1 (MQ=255)
                                                        cacCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGTGAGGGGTGAGGGTGGTCGGTACCGATGCCTg                                                                                                                      >  7:192282/1‑120 (MQ=255)
                                                        cacCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGTGAGGGGTGAGGGTGGTCGGTACCGATGCCTg                                                                                                                      >  7:187352/1‑120 (MQ=255)
                                                                    cAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCa                                                                             <  7:233305/149‑1 (MQ=255)
                                                                       cTTGTTCGGCCAGCCTGGCTATCTCGAAGGCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTTAGGGTGGTCGGTACCCATCCGTGGAGTTGGGATGCGCCATTCAAGTACACCCGCGAACGCTACAGt                                                                           >  1:21042/1‑149 (MQ=255)
                                                                           ttCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGaaa                                                                       <  1:2769/148‑1 (MQ=255)
                                                                                      tGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATACCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGCCGCCt                                                           >  3:14254/1‑149 (MQ=255)
                                                                                                gAA‑GCCGTAGTGGGGATGGGCCGCGAGGCTACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCa                                                     <  3:91849/145‑1 (MQ=255)
                                                                                                gAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCAttt                                                  <  1:134321/148‑1 (MQ=255)
                                                                                                gAA‑GCCGGAGGGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCAATCAAGTTCACCCGCGAACGCTTAAGTGAAACGGGCGACGCCTCGATCa                                                     <  3:96552/145‑1 (MQ=255)
                                                                                                                  gggCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCg                                >  7:180852/1‑149 (MQ=255)
                                                                                                                   ggCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCgg                               <  3:101799/149‑1 (MQ=255)
                                                                                                                    gCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCggg                              >  6:279035/1‑149 (MQ=255)
                                                                                                                     ccGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGc                             <  4:319956/149‑1 (MQ=255)
                                                                                                                       gcgAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGt                           >  3:19491/1‑149 (MQ=255)
                                                                                                                        cgAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTg                          <  8:180852/149‑1 (MQ=255)
                                                                                                                            gCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGatat                      <  5:279035/149‑1 (MQ=255)
                                                                                                                               aCCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGATATCgg                   >  3:78660/1‑149 (MQ=255)
                                                                                                                                       ccTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGATATCGGCTATggg            >  5:313648/1‑148 (MQ=255)
                                                                                                                                        ctgctgGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCgg                               >  6:98207/1‑128 (MQ=255)
                                                                                                                                        ctgctgGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCgg                               <  5:98207/128‑1 (MQ=255)
                                                                                                                                            tgGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGATATCGGCTATGGGCAGATg      <  4:19491/149‑1 (MQ=255)
                                                                                                                                                gCGAGGTGTGAGGGTGGGCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGATATCGGCTATGGGCAGATGgaga  <  2:21042/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAA‑GCCGGAGTGGGGATGGGCCGCGAGGCGACCTTGTACCTGCTGGAGCGAGGGGTGAGGGTGGTCGGTACCGATGCCTGGAGTTGGGATGCGCCATTCAAGTTCACCCGCGAACGCTTCAGTGAAACGGGCGACGCCTCGATCATTTGGGAGGGCCATAAGGCCGGGCGTGATATCGGCTATGGGCAGATGGAGA  >  NC_002947/3640349‑3640640

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