Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,377,197 C→T 7.5% R125R (CGG→CGA)  PP_3852 ← BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,377,1970CT 7.5% 152.0 / 2.0 53R125R (CGG→CGA) PP_3852BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

CGCGCCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCGGGAATAA  >  NC_002947/4377055‑4377341
                                                                                                                                              |                                                                                                                                                
cgcgcCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGc                                                                                                                                               <  5:213595/146‑1 (MQ=255)
cgcgcCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGa                                                                                                                                            >  6:1952/1‑149 (MQ=255)
cgcgcCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGa                                                                                                                                            <  8:167314/149‑1 (MQ=255)
cgcgcCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGa                                                                                                                                            >  1:71291/1‑149 (MQ=255)
 gcgcCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGAt                                                                                                                                           <  2:278051/149‑1 (MQ=255)
     cATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGaa                                                                                                                                       <  5:1952/149‑1 (MQ=255)
         cGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCt                                                                                                                                    <  6:221654/148‑1 (MQ=255)
          ggTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTg                                                                                                                                   >  4:308583/1‑148 (MQ=255)
          ggTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTg                                                                                                                                   >  8:119595/1‑148 (MQ=255)
           gTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTgg                                                                                                                                  >  7:185548/1‑148 (MQ=255)
                       cTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGa                                                                                                                     <  1:275569/149‑1 (MQ=255)
                                   cTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGaaa                                                                                                                               >  5:300531/1‑127 (MQ=255)
                                   cTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGaaa                                                                                                                               <  6:300531/127‑1 (MQ=255)
                                      gCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGCCTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTg                                                                                                      >  4:138082/1‑149 (MQ=255)
                                       cGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTgg                                                                                                     <  7:119595/149‑1 (MQ=255)
                                       cGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTgg                                                                                                     <  8:230579/149‑1 (MQ=255)
                                       cGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTgg                                                                                                     <  5:164378/149‑1 (MQ=255)
                                              ggccgtggccTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTgg                                                                                              >  5:373877/1‑149 (MQ=255)
                                              ggccgtggccTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTgg                                                                                              >  8:229022/1‑149 (MQ=255)
                                              ggccgtggccTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGACGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTg                                                                                               >  7:39102/1‑148 (MQ=255)
                                                ccgtggccTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGct                                                                                            >  1:41240/1‑149 (MQ=255)
                                                     gccTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTAt                                                                                        <  6:373877/148‑1 (MQ=255)
                                                       cTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGc                                                                                     <  8:185548/149‑1 (MQ=255)
                                                         ggTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCt                                                                                   >  3:220660/1‑149 (MQ=255)
                                                         ggTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCt                                                                                   <  4:240426/149‑1 (MQ=255)
                                                              gCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCgg                                                                              >  5:397133/1‑149 (MQ=255)
                                                               cACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGTTATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGgc                                                                             <  3:308583/149‑1 (MQ=255)
                                                               cACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGgc                                                                             >  8:302616/1‑149 (MQ=255)
                                                                    tttGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAg                                                                        >  6:202587/1‑149 (MQ=255)
                                                                    tttGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAg                                                                        >  1:209746/1‑149 (MQ=255)
                                                                        cTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGct                                                                                            >  8:360188/1‑125 (MQ=255)
                                                                        cTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGct                                                                                            <  7:360188/125‑1 (MQ=255)
                                                                        cTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCgg                                                                     <  5:202587/148‑1 (MQ=255)
                                                                          aGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGt                                                                   >  7:260429/1‑148 (MQ=255)
                                                                             cGGTAGAAAGGGTCGTCCGTGTTTCGCTCGAGGCGCGCCGGGTAAGCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTCCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCAtt                                                               <  7:67790/149‑1 (MQ=255)
                                                                                    aaGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCt                                                        <  2:41240/149‑1 (MQ=255)
                                                                                            tcCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACaa                                                >  1:223526/1‑149 (MQ=255)
                                                                                                tgtTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATaca                                            <  4:220660/149‑3 (MQ=255)
                                                                                                      gCTCGATGCGCGCCGGGGAATCCCCATAGACGCGTGCCGGTCGCTGAACGAAGCTGGAAACGATATTGGGCTCAGTAGTGCCTTGGCGAATGGCTCTATCGCCACCCGGCGCAAGCGGGTGATTTTTGGCTAAATACAATAGAgcgcag                                      >  2:98351/1‑147 (MQ=255)
                                                                                                      gCTCGATGCGCGCCCGGGAATCCACATAGACGCGTGCCGGTCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCCGGGAATTTTTCGCTAAATACAATATAGCGCTc                                      >  4:266553/1‑149 (MQ=255)
                                                                                                              cgcgcCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGc                              <  7:302616/149‑1 (MQ=255)
                                                                                                               gcgcCGGGGAATCCACATAGACGCGTGCCGGCCGCTGAGCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCa                             <  3:138082/149‑1 (MQ=255)
                                                                                                                 gcCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAgg                           >  3:311494/1‑149 (MQ=255)
                                                                                                                   cGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGc                          >  4:104175/1‑148 (MQ=255)
                                                                                                                          tCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTc                  <  6:397133/149‑1 (MQ=255)
                                                                                                                           ccACATAGACGCGTGCCGGCCGCTGATAGAAGCTGGAAACGATATTCGACTCTGTAGTGCCATGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGACATATTTCGCTAAATACAATAGAGCGCACCCGTTGGCAGGCGAGCTGTc                  >  1:316472/1‑148 (MQ=255)
                                                                                                                            cacaTAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGa                >  6:380779/1‑149 (MQ=255)
                                                                                                                                tAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTc            >  2:368701/1‑149 (MQ=255)
                                                                                                                                   aCGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATagag                                           <  8:354515/115‑1 (MQ=255)
                                                                                                                                    cgggTGCCGGTCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCgg        <  3:266553/146‑1 (MQ=255)
                                                                                                                                    cgcgTGCCGGTCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCgg        <  1:98351/149‑1 (MQ=255)
                                                                                                                                    cgcgTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCgg        <  3:104175/149‑1 (MQ=255)
                                                                                                                                     gcgTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCgg        >  2:24706/1‑148 (MQ=255)
                                                                                                                                          ccggccgCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCGGGAATaa  <  7:229022/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGCGCCATGCGGTCCAGTATCGGCTGGGACCGGTTCTGGCGATTGTGGCCGTGGCCTGGTACGCACGATTTGCTAGGCGGTAGAAAGGGTCGTCCGTGTTTCGCTCGATGCGCGCCGGGGAATCCACATAGACGCGTGCCGGCCGCTGATCGAAGCTGGAAACGATATTGGACTCTGTAGTGCCTTGGCGATTGGCTCTATCGCCTCCCGGCGCAAGCGGGTCATTTTTCGCTAAATACAATAGAGCGCTCCCGTTGGCAGGCGAGCTGTCGACGTCGCGGGAATAA  >  NC_002947/4377055‑4377341

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