Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,797,339 A→G 7.7% D131G (GAT→GGT)  PP_2449 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,797,3390AG 7.7% 145.6 / 2.3 52D131G (GAT→GGT) PP_2449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base G (2/2);  total (29/23)
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.97e-01

GTCGACGATACCGGCAACAACCTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACC  >  NC_002947/2797191‑2797481
                                                                                                                                                    |                                                                                                                                              
gtcgACGATACCGGCAACAACCTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGa                                                                                                                                                >  8:321784/1‑149 (MQ=255)
gtcgACGATACCGGCAACAACCTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGa                                                                                                                                                <  8:124561/149‑1 (MQ=255)
   gacgaTACCGGCAACAACCTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGatat                                                                                                                                             >  5:169578/1‑149 (MQ=255)
   gacgaTACCGGCAACAACCTGCGTATCCGGGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGat                                                                                                                                               <  3:160413/147‑1 (MQ=255)
       aTACCGGCAACAACCTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACg                                                                                                                                         <  3:294241/149‑1 (MQ=255)
                    ccTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACg                                                                                                                                         <  4:286756/136‑1 (MQ=255)
                    ccTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACg                                                                                                                                         >  3:286756/1‑136 (MQ=255)
                      tGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGGTATGACGGCGCGCCTGATCGa                                                                                                                           >  5:143674/1‑148 (MQ=255)
                      tGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGGTATGACGGCGCGCCTGATCGa                                                                                                                           >  5:137772/1‑148 (MQ=255)
                       gCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGaa                                                                                                                          >  4:330610/1‑148 (MQ=255)
                        cGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGaaaa                                                                                                                        <  7:125061/149‑1 (MQ=255)
                        cGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGaa                                                                                                                          >  6:166089/1‑147 (MQ=255)
                          tATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAgg                                                                                                                      <  2:32218/149‑1 (MQ=255)
                          tATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAgg                                                                                                                      <  5:166089/149‑1 (MQ=255)
                               gTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATc                                                                                                                 >  5:348930/1‑149 (MQ=255)
                                tGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATcc                                                                                                                >  2:100480/1‑149 (MQ=255)
                                    ctctGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCctc                                                                                                            >  3:71375/1‑149 (MQ=255)
                                     tctGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCctct                                                                                                           >  1:223290/1‑149 (MQ=255)
                                       tGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCGTc                                                                                                         >  8:128101/1‑149 (MQ=255)
                                         cGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCg                                                                                                        <  6:169578/148‑1 (MQ=255)
                                          ggTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGt                                                                                                      <  7:321784/149‑1 (MQ=255)
                                                  cGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGaaaa                                                                                                                        <  1:57420/123‑1 (MQ=255)
                                                  cGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGaaaa                                                                                                                        >  2:57420/1‑123 (MQ=255)
                                                      cTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATcccc                                                                                                              <  4:17649/129‑1 (MQ=255)
                                                      cTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATcccc                                                                                                              >  3:17649/1‑129 (MQ=255)
                                                           gcgcCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGc                                                                                      >  4:111303/1‑148 (MQ=255)
                                                                gTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCg                                                                                <  2:45693/149‑1 (MQ=255)
                                                                    cGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCggtgg                                                                            <  1:349332/149‑1 (MQ=255)
                                                                     gAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCggtggt                                                                           >  2:160296/1‑149 (MQ=255)
                                                                           cTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTgctgct                                                                     <  3:330610/149‑1 (MQ=255)
                                                                             ggCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTgg                                                                   >  1:190179/1‑149 (MQ=255)
                                                                              gCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTggg                                                                  <  7:72040/149‑1 (MQ=255)
                                                                               cAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTgggg                                                                 <  8:148379/149‑1 (MQ=255)
                                                                                aGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATccc                                                                                                               >  8:18906/1‑102 (MQ=255)
                                                                                aGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATccc                                                                                                               <  7:18906/102‑1 (MQ=255)
                                                                                aGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGc                                                                <  1:160296/149‑1 (MQ=255)
                                                                                    cAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGtt                                                            >  1:126193/1‑149 (MQ=255)
                                                                                     aGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCTAGCCGGTGGTGCTGCTGGGGCCGTTc                                                           <  1:202276/149‑1 (MQ=255)
                                                                                     aGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTc                                                           >  4:88414/1‑149 (MQ=255)
                                                                                      gCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTcc                                                          <  1:100480/149‑1 (MQ=255)
                                                                                                     gTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGaagaag                                           <  6:348930/149‑1 (MQ=255)
                                                                                                      tGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGGTATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAgc                                          <  6:137772/149‑1 (MQ=255)
                                                                                                      tGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGGTATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAgc                                          <  6:143674/149‑1 (MQ=255)
                                                                                                         gATGTGGAAGTGGTACAGGGCGCCAACACGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCATCTTCGGTGGAGTAGCGTGACGGCAAGCCGTTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAgcg                                         >  1:365786/1‑147 (MQ=255)
                                                                                                              ggAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgc                                  >  8:253877/1‑149 (MQ=255)
                                                                                                              ggAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgc                                  >  3:364918/1‑149 (MQ=255)
                                                                                                                aaGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgcg                                 >  4:367732/1‑148 (MQ=255)
                                                                                                                     gTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGa                           >  5:71697/1‑149 (MQ=255)
                                                                                                                     gTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAAGGa                           <  3:111303/149‑1 (MQ=255)
                                                                                                                      tACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGAc                          >  4:4581/1‑149 (MQ=255)
                                                                                                                             cgcCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGa                   >  8:288217/1‑149 (MQ=255)
                                                                                                                               ccAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATg                 <  3:367732/149‑1 (MQ=255)
                                                                                                                                cAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATgg                >  2:192128/1‑149 (MQ=255)
                                                                                                                                      gggCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGc          >  4:85466/1‑149 (MQ=255)
                                                                                                                                        gcgcACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTg        >  5:102595/1‑149 (MQ=255)
                                                                                                                                              tCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTCGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTAcc  >  3:210173/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTCGACGATACCGGCAACAACCTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACATGGGCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACC  >  NC_002947/2797191‑2797481

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