Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,486,941 G→T 7.7% P357H (CCC→CAC)  PP_3093 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,486,9410GT 7.7% 149.1 / 3.4 52P357H (CCC→CAC) PP_3093hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base T (2/2);  total (30/22)
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.52e-01

GTTGTCCAGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGTTGTCGAGCATCACCCGGATC  >  NC_002947/3486793‑3487087
                                                                                                                                                     |                                                                                                                                                  
gTTGTCCAGGCGCACCGGACTCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCgg                                                                                                                                                    <  5:202102/149‑1 (MQ=255)
gTTGTCCAGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCgg                                                                                                                                                    >  2:32705/1‑149 (MQ=255)
    tCCAGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCAt                                                                                                                                                >  2:356482/1‑149 (MQ=255)
     ccAGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATg                                                                                                                                               >  4:99893/1‑149 (MQ=255)
     ccAGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGCAT‑GCCGCCGTGCacg                                                                                                                                               >  2:25166/1‑147 (MQ=255)
      cAGGCGCACCGGAATCGCGTTCGGCAGCACCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCATCCGGCGAACCGACCTTGAAGGCACGGCGCACCAGCTCAAGCAGGCT‑GGCGCAGGGCATGt                                                                                                                                              >  2:50213/1‑149 (MQ=255)
       aGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGTGCacgcc                                                                                                                                             >  2:231139/1‑145 (MQ=255)
           gcACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGCCATGGCGCTGTTCACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGCGGCACCAGCTCCAGCAGGCA‑GGCGCCGCGCATGTCGGcc                                                                                                                                         >  6:82582/1‑148 (MQ=255)
              ccGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTg                                                                                                                                       >  1:138121/1‑148 (MQ=255)
                    tCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGTCACCCGGCATGGCGGTGGTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GACGCCGTGCACGCCGACATGCCCGGCg                                                                                                                                <  1:25166/149‑1 (MQ=255)
                     cgcgTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGa                                                                                                                               >  6:86299/1‑149 (MQ=255)
                      gcgTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCTAGGCGCCGTGCACGCCGACATGCCCGGCg                                                                                                                                <  1:231139/148‑1 (MQ=255)
                            ggCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTGCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGaa                                                                                                                        <  8:369886/149‑1 (MQ=255)
                            ggCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTGCGGCGAACAGACCTTGAAGGCACGGGGCAGCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGTCGATATAGa                                                                                                                         <  5:412568/148‑1 (MQ=255)
                            ggCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGaa                                                                                                                        <  5:86299/149‑1 (MQ=255)
                            ggCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGaa                                                                                                                        <  1:393693/149‑1 (MQ=255)
                              cAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGt                                                                                                                      >  2:319303/1‑149 (MQ=255)
                                  cgcgcgGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGc                                                                                                                  >  8:251202/1‑149 (MQ=255)
                                  cgcgcgGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGc                                                                                                                  >  6:25569/1‑149 (MQ=255)
                                      cgGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCg                                                                                                              <  1:50213/149‑1 (MQ=255)
                                         cATGATGGAGGGGGACACCCGGCATGGCGCTGATGACCACCACTTCAATGTCTTCCGCCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGcc                                                                                                           <  6:417901/149‑1 (MQ=255)
                                              ttGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGgcg                                                                                                      >  1:36311/1‑149 (MQ=255)
                                                gAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGcgcg                                                                                                    >  2:393129/1‑149 (MQ=255)
                                                       acacCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTg                                                                                             <  1:32705/149‑1 (MQ=255)
                                                           ccGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGt                                                                                         >  5:34399/1‑149 (MQ=255)
                                                                    cgcTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGc                                                                                <  2:138121/149‑1 (MQ=255)
                                                                           gACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGc                                                                          >  1:10514/1‑148 (MQ=255)
                                                                              caccacTTCAATGTCTTCCGGCGAACCGACGTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATg                                                                      >  1:34901/1‑149 (MQ=255)
                                                                              caccacTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATg                                                                      >  7:336839/1‑149 (MQ=255)
                                                                                ccacTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGc                                                                    >  8:176681/1‑149 (MQ=255)
                                                                                  acTTCTATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCtt                                                                  <  5:82582/149‑1 (MQ=255)
                                                                                   cTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTg                                                                 <  6:34399/149‑1 (MQ=255)
                                                                                       aaTGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCt                                                             >  1:261949/1‑149 (MQ=255)
                                                                                         tGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTgg                                                           <  5:25569/149‑1 (MQ=255)
                                                                                         tGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTgg                                                           <  1:393129/149‑1 (MQ=255)
                                                                                             ttCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGt                                                         >  2:355015/1‑147 (MQ=255)
                                                                                                  gCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGtt                                                   >  4:159401/1‑148 (MQ=255)
                                                                                                    gAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTtgat                                                <  3:99893/149‑1 (MQ=255)
                                                                                                             ttGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCa                                       >  3:138360/1‑149 (MQ=255)
                                                                                                               gAAGGCACGGGGCACCAGCTCCTGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCTATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCAc                                      <  7:251202/148‑1 (MQ=255)
                                                                                                                  ggCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGt                                  <  2:30727/149‑1 (MQ=255)
                                                                                                                             ccAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCAGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCacgacg                       >  2:308274/1‑149 (MQ=255)
                                                                                                                              cAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAgg                                                     >  4:188586/1‑118 (MQ=255)
                                                                                                                              cAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAgg                                                     <  3:188586/118‑1 (MQ=255)
                                                                                                                               aGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGtt                     >  8:359530/1‑149 (MQ=255)
                                                                                                                               aGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGtt                     <  2:261949/149‑1 (MQ=255)
                                                                                                                                 cTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGGTGATGACAATGCACTGGTTGGGCACGACGTGGt                   >  6:65419/1‑149 (MQ=255)
                                                                                                                                  tCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGTGATGGCTTGGTCTGGGTCAGGTTGATGACAAt                                          <  8:80342/125‑1 (MQ=255)
                                                                                                                                  tCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGTGATGGCTTGGTCTGGGTCAGGTTGATGACAAt                                          >  7:80342/1‑125 (MQ=255)
                                                                                                                                  tCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGTTGTc                  >  5:296258/1‑149 (MQ=255)
                                                                                                                                   ccAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGTTGTCg                 >  6:172213/1‑149 (MQ=255)
                                                                                                                                                 gcCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGTTGTCGAGCATCACCCGGa    <  2:10514/149‑1 (MQ=255)
                                                                                                                                                   cGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGTTGTCGAGCATCACCCGGATc  <  8:336839/149‑1 (MQ=255)
                                                                                                                                                   cGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCAAGTGCTATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGAAGTTGTCGAGCATCACCCGGAt   >  6:227808/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTTGTCCAGGCGCACCGGAATCGCGTTCGGCAGCCGCGCGGCATGATTGAGGGTGACACCCGGCATGGCGCTGTTGACCACCACTTCAATGTCTTCCGGCGAACCGACCTTGAAGGCACGGGGCACCAGCTCCAGCAGGCT‑GGCGCCGGGCATGTCGGCGTGCACGGCGATATAGAAGTCGGCTTCGGCCAGGCGCGGGTCGTGCAGTTGGCCTTGCCAGTGCGATGGCTTGGTCTGGGTCAGGTTGATGACAATGCACTGGTTGGGCACGACGTTGTCGAGCATCACCCGGATC  >  NC_002947/3486793‑3487087

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