Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,416,803 G→A 7.0% N146N (AAC→AAT)  PP_3029 ← DNA cytosine methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,416,8030GA 7.0% 169.2 / 2.1 57N146N (AAC→AAT) PP_3029DNA cytosine methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (2/2);  total (27/30)
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.89e-01

ACCAATGCGTCCATCGGTACCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTG  >  NC_002947/3416657‑3416941
                                                                                                                                                  |                                                                                                                                          
aCCAATGCGTCCATCGGTACCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGtt                                                                                                                                          >  3:270056/1‑149 (MQ=255)
  cAATGCGTCCATCGGTACCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGtt                                                                                                                                          >  5:10802/1‑147 (MQ=255)
   aaTGCGTCCATCGGTACCCACTGCCCGATGAAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGGGGTGAGGCAGTCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGGTTGCCTTCGTCGTTCAg                                                                                                                                       <  2:192428/149‑1 (MQ=255)
                 tACCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGtt                                                                                                                                          <  4:313089/132‑1 (MQ=255)
                 tACCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGtt                                                                                                                                          >  3:313089/1‑132 (MQ=255)
                  aCCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTg                                                                                                                        >  5:338945/1‑149 (MQ=255)
                   cccACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTgg                                                                                                                       <  5:60872/149‑1 (MQ=255)
                   cccACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTg                                                                                                                        >  8:354063/1‑148 (MQ=255)
                   cccACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTg                                                                                                                        <  4:305039/148‑1 (MQ=255)
                   cccACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTg                                                                                                                        <  1:251000/148‑1 (MQ=255)
                     cACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGaa                                                                                                                     <  4:270056/149‑1 (MQ=255)
                          cccGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAActc                                                                                                                  <  7:349082/147‑1 (MQ=255)
                               tGCAGTCCATGACCTAGCTCATGTCGGTATTATCGGCAGGCAAGCCAGTGGTGATGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTTGAACTCCTCAACg                                                                                                           <  4:375612/149‑1 (MQ=255)
                                    tCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTtctc                                                                                                      <  8:124925/149‑1 (MQ=255)
                                       aTGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCAt                                                                                                   <  8:201652/149‑1 (MQ=255)
                                       aTGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGCTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCAt                                                                                                   <  6:10802/149‑1 (MQ=255)
                                          aCCTCGCTCATGTCGGGATGATCGGAAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATca                                                                                                 <  6:122487/148‑1 (MQ=255)
                                            cTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCtt                                                                                                                         >  3:108176/1‑122 (MQ=255)
                                            cTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCtt                                                                                                                         <  4:108176/122‑1 (MQ=255)
                                            cTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATcag                                                                                                >  4:154505/1‑147 (MQ=255)
                                                cTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATcagcagcag                                                                                          <  5:219649/149‑1 (MQ=255)
                                                 tCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAgg                                                                                         >  2:216545/1‑149 (MQ=255)
                                                   aTGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCGCCAGGCTTGCAACTCCTCAACGTTCTCCATCAGCAGCAGGCg                                                                                       >  2:152452/1‑149 (MQ=255)
                                                   aTGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCg                                                                                       >  4:7235/1‑149 (MQ=255)
                                                    tGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGt                                                                                      <  6:338945/149‑1 (MQ=255)
                                                         ggATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAActc                                                                                                                  >  1:197358/1‑116 (MQ=255)
                                                         ggATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAActc                                                                                                                  <  2:197358/116‑1 (MQ=255)
                                                          gatgatCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCgg                                                                                >  6:377319/1‑149 (MQ=255)
                                                          gatgatCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCgg                                                                                <  7:354063/149‑1 (MQ=255)
                                                               tCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACACATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTgg                                                                            >  4:62590/1‑148 (MQ=255)
                                                                cGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGc                                                                           >  1:242827/1‑148 (MQ=255)
                                                                 ggcaggcaAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTCCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCa                                                                          >  7:182941/1‑148 (MQ=255)
                                                                  gcaggcaAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCaaa                                                                        >  5:342847/1‑149 (MQ=255)
                                                                  gcaggcaAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCATTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCaaa                                                                        >  2:296238/1‑149 (MQ=255)
                                                                  gcaggcaAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCATTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTGCTCCATCAGCAGCAGGCGTGGGCGGGTGGCaaa                                                                        >  2:301106/1‑149 (MQ=255)
                                                                  gcaggcaAGCCAGTGGTGAGCCAGTCGACAAATGCCTTGAAGGTGCGACCCATTTCAGACTTGATCGGATTGCCTTCGTCGTTCAGCTGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCaaa                                                                        <  5:26752/149‑1 (MQ=255)
                                                                  gcaggcaAGACAGTGGTGAGGAAGGCGACAAATGCCTTGAAGGTGCTACGCATTTCCGACTTGATCGGCTTGCCTTCGTAGTTCAGCGGCCCCCAGGATTGGAACTCATCAACGTTCTCCATCAGCAGCAGGAGTGGGCGGGTGGCcaa                                                                        >  6:26752/1‑149 (MQ=255)
                                                                        aaGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGGGGGTGGCAAACACCCa                                                                  <  5:238014/149‑1 (MQ=255)
                                                                              gtggtgAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAAc                                                            <  3:7235/149‑1 (MQ=255)
                                                                                    aggcaggcGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACaa                                                       <  1:136462/148‑1 (MQ=255)
                                                                                        aggcGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCATAACCCAt                                                  <  3:62590/149‑1 (MQ=255)
                                                                                        aggcGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCATTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCAt                                                  <  1:296238/149‑1 (MQ=255)
                                                                                        aggcGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCATTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCAt                                                  <  1:301106/149‑1 (MQ=255)
                                                                                           cGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCGCCATCAGCAGCAGGCGTGGGCGGGTGGCAACCAGCCAGTGACCCACAACCCATgg                                                >  8:316809/1‑148 (MQ=255)
                                                                                               aaaTGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGCGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCagc                                            >  6:71780/1‑148 (MQ=255)
                                                                                                      ttGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCAc                                    >  5:263162/1‑149 (MQ=255)
                                                                                                               cGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACt                           <  5:377319/149‑1 (MQ=255)
                                                                                                                gACCAATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTg                          >  5:173970/1‑149 (MQ=255)
                                                                                                                 aCCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTg                                                               >  6:137353/1‑111 (MQ=255)
                                                                                                                 aCCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTg                                                               <  5:137353/111‑1 (MQ=255)
                                                                                                                    cATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTg                      <  3:130841/149‑1 (MQ=255)
                                                                                                                           gACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGcgcc               <  6:263162/149‑1 (MQ=255)
                                                                                                                           gACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGcgcc               >  1:7977/1‑149 (MQ=255)
                                                                                                                                gATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCcct          >  7:151626/1‑149 (MQ=255)
                                                                                                                                   cGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCccttgg       <  5:162295/149‑1 (MQ=255)
                                                                                                                                     gCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACAccc                                                                   >  6:10141/1‑87 (MQ=255)
                                                                                                                                         gCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCccttggccttg  <  1:152452/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
ACCAATGCGTCCATCGGTACCCACTGCCCGATGCAGTCCATGACCTCGCTCATGTCGGGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTG  >  NC_002947/3416657‑3416941

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