Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,667,937 C→T 7.4% noncoding (70/97 nt) PP_mr59 → SAH_riboswitch regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,667,9370CT 7.4% 148.1 / 2.0 54noncoding (70/97 nt)PP_mr59SAH_riboswitch regulatory region
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base T (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

AACCTTGGGCAAGCCGGCGGGGGCAGCTATAATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCC  >  NC_002947/5667794‑5668082
                                                                                                                                               |                                                                                                                                                 
aaCCTTGGGCAAGCCGGCGGGGGCAGCTATAATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAgcg                                                                                                                                              <  6:176775/149‑1 (MQ=255)
    ttGGGCAAGCCGGCGGGGGCAGCTATAATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCtt                                                                                                                                           >  6:302928/1‑148 (MQ=255)
                          cTATAATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGc                                                                                                                    <  5:302928/149‑1 (MQ=255)
                           tataATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGca                                                                                                                   >  5:72491/1‑149 (MQ=255)
                            ataATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCAAGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTACCGCGAGCGCTTAACGCACAACGGCCCCCATTCGca                                                                                                                   >  2:166430/1‑148 (MQ=255)
                                    cgcCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACg                                                                                                            <  3:151540/147‑1 (MQ=255)
                                    cgcCGTTTTACATGGGTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTAc                                                                                                             <  5:73562/146‑1 (MQ=255)
                                     gcCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCTCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAAt                                                                                                         >  5:59277/1‑149 (MQ=255)
                                        gTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGGTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGa                                                                                                      <  1:284172/149‑1 (MQ=255)
                                         ttttACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGTGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAg                                                                                                     <  6:364912/149‑1 (MQ=255)
                                            tACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGct                                                                                                  >  6:271141/1‑149 (MQ=255)
                                             aCATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGct                                                                                                  <  7:169690/148‑1 (MQ=255)
                                             aCATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTACCGCGAGCGCTTAACGCACAACGGCGCCCATTCGGACATTACGAATGGAGGctc                                                                                                 >  1:278685/1‑149 (MQ=255)
                                                 ggTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCAt                                                                                             >  2:193671/1‑149 (MQ=255)
                                                       ccACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAgcgc                                                                                       >  8:284175/1‑149 (MQ=255)
                                                        cacaGCGGTGCCATGCCCGCGACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGTCGCCAATTCGGACATTACGAATGGAGGCTCTCATGAGCGCt                                                                                      >  2:236744/1‑149 (MQ=255)
                                                         acaGAGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCAGATGGGGCGTTGCTAGATCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACAGTACGAATGGCGGCTCTCATGCGCGCTg                                                                                     >  7:233370/1‑149 (MQ=255)
                                                             cGGTCCCATGCCCGCCACCTGTCCGAGGTGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCaaa                                                                                 <  1:236744/149‑1 (MQ=255)
                                                                  ccATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCTCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGc                                                                            <  6:59277/149‑1 (MQ=255)
                                                                  ccATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGc                                                                            <  1:193671/149‑1 (MQ=255)
                                                                   cATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGcc                                                                           >  8:209861/1‑149 (MQ=255)
                                                                    aTGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAGCATGCct                                                                          >  3:99167/1‑149 (MQ=255)
                                                                    aTGCCCGCCACCTGTCCGAGGGGAGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCct                                                                          >  5:139332/1‑149 (MQ=255)
                                                                        ccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGTCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctgctg                                                                      >  2:305795/1‑149 (MQ=255)
                                                                        ccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCt                                                                                      >  4:235060/1‑133 (MQ=255)
                                                                        ccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCt                                                                                      <  3:235060/133‑1 (MQ=255)
                                                                        ccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctgctg                                                                      >  5:353753/1‑149 (MQ=255)
                                                                               cTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAATCATGCCTGCTGGTTTTAc                                                               <  4:99167/149‑1 (MQ=255)
                                                                                         ggCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACaaa                                                     <  7:284175/149‑1 (MQ=255)
                                                                                            gcTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTc                                                  <  2:14993/149‑1 (MQ=255)
                                                                                            gcTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTc                                                  >  8:276388/1‑149 (MQ=255)
                                                                                            gcTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTc                                                  >  5:323477/1‑149 (MQ=255)
                                                                                                cagcagGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGcc                                               <  6:353753/148‑1 (MQ=255)
                                                                                                cagcagGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGcc                                               <  5:271141/148‑1 (MQ=255)
                                                                                                cagcagGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCg                                              <  6:139332/149‑1 (MQ=255)
                                                                                                 agcagGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGTTTTTACCGATTACAAAGTCGCCGa                                             <  2:278685/149‑1 (MQ=255)
                                                                                                 agcagGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGa                                             <  7:276388/149‑1 (MQ=255)
                                                                                                           gCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTg                                   >  6:85680/1‑149 (MQ=255)
                                                                                                            gcTGTCAGGCTCGGTTGGGGCGTTGCTCGCTCCCGCGAGCGATTAACTCTCAACGGCTCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACGGATTACAAAGTCGCCGACATCTCCCTgg                                  <  1:166430/148‑1 (MQ=255)
                                                                                                             cTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCATATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGgc                                 >  3:90568/1‑149 (MQ=255)
                                                                                                              tGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGCGCGCGTAACGCACAACGGCGCCCATTCGCCCATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTAAAAAGTCGCCGACATCTCCCTGgcc                                >  3:3633/1‑149 (MQ=255)
                                                                                                                tCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGgcc                                <  6:72491/147‑1 (MQ=255)
                                                                                                                    gCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTgg                          <  7:209861/149‑1 (MQ=255)
                                                                                                                     cTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCATATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTggg                         <  4:90568/149‑1 (MQ=255)
                                                                                                                              ggCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTcgc                 >  7:114769/1‑148 (MQ=255)
                                                                                                                                         ctcCCGTGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTg                           >  3:269226/1‑127 (MQ=255)
                                                                                                                                         ctcCCGTGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTg                           >  3:264496/1‑127 (MQ=255)
                                                                                                                                         ctcCCGTGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTg                           <  4:264496/127‑1 (MQ=255)
                                                                                                                                         ctcCCGTGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTg                           <  4:269226/127‑1 (MQ=255)
                                                                                                                                         ctcCCGCGAGCGCCTAACGCACAACCGCGCCCATTCGCACATTACGAATGGAGGCCCTCATGAGCGCGGCACACATGCCTGCGGGTTTCACCGCGTACAACGTCCCCGACATCTCCCTCGCCGCCTGGGGCCGCCGCGAAACCATCAtt     <  8:233370/149‑1 (MQ=255)
                                                                                                                                            ccGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCt                            >  2:116751/1‑123 (MQ=255)
                                                                                                                                            ccGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGcc  >  7:198660/1‑149 (MQ=255)
                                                                                                                                            ccGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGcc  >  2:326357/1‑149 (MQ=255)
                                                                                                                                            ccGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAAAGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCt                            <  1:116751/123‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
AACCTTGGGCAAGCCGGCGGGGGCAGCTATAATCTGCGCCGTTTTACATGGTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCC  >  NC_002947/5667794‑5668082

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