Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,563,519 G→A 7.7% A421T (GCC→ACC)  PP_1372 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,563,5190GA 7.7% 149.8 / 2.3 52A421T (GCC→ACC) PP_1372hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (2/2);  total (25/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 = 2.43e-01

GAGAAGGGCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCAAGTCGGCAGTGCGCCAGGC  >  NC_002947/1563372‑1563664
                                                                                                                                                   |                                                                                                                                                 
gagaAGGGCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGc                                                                                                                                                  >  7:35681/1‑149 (MQ=255)
gagaAGGGCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGc                                                                                                                                                  <  4:178858/149‑1 (MQ=255)
gagaAGGGCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGc                                                                                                                                                  <  4:174773/149‑1 (MQ=255)
       gCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTTCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCt                                                                                                                                                <  1:245646/144‑1 (MQ=255)
       gCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCt                                                                                                                                                >  2:245646/1‑144 (MQ=255)
              ggCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGa                                                                                                                                    <  4:258156/149‑1 (MQ=255)
                 cATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCTGTGGACCGTGAGTCAGCCTATGAAATGCTGAGcc                                                                                                                                 <  5:160930/149‑1 (MQ=255)
                          gagGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAgg                                                                                                                        >  2:204793/1‑149 (MQ=255)
                              gTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTTATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGa                                                                                                                    <  1:61631/149‑1 (MQ=255)
                               tGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAg                                                                                                                   >  7:6731/1‑149 (MQ=255)
                                 cTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGc                                                                                                                      <  5:259301/144‑1 (MQ=255)
                                 cTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGc                                                                                                                      >  6:259301/1‑144 (MQ=255)
                                    gTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCgg                                                                                                               <  4:108575/148‑1 (MQ=255)
                                        cACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAACCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGaa                                                                                                          >  4:57504/1‑149 (MQ=255)
                                        cACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAACCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGaa                                                                                                          >  4:53154/1‑149 (MQ=255)
                                         aCCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCgg                                                                                                               >  8:26674/1‑143 (MQ=255)
                                         aCCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCgg                                                                                                               <  7:26674/143‑1 (MQ=255)
                                         aCCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAAc                                                                                                         >  6:94198/1‑149 (MQ=255)
                                                 aGTCGCGATGCGGGCCACGGCGTGCGGCGGAGCGCAGTGCGCGGATGGCCGCGTCGCCACTGGCGGGGCGCTATGACACACCGGGGCACCGGGAGGCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGcgccg                                                                                                 <  5:187182/149‑1 (MQ=255)
                                                            gggCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGc                                                                                      <  6:268180/149‑1 (MQ=255)
                                                            gggCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTGCGCCACTGGCGGGGCGCTATGACAACCCGGTGGACCGGGAGTCAGCCTATGAACTGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACAGAAGGc                                                                                      <  5:258182/149‑1 (MQ=255)
                                                             ggCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCt                                                                                                                                      >  5:248264/1‑100 (MQ=255)
                                                             ggCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCt                                                                                                                                      <  6:248264/100‑1 (MQ=255)
                                                             ggCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAAc                                                                                                         <  2:44755/129‑1 (MQ=255)
                                                             ggCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATAGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAAc                                                                                                         >  1:44755/1‑129 (MQ=255)
                                                                cACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGAt                                                                                  >  6:229135/1‑149 (MQ=255)
                                                                  cTGTGTGCGGATGAGTGTAGTGCGCTGATTGCCGGTTCGCCTCTGGGGGGGCGCTATGACATTCCGGTTGACCGTGAGTCAGCCTATGAAATGCTGAGACAGCGTAAGGGCGAGCCGGTg                                                                                                             <  8:263931/120‑1 (MQ=255)
                                                                  cTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCAACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGATCCGGTg                                                                                                             >  7:263931/1‑120 (MQ=255)
                                                                          ggCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGTACCGTGAGTCAACCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGtt                                                                        <  3:53154/149‑1 (MQ=255)
                                                                          ggCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAACCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGtt                                                                        <  3:57504/149‑1 (MQ=255)
                                                                                  gcAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGa                                                                 >  2:108343/1‑148 (MQ=255)
                                                                                         gcTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCgg                                                         >  1:50032/1‑149 (MQ=255)
                                                                                             aTTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGa                                                     <  1:204793/149‑1 (MQ=255)
                                                                                               tGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGt                                                   >  8:65038/1‑149 (MQ=255)
                                                                                                   gCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCtt                                               >  4:38031/1‑149 (MQ=255)
                                                                                                      tCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTg                                              >  4:163394/1‑147 (MQ=255)
                                                                                                        gCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCaga                                          <  3:38031/149‑1 (MQ=255)
                                                                                                         ccACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCagag                                         <  1:108343/149‑1 (MQ=255)
                                                                                                          cACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGt                                        >  2:129642/1‑149 (MQ=255)
                                                                                                           aCTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTg                                       <  5:229135/149‑1 (MQ=255)
                                                                                                                  gggCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCggg                                <  7:65038/149‑1 (MQ=255)
                                                                                                                  gggCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCggg                                <  8:35681/149‑1 (MQ=255)
                                                                                                                   ggCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCgggg                               <  3:163394/149‑1 (MQ=255)
                                                                                                                        tATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGAGAAGGCGGGt                                                       <  8:78505/120‑1 (MQ=255)
                                                                                                                        tATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGt                                                       >  7:78505/1‑120 (MQ=255)
                                                                                                                          tGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGa                        >  7:244420/1‑149 (MQ=255)
                                                                                                                              taaCCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCaa                    <  5:215908/148‑1 (MQ=255)
                                                                                                                               aaCCGGTCGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGAAACGGCGAGCCGGAGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGACTTCTTGCAGAGTGCAGCGGGGCACGCGATc                      >  4:377/1‑146 (MQ=255)
                                                                                                                                aCCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCaa                    >  6:215908/1‑147 (MQ=255)
                                                                                                                                aCCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCAAGt                  >  8:283424/1‑149 (MQ=255)
                                                                                                                                     tGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCAAGTCGGCa             <  8:6731/149‑1 (MQ=255)
                                                                                                                                                tCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCAAGTCGGCAGTGCGCCAGGc  >  3:3511/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAGAAGGGCACGCCGGCCATGGTGCAGAGGGTGCTGGTTGCACCGCCGCAGTCGCGTATCGGGCCACTGAGTGCGGCTGAGCGCAGTGCGCTGATTGCCGCTTCGCCACTGGCGGGGCGCTATGACAAACCGGTGGACCGTGAGTCAGCCTATGAAATGCTGAGCCAGCGTAAGGGCGAGCCGGTGGAACCGGCGCCGCAACCGAAGGCGGATGAAGAGAGTTTTGCCGACAAGGCGGGTGAGTTCTTGCAGAGTGCAGCGGGGCAGGCGATCAAGTCGGCAGTGCGCCAGGC  >  NC_002947/1563372‑1563664

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