Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 432,016 A→G 7.8% C100R (TGC→CGC) ‡ PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947432,0160AG 7.8% 143.4 / 2.6 51C100R (TGC→CGC) ‡PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (26/21);  new base G (2/2);  total (28/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 = 9.20e-01

GGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGAAG  >  NC_002947/431875‑432164
                                                                                                                                              |                                                                                                                                                    
gggCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCg                                                                                                                                               >  8:95248/1‑149 (MQ=255)
 ggCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGc                                                                                                                                              >  2:199928/1‑149 (MQ=255)
 ggCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGc                                                                                                                                              >  6:304925/1‑149 (MQ=255)
   cAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCa                                                                                                                                                      <  6:15950/139‑1 (MQ=255)
   cAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCa                                                                                                                                                      >  5:15950/1‑139 (MQ=255)
   cAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTg                                                                                                                                            >  4:230645/1‑149 (MQ=255)
              aaCACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGaa                                                                                                                                 <  3:21277/149‑1 (MQ=255)
                  caTCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCTGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTg                                                                                                                             <  7:45546/149‑1 (MQ=255)
                         cATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGa                                                                                                                      <  7:187638/149‑1 (MQ=255)
                                       ggATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGc                                                                                                         >  2:252548/1‑148 (MQ=255)
                                            aaTGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCTAGGTg                                                                                                   <  7:141905/149‑1 (MQ=255)
                                                  gCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcg                                                                                             >  8:173382/1‑149 (MQ=255)
                                                           aGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGTGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCaa                                                                                     >  1:72164/1‑149 (MQ=255)
                                                             ttGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCa                                                                                  >  1:137504/1‑149 (MQ=255)
                                                                caggcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCAcgc                                                                               >  1:150163/1‑149 (MQ=255)
                                                                   gcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGc                                                                            <  3:40161/149‑1 (MQ=255)
                                                                          gtgtCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCtt                                                                     >  4:190521/1‑149 (MQ=255)
                                                                                 ccccGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTg                                                                                                           <  5:266623/104‑1 (MQ=255)
                                                                                 ccccGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTg                                                                                                           >  6:266623/1‑104 (MQ=255)
                                                                                   ccGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATg                                                              >  3:252599/1‑147 (MQ=255)
                                                                                   ccGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTa                                                            <  1:199928/149‑1 (MQ=255)
                                                                                      gCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGaa                                                         <  3:230645/149‑1 (MQ=255)
                                                                                      gCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGaa                                                         <  5:304925/149‑1 (MQ=255)
                                                                                       ccAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAg                                                        >  2:155291/1‑149 (MQ=255)
                                                                                            tGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTc                                                   <  7:95248/149‑1 (MQ=255)
                                                                                                 ggTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCg                                              <  2:72164/149‑1 (MQ=255)
                                                                                                            gCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCgg                                     >  8:13376/1‑148 (MQ=255)
                                                                                                              ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                                  >  4:303118/1‑149 (MQ=255)
                                                                                                              ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                                  >  2:37054/1‑149 (MQ=255)
                                                                                                              ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                                  <  7:173382/149‑1 (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                           >  6:177599/1‑148 (MQ=255)
                                                                                                                      ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                           >  5:156617/1‑148 (MQ=255)
                                                                                                                      ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                           >  5:114664/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                       >  5:308251/2‑141 (MQ=255)
                                                                                                                                 accAGCGGGGAGTGCTTCCCGCTGGCGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATg                       <  6:308251/140‑1 (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        >  7:96520/1‑149 (MQ=255)
                                                                                                                                        ggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCg        >  3:106403/1‑149 (MQ=255)
                                                                                                                                             cACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGaa   <  3:70072/149‑1 (MQ=255)
                                                                                                                                              aCTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGAAg  >  7:100825/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGTAGAAG  >  NC_002947/431875‑432164

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: