Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,639,903 C→G 8.3% A291A (GCC→GCG)  PP_3206 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,639,9030CG 8.3% 141.1 / 2.9 48A291A (GCC→GCG) PP_3206hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base G (2/2);  total (22/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 = 8.73e-01

TTGCATTTGACCATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAACACACGTCGCTTGCTGGACCTGTCTGT  >  NC_002947/3639763‑3640051
                                                                                                                                            |                                                                                                                                                    
ttGCATTTGACCATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCGGGCGAAcga                                                                                                                                             <  2:156644/150‑3 (MQ=255)
       tGACCATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCt                                                                                                                                      <  1:926612/150‑1 (MQ=255)
         aCCATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCCTCTGc                                                                                                                                    >  1:243172/1‑150 (MQ=255)
           cATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCGGGCGAAcga                                                                                                                                             >  1:1331686/1‑137 (MQ=255)
           cATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACATCGTTCGCGGGCGAAcga                                                                                                                                             <  2:1331686/139‑3 (MQ=255)
                  cAGGTGCTGACCGGGCTGGCGGCTGTTGCTCGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGTCTGTTTGCCCGCTACCCGCCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCGTCTGCCACGACGGt                                                                                                                           <  2:1198756/150‑1 (MQ=255)
                          gACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTc                                                                                                                                       >  1:369784/1‑130 (MQ=255)
                          gACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTc                                                                                                                                       <  2:369784/130‑1 (MQ=255)
                              gggCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTggccggc                                                                                                               <  2:360491/150‑1 (MQ=255)
                               ggCTGGCGGCTGTTGTTAGCGCGAACCCCCGAAATCGCACCGCCTTCGCCCCCTATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTggccggcc                                                                                                              <  2:513714/150‑1 (MQ=255)
                                     cggcTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCtgctgc                                                                                                        <  2:1010998/150‑1 (MQ=255)
                                      ggcTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCg                                                                                                       >  1:683474/1‑150 (MQ=255)
                                          gttggcggCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGTTCTCCAGCCCGAAAGCCTCTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGGTTCTGCCACGACGTCGCGGTGGCAGGAATGCTGCCTAAt                                                                                                   <  2:54896/144‑1 (MQ=255)
                                             gCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGcc                                                                                                >  1:86549/1‑150 (MQ=255)
                                                  cctcgCCAGCCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGa                                                                                           <  1:1151737/146‑1 (MQ=255)
                                                          aCCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTgccgc                                                                                   >  1:1216687/1‑150 (MQ=255)
                                                                  cgcATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCcac                                                                           >  2:315083/1‑150 (MQ=255)
                                                                          cTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTcg                                                                   >  1:276435/1‑150 (MQ=255)
                                                                                ccccGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCa                                                                                                                        <  1:397912/91‑1 (MQ=255)
                                                                                ccccGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCa                                                                                                                        >  2:397912/1‑91 (MQ=255)
                                                                                ccccGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCcgccgc                                                                                >  2:990372/1‑131 (MQ=255)
                                                                                ccccGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCcgccgc                                                                                <  1:990372/131‑1 (MQ=255)
                                                                                  ccGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATg                                                                                                  >  1:1234617/1‑111 (MQ=255)
                                                                                  ccGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATg                                                                                                  <  2:1234617/111‑1 (MQ=255)
                                                                                      tCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACtt                                                                                        >  1:227282/1‑117 (MQ=255)
                                                                                      tCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACtt                                                                                        <  2:227282/117‑1 (MQ=255)
                                                                                       ccccAGCGCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCGCTCGACTGGCCGCCGCGCCACACCCGTCTCGCGGCCATCCGc                                                      >  1:1274133/1‑150 (MQ=255)
                                                                                          cAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCaaa                                                   <  1:315083/150‑1 (MQ=255)
                                                                                              tCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTc                                                                                             >  2:776627/1‑104 (MQ=255)
                                                                                              tCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTc                                                                                             <  1:776627/104‑1 (MQ=255)
                                                                                                     ccTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCacac                                                                         >  2:882254/1‑117 (MQ=255)
                                                                                                     ccTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCacac                                                                         <  1:882254/117‑1 (MQ=255)
                                                                                                        gTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATg                                     >  1:1151849/1‑150 (MQ=255)
                                                                                                           tGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGTCCTTCCGCAAAGGAATCCTCGTATGAGc                                  <  1:616612/150‑1 (MQ=255)
                                                                                                              ccGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTg                               <  2:833794/150‑1 (MQ=255)
                                                                                                                gCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACTCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGaa                             <  2:1031685/150‑1 (MQ=255)
                                                                                                                gCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGaa                             <  2:1216687/150‑1 (MQ=255)
                                                                                                                gCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGaa                             <  2:86549/150‑1 (MQ=255)
                                                                                                                     ccACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAacacac                        >  1:166568/1‑150 (MQ=255)
                                                                                                                     ccACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAacacac                        >  1:1115964/1‑150 (MQ=255)
                                                                                                                       aCCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGcg                                                            <  2:1078031/112‑2 (MQ=255)
                                                                                                                       aCCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGcg                                                            >  1:1078031/1‑111 (MQ=255)
                                                                                                                                tACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAacacac                        >  2:378008/1‑139 (MQ=255)
                                                                                                                                tACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAacacac                        <  1:378008/139‑1 (MQ=255)
                                                                                                                                  cccagTTCGCGGGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAacaca                         >  2:933284/5‑136 (MQ=255)
                                                                                                                                  ccTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAg                                   <  2:257065/126‑1 (MQ=255)
                                                                                                                                  ccTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAg                                   >  1:257065/1‑126 (MQ=255)
                                                                                                                                         ggCCCGCGAACTGGGCTTCTGCCACGCAGGTCACGTAGCCGGCATTCTTCTAAATGGCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAACACATGTCGCTTGCTGGACctgtct    <  2:233956/149‑1 (MQ=255)
                                                                                                                                           cccGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAACACACGTCGCTTGCTGGACctgtctgt  <  2:276435/150‑1 (MQ=255)
                                                                                                                                           cccGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAACACACGTCGCTTGCTGGACctgtctgt  <  2:807002/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
TTGCATTTGACCATCGCCCAGGTGCTGACCGGGCTGGCGGCTGTTGCTGGCGAGGCCAACCGTGAGCGCATCGCCTTCGCCCCCGATCCCCAGCTCGAAAGCCTGTTTGGCCGCTACCCACCGTTGCGTACCTCGTTCGCCCGCGAACTGGGCTTCTGCCACGACGGTTCAGTGGCCGGCCTGCTGCGTAATGCCCTCGACTTGCCGCCGCGCCACACCCGTCGCGCGGCCATCCGCAAAGGAATCCTCGTATGAGCTTGAACACACGTCGCTTGCTGGACCTGTCTGT  >  NC_002947/3639763‑3640051

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑