Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 333,177 G→T 7.8% A106S (GCC→TCC)  PP_5435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947333,1770GT 7.8% 152.4 / 3.7 51A106S (GCC→TCC) PP_5435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base T (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAT  >  NC_002947/333032‑333325
                                                                                                                                                 |                                                                                                                                                    
cATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGcc                                                                                                                                                    >  6:190082/1‑148 (MQ=255)
     cccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                              >  6:260919/1‑149 (MQ=255)
      ccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                              >  3:40668/1‑148 (MQ=255)
      ccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                              >  7:18560/1‑148 (MQ=255)
          tCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACa                                                                                                                                          >  6:21725/1‑148 (MQ=255)
            aaaaCTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAAt                                                                                                                                       >  8:4815/1‑149 (MQ=255)
             aaaCTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATg                                                                                                                                      <  1:4754/149‑1 (MQ=255)
              aaCTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATgg                                                                                                                                     <  5:154341/149‑1 (MQ=255)
                 tCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTg                                                                                                                                   >  5:18938/1‑148 (MQ=255)
                 tCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCCCCTGATCGCGTGGCTCAATGCACGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCCCGCCGAGTACGACAAATGGTGc                                                                                                                                  <  6:284674/149‑1 (MQ=255)
                   aGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGc                                                                                                                                  <  7:138220/147‑1 (MQ=255)
                   aGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGc                                                                                                                                  <  4:203498/147‑1 (MQ=255)
                   aGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCa                                                                                                                                 <  8:18560/148‑1 (MQ=255)
                      cccGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCGGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGaa                                                                                                                              <  4:335585/148‑1 (MQ=255)
                        cGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCg                                                                                                                           <  6:18938/149‑1 (MQ=255)
                         gCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGa                                                                                                                          >  7:235435/1‑149 (MQ=255)
                                  aaaGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTa                                                                                                                 <  7:301960/149‑1 (MQ=255)
                                  aaaGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGACTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTa                                                                                                                 <  7:311488/149‑1 (MQ=255)
                                                      cGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGa                                                                                               <  8:235435/147‑1 (MQ=255)
                                                        tgtTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATAtt                                                                                           <  6:233458/149‑1 (MQ=255)
                                                         gtTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGCACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATAttt                                                                                          <  3:249301/149‑1 (MQ=255)
                                                              atatCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAAtt                                                                                     <  7:4815/149‑1 (MQ=255)
                                                               tatCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTc                                                                                    <  5:21725/149‑1 (MQ=255)
                                                                   aaGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                                  >  1:380422/1‑147 (MQ=255)
                                                                   aaGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACa                                                                                <  2:380422/149‑1 (MQ=255)
                                                                    aGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCcgc                                                                                                   <  6:354516/129‑1 (MQ=255)
                                                                    aGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAg                                                                               <  5:190082/149‑1 (MQ=255)
                                                                     gAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAgg                                                                              >  7:233350/1‑149 (MQ=255)
                                                                        aaaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCa                                                                                        <  4:71211/136‑1 (MQ=255)
                                                                        aaaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCa                                                                                        >  3:71211/1‑136 (MQ=255)
                                                                         aaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGa                                                                           >  8:180448/1‑148 (MQ=255)
                                                                           cGAATGCAGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGacaa                                                                        <  4:40668/149‑1 (MQ=255)
                                                                                      tGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                                  <  2:204160/128‑1 (MQ=255)
                                                                                      tGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                                  >  1:204160/1‑128 (MQ=255)
                                                                                                     aTGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                <  2:264261/147‑1 (MQ=255)
                                                                                                     aTGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                >  1:264261/1‑147 (MQ=255)
                                                                                                               gggTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGa                                    <  5:260919/149‑1 (MQ=255)
                                                                                                                 gTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCGACTGGACTAGATg                                  >  8:105210/1‑149 (MQ=255)
                                                                                                                   aCTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                <  1:175019/133‑1 (MQ=255)
                                                                                                                   aCTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                >  2:175019/1‑133 (MQ=255)
                                                                                                                       tctACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAAAGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCg                                                               >  6:335947/1‑114 (MQ=255)
                                                                                                                            aTCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTg                                                                                                         >  8:107035/1‑67 (MQ=255)
                                                                                                                               aaGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTaa                     <  7:180448/148‑1 (MQ=255)
                                                                                                                                   ccATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTaa                     >  2:124753/1‑144 (MQ=255)
                                                                                                                                   ccATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTaa                     <  1:124753/144‑1 (MQ=255)
                                                                                                                                   accatttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATg                                  <  4:261201/127‑1 (MQ=255)
                                                                                                                                   accatttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATg                                  >  3:261201/5‑131 (MQ=255)
                                                                                                                                       tttGTCGTACTCGGCGTGCGACAAATGGTTCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAAt                                                                                      <  1:268393/75‑1 (MQ=25)
                                                                                                                                       tttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAAt                                                                                      >  2:268393/1‑75 (MQ=37)
                                                                                                                                            cGCACGCCGAGTACGACAAATTGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGc       >  6:138131/1‑149 (MQ=255)
                                                                                                                                             gCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATATCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCAATGCt      >  4:47143/1‑149 (MQ=255)
                                                                                                                                                 gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTggg                                                                                                       >  4:273385/1‑48 (MQ=255)
                                                                                                                                                 gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTggg                                                                                                       <  3:273385/48‑1 (MQ=255)
                                                                                                                                                 gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAt  >  6:380132/1‑149 (MQ=255)
                                                                                                                                                 gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAt  >  8:53262/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAT  >  NC_002947/333032‑333325

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