Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 432,015 C→T 7.7% C100Y (TGC→TAC) ‡ PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947432,0150CT 7.7% 150.7 / 2.1 52C100Y (TGC→TAC) ‡PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base T (2/2);  total (29/23)
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

ACCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGAA  >  NC_002947/431867‑432163
                                                                                                                                                     |                                                                                                                                                    
acCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGc                                                                                                                                                      >  8:141905/1‑149 (MQ=255)
        gggCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCg                                                                                                                                              >  8:95248/1‑149 (MQ=255)
         ggCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGc                                                                                                                                             >  2:199928/1‑149 (MQ=255)
         ggCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGc                                                                                                                                             >  6:304925/1‑149 (MQ=255)
           cAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCa                                                                                                                                                     <  6:15950/139‑1 (MQ=255)
           cAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCa                                                                                                                                                     >  5:15950/1‑139 (MQ=255)
           cAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTg                                                                                                                                           >  4:230645/1‑149 (MQ=255)
                      aaCACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGaa                                                                                                                                <  3:21277/149‑1 (MQ=255)
                          caTCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCTGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTg                                                                                                                            <  7:45546/149‑1 (MQ=255)
                                 cATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGa                                                                                                                     <  7:187638/149‑1 (MQ=255)
                                               ggATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGc                                                                                                        >  2:252548/1‑148 (MQ=255)
                                                    aaTGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCTAGGTg                                                                                                  <  7:141905/149‑1 (MQ=255)
                                                          gCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcg                                                                                            >  8:173382/1‑149 (MQ=255)
                                                                   aGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGTGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCaa                                                                                    >  1:72164/1‑149 (MQ=255)
                                                                     ttGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCa                                                                                 >  1:137504/1‑149 (MQ=255)
                                                                        caggcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCAcgc                                                                              >  1:150163/1‑149 (MQ=255)
                                                                           gcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGc                                                                           <  3:40161/149‑1 (MQ=255)
                                                                                  gtgtCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCtt                                                                    >  4:190521/1‑149 (MQ=255)
                                                                                         ccccGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTg                                                                                                          >  6:266623/1‑104 (MQ=255)
                                                                                         ccccGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTg                                                                                                          <  5:266623/104‑1 (MQ=255)
                                                                                           ccGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATg                                                             >  3:252599/1‑147 (MQ=255)
                                                                                           ccGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTa                                                           <  1:199928/149‑1 (MQ=255)
                                                                                              gCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGaa                                                        <  3:230645/149‑1 (MQ=255)
                                                                                              gCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGaa                                                        <  5:304925/149‑1 (MQ=255)
                                                                                               ccAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAg                                                       >  2:155291/1‑149 (MQ=255)
                                                                                                    tGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTc                                                  <  7:95248/149‑1 (MQ=255)
                                                                                                         ggTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCg                                             <  2:72164/149‑1 (MQ=255)
                                                                                                                    gCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCgg                                    >  8:13376/1‑148 (MQ=255)
                                                                                                                      ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                                 >  2:37054/1‑149 (MQ=255)
                                                                                                                      ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                                 <  7:173382/149‑1 (MQ=255)
                                                                                                                      ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                                 >  4:303118/1‑149 (MQ=255)
                                                                                                                       cGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAgg                                <  2:150163/149‑1 (MQ=255)
                                                                                                                            aCCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGttt                             <  7:13376/147‑1 (MQ=255)
                                                                                                                            aCCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGa                           >  6:166052/1‑149 (MQ=255)
                                                                                                                             cccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                          <  1:155291/149‑1 (MQ=255)
                                                                                                                              ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGccc                                         <  8:91531/133‑1 (MQ=255)
                                                                                                                              ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGccc                                         >  7:91531/1‑133 (MQ=255)
                                                                                                                              ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                          >  5:114664/1‑148 (MQ=255)
                                                                                                                              ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                          >  5:156617/1‑148 (MQ=255)
                                                                                                                              ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                          >  6:177599/1‑148 (MQ=255)
                                                                                                                               cagcagGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGa                        >  3:182360/1‑149 (MQ=255)
                                                                                                                                agcagGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGAt                       >  4:70072/1‑149 (MQ=255)
                                                                                                                                       cTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCa                >  2:224581/1‑149 (MQ=255)
                                                                                                                                       cTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCa                <  3:303118/149‑1 (MQ=255)
                                                                                                                                         accAGCGGGGAGTGCTTCCCGCTGGCGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATg                      <  6:308251/140‑1 (MQ=255)
                                                                                                                                         accAGCGGGGAGTGCTTCCCGCTGGCGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATg                      >  5:308251/2‑141 (MQ=255)
                                                                                                                                          ccAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCAc               <  5:166052/147‑1 (MQ=255)
                                                                                                                                           cAGCGGGGAGTGCTTCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCg             >  2:277558/1‑148 (MQ=255)
                                                                                                                                            aGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCa           <  6:114664/149‑1 (MQ=255)
                                                                                                                                             gCGGGGAGTGCTTCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGc                                                    <  6:45627/107‑1 (MQ=255)
                                                                                                                                               gggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGc        <  1:37054/149‑1 (MQ=255)
                                                                                                                                                ggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCg       >  3:106403/1‑149 (MQ=255)
                                                                                                                                                ggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCg       >  7:96520/1‑149 (MQ=255)
                                                                                                                                                     cACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGaa  <  3:70072/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCG‑GGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGAA  >  NC_002947/431867‑432163

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