Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,665,250 G→C 7.0% R206R (CGC→CGG)  PP_1459 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,665,2500GC 7.0% 178.2 / 3.2 57R206R (CGC→CGG) PP_1459CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base C (2/2);  total (28/29)
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.36e-01

CTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTCG  >  NC_002947/1665107‑1665397
                                                                                                                                               |                                                                                                                                                   
cTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGtt                                                                                                                                                 <  7:267170/148‑1 (MQ=255)
   tgCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGtt                                                                                                                                                 >  8:267170/1‑145 (MQ=255)
         gcTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACATGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGATGACGATCAGCTGATCGACGATGCGTTCGATCGGCtc                                                                                                                                       >  7:333062/1‑149 (MQ=255)
              cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCCTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGa                                                                                                                                  <  6:319190/149‑1 (MQ=255)
               aGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAg                                                                                                                                 >  4:43824/1‑149 (MQ=255)
                 ggTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgca                                                                                                                               <  6:154579/149‑3 (MQ=255)
                   tCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTg                                                                                                                             <  2:321786/149‑1 (MQ=255)
                   tCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAg                                                                                                                                 >  1:21419/1‑145 (MQ=255)
                   tCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAg                                                                                                                                 <  2:21419/145‑1 (MQ=255)
                       cTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGc                                                                                                                         >  5:312954/1‑149 (MQ=255)
                           gCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATc                                                                                                                     >  2:322564/1‑149 (MQ=255)
                            ccTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATcc                                                                                                                    <  5:80187/149‑1 (MQ=255)
                            ccTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATc                                                                                                                     >  8:176918/1‑148 (MQ=255)
                             cTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCTGCTCGTCGAGGTTGATGCCATTCa                                                                                                                   <  8:157314/149‑1 (MQ=255)
                             cTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCa                                                                                                                   >  5:253926/1‑149 (MQ=255)
                                 ggCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACtt                                                                                                                <  4:166132/148‑1 (MQ=255)
                                 ggCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTc                                                                                                               >  3:1019/1‑149 (MQ=255)
                                  gcagcaGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCg                                                                                                              <  1:337080/149‑1 (MQ=255)
                                                   aGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGAt                                                                                             >  5:235925/1‑149 (MQ=255)
                                                    gCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATg                                                                                            <  6:235925/149‑1 (MQ=255)
                                                         gtgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGtt                                                                                       <  7:176918/149‑1 (MQ=255)
                                                          tgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGtt                                                                                       >  1:134154/1‑148 (MQ=255)
                                                                ggATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGt                                                                                 <  6:253926/148‑1 (MQ=255)
                                                                   tCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATc                                                                             <  2:134154/149‑1 (MQ=255)
                                                                       ttCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGATGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTc                                                                                                               <  6:27288/111‑1 (MQ=255)
                                                                       ttCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTAGAAGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTc                                                                                                               >  5:27288/1‑111 (MQ=255)
                                                                         cGAACTGGTTGATTTCGTTGTGTTATACCGGCAGGCGGGTGTGGCTGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTc                                                                        <  8:333062/148‑1 (MQ=255)
                                                                           aCCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGAGTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAg                                                                     >  1:255439/1‑149 (MQ=255)
                                                                            ccTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                    >  1:267314/1‑149 (MQ=255)
                                                                                   gATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGtt                                                                                                                              <  1:132532/84‑1 (MQ=255)
                                                                                   gATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGtt                                                                                                                              >  2:132532/1‑84 (MQ=255)
                                                                                   gATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                    <  8:26579/142‑1 (MQ=255)
                                                                                   gATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                    >  7:26579/1‑142 (MQ=255)
                                                                                         gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCg                                                        <  3:43824/148‑1 (MQ=255)
                                                                                              ggTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGa                                                                                              <  1:184237/105‑1 (MQ=255)
                                                                                              ggTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGa                                                                                              >  2:184237/1‑105 (MQ=255)
                                                                                               gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAg                                                   >  6:154489/1‑147 (MQ=255)
                                                                                               gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGc                                                  <  6:312954/148‑1 (MQ=255)
                                                                                               gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCa                                                 <  1:322564/149‑1 (MQ=255)
                                                                                               gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTAGTCGAGGTTGATGACATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAg                                                   >  5:2559/1‑147 (MQ=255)
                                                                                                     cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGa                                                                               >  8:299721/1‑113 (MQ=255)
                                                                                                     cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGa                                                                               <  7:299721/113‑1 (MQ=255)
                                                                                                     cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCt                                           >  4:286577/1‑149 (MQ=255)
                                                                                                       gCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTgg                                         <  4:175237/149‑1 (MQ=255)
                                                                                                           gCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTgg                                         <  4:292784/145‑1 (MQ=255)
                                                                                                           gCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTgg                                         >  3:292784/1‑145 (MQ=255)
                                                                                                           gCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAg                                     >  8:220003/1‑149 (MQ=255)
                                                                                                             gggTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgc                                   >  7:261229/1‑149 (MQ=255)
                                                                                                               gtgtGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgcgc                                 >  8:150040/1‑149 (MQ=255)
                                                                                                                        cTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTg                        <  2:255439/149‑1 (MQ=255)
                                                                                                                          gACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGtt                      <  7:150040/149‑1 (MQ=255)
                                                                                                                           aCGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGt                       <  7:220003/147‑1 (MQ=255)
                                                                                                                                 aGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTc               <  8:261229/149‑1 (MQ=255)
                                                                                                                                  gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCa              >  4:177942/1‑149 (MQ=255)
                                                                                                                                  gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCa              >  8:106666/1‑149 (MQ=255)
                                                                                                                                      ctCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg           >  1:135321/1‑148 (MQ=255)
                                                                                                                                              tGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGtcgtcg  <  2:267314/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTCG  >  NC_002947/1665107‑1665397

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