Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,601,674 C→T 7.5% R131C (CGC→TGC)  PP_3177 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,601,6740CT 7.5% 152.9 / 2.5 53R131C (CGC→TGC) PP_3177hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base T (2/2);  total (24/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 = 7.79e-01

GATCGTCGAACTGGATACCGACGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAGCCCGTCAGTTGGGC  >  NC_002947/3601526‑3601820
                                                                                                                                                    |                                                                                                                                                  
gATCGTCGAACTGGATACCGACGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTgc                                                                                                                                                    <  8:176560/149‑1 (MQ=255)
    gtcgAACTGGATACCGACGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCAt                                                                                                                                                >  6:277408/1‑149 (MQ=255)
             gATACCGACGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGAccgccg                                                                                                                                        >  1:310798/1‑148 (MQ=255)
             gATACCGACGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGAccgccg                                                                                                                                        >  1:53535/1‑148 (MQ=255)
                    aCGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAAtgctg                                                                                                                                <  2:310798/149‑1 (MQ=255)
                     cgcgGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAAtgctgc                                                                                                                               <  8:339935/149‑1 (MQ=255)
                        cgCCATGTGCCTGAAGGTCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCTATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGTCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTgg                                                                                                                            <  4:285032/148‑1 (MQ=255)
                         gCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGc                                                                                                                           <  1:137332/149‑1 (MQ=255)
                              gTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAAc                                                                                                                      <  2:53535/149‑1 (MQ=255)
                                    gAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCAt                                                                                                                                                >  4:315292/1‑117 (MQ=255)
                                    gAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCAt                                                                                                                                                <  3:315292/117‑1 (MQ=255)
                                       gggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGTGCATGACCGCCGAAATGCTGCTGGCCCAACGCAAtt                                                                                                                <  6:8801/146‑1 (MQ=255)
                                       gggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGTGCATGACCGCCGAAATGCTGCTGGCCCAACGCAAtt                                                                                                                <  6:184242/146‑1 (MQ=255)
                                       gggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGTGCATGACCGCCGAAATGCTGCTGGCCCAACGCAAtt                                                                                                                >  5:8801/1‑146 (MQ=255)
                                       gggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGTGCATGACCGCCGAAATGCTGCTGGCCCAACGCAAtt                                                                                                                >  5:184242/1‑146 (MQ=255)
                                       gggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTa                                                                                                              >  6:45054/1‑148 (MQ=255)
                                        ggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCACAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGAccgc                                                                                                                                          <  6:70332/119‑1 (MQ=255)
                                        ggCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCACAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGAccgc                                                                                                                                          >  5:70332/1‑119 (MQ=255)
                                          cGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAgg                                                                                                          <  2:176654/149‑1 (MQ=255)
                                                 cccGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCa                                                                                                   <  1:42673/149‑1 (MQ=255)
                                                  ccGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAg                                                                                                  <  5:277408/149‑1 (MQ=255)
                                                      cGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAAt                                                                                              <  3:257334/149‑1 (MQ=255)
                                                       gTGAACCTGCCATTGATGCGCGATGGGCGGCTGATGGGCGTGTTCGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATg                                                                                             <  4:30419/149‑1 (MQ=255)
                                                           aCCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCg                                                                                          >  4:274656/1‑148 (MQ=255)
                                                              tGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAgg                                                                                      >  4:50511/1‑149 (MQ=255)
                                                               gCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGaa                                                                                        <  6:203851/146‑1 (MQ=255)
                                                               gCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGaa                                                                                        >  5:203851/1‑146 (MQ=255)
                                                                ccATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCa                                                                                    <  5:139553/149‑1 (MQ=255)
                                                                  aTTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAgc                                                                                  >  1:141427/1‑149 (MQ=255)
                                                                      aTGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCtg                                                                              >  4:339028/1‑149 (MQ=255)
                                                                       tGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCtgt                                                                             <  2:236785/149‑1 (MQ=255)
                                                                         cTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTga                                                                           <  5:63977/149‑1 (MQ=255)
                                                                            gATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTgacga                                                                        >  6:129599/1‑149 (MQ=255)
                                                                                 gcggcTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGAc                                                                                                     >  8:91887/1‑115 (MQ=255)
                                                                                 gcggcTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGAc                                                                                                     <  7:91887/115‑1 (MQ=255)
                                                                                    gcTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTgacgac                                                                       <  8:169731/142‑1 (MQ=255)
                                                                                    gcTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTgacgac                                                                       >  7:169731/1‑142 (MQ=255)
                                                                                    gcTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTgg                                                                <  6:149675/149‑1 (MQ=255)
                                                                                       gATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGgcgc                                                             >  7:298540/1‑149 (MQ=255)
                                                                                        aTGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGgcgc                                                             <  5:129599/148‑1 (MQ=255)
                                                                                          gggCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTAc                                                          >  7:201893/1‑149 (MQ=255)
                                                                                                            tGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGt                                         >  6:92821/1‑148 (MQ=255)
                                                                                                              gCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTGAc                                      >  8:190753/1‑149 (MQ=255)
                                                                                                               cGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACt                                     <  5:45054/149‑1 (MQ=255)
                                                                                                                     gcagcaGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTcgccgc                               >  1:333306/1‑149 (MQ=255)
                                                                                                                         cagGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCtt                           >  4:29744/1‑149 (MQ=255)
                                                                                                                         cagGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCtt                           <  8:201893/149‑1 (MQ=255)
                                                                                                                         cagGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCtt                           <  3:274656/149‑1 (MQ=255)
                                                                                                                           ggTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGc                         <  3:339028/149‑1 (MQ=255)
                                                                                                                                    cTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAg                >  4:114823/1‑149 (MQ=255)
                                                                                                                                      aTGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAGcc              >  6:41645/1‑149 (MQ=255)
                                                                                                                                         gCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAGCCCGt           <  8:298540/149‑1 (MQ=255)
                                                                                                                                                 gTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAGCCCGTCAGTTgg    <  5:41645/148‑1 (MQ=255)
                                                                                                                                                  tGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAGCCCGTCAGTTGGGc  <  5:92821/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATCGTCGAACTGGATACCGACGCGGCCAAGTGCCTGAAGGGCGTGCAGCCCGGCGTGAACCTGCCATTGATGCTCGATGGGCGGCTGATGGGCGTGTTGGGCCTGACTGGCGACCCGCAGCAGGTGCGCACCTATGGCGAGCTGGTGCGCATGACCGCCGAAATGCTGCTGGCCCAACGCAATTTACAGGCTGACCAGCAATGGCGAAGGCAGCGCTGTGACGACCTGCTGGCGCTACTGCTCGGCGGCAGCGGTAACTCGCCGCGCTTGCTGGACGAAGCCCGTCAGTTGGGC  >  NC_002947/3601526‑3601820

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