Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,138,575 A→C 7.0% D25A (GAC→GCC)  PP_3643 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,138,5750AC 7.0% 161.4 / 2.6 54D25A (GAC→GCC) PP_3643oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base C (2/2);  total (26/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 = 9.99e-01

CGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGT  >  NC_002947/4138435‑4138721
                                                                                                                                             |                                                                                                                                                  
cGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGgacgacgac                                                                                                                                             <  7:142839/148‑1 (MQ=255)
     taCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCg                                                                                                                                       <  8:163389/149‑1 (MQ=255)
     taCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCg                                                                                                                                       <  6:168876/149‑1 (MQ=255)
        cATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGgtc                                                                                                                                    >  2:226365/1‑149 (MQ=255)
               cacaCCCACGATACAAACCGCCGACCACCGCAGCCCTTCAGACGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGcc                                                                                                                             <  6:244383/149‑1 (MQ=255)
                  aTCAATGATAAAATACTCGTACCACCGCAGCCCTTCAGATGTACTCGCGATGATCGAGGCAACTGTATAGAAGATCGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCT‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGa                                                                                                                          <  5:268244/149‑1 (MQ=255)
                  aTCAACGATAAAAACAGCGGACAACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACAC‑ACCAGGACGACGACGGCTCGGTCGTCGGCCTGa                                                                                                                          >  2:272372/1‑149 (MQ=255)
                   tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACGCCCGATGATCGACGACACTGTCTCCACGCACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAAGCCGACGACGGCTCGGTCGTAGGCCTGAc                                                                                                                         >  5:231984/1‑149 (MQ=255)
                     aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCg                                                                                                                       >  1:190943/1‑149 (MQ=255)
                     aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCg                                                                                                                       >  2:218335/1‑149 (MQ=255)
                      aCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGCCGACGACCGAGCCGTCGTCGTCCTgggtgg                                                                                                                      <  3:296965/149‑6 (MQ=38)
                               aCCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCtt                                                                                                             <  1:226365/149‑1 (MQ=255)
                                 cgcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCt                                                                                                                            >  5:13622/1‑132 (MQ=255)
                                 cgcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCt                                                                                                                            <  6:13622/132‑1 (MQ=255)
                                 cgcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCACGACGACGACGGCACGGTCGTCGGACTGACCGCCAGCCCCTtca                                                                                                           >  2:175330/1‑148 (MQ=255)
                                  gcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTAGAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtctc                                                                                                          <  1:264285/149‑1 (MQ=255)
                                      accaccGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGACTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCtcctcc                                                                                                      <  3:172871/149‑1 (MQ=255)
                                           cGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTc                                                                                                 >  3:296946/1‑149 (MQ=255)
                                                cccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCATCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTTACCGCCAGCGCCTTCTACTCCCTGTCGATg                                                                                             >  7:250227/1‑148 (MQ=255)
                                                 ccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGCCGACGACCGAGCCGTCGTCGTCCt                                                                                                                            <  2:17372/116‑1 (MQ=38)
                                                 ccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGCCGACGACCGAGCCGTCGTCGTCCt                                                                                                                            >  1:17372/1‑116 (MQ=38)
                                                      agagGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgc                                                                                      <  4:296946/149‑1 (MQ=255)
                                                           tACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTg                                                                                 <  2:190943/149‑1 (MQ=255)
                                                             cTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGt                                                                               >  3:188848/1‑149 (MQ=255)
                                                                ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCt                                                                            >  6:259649/1‑149 (MQ=255)
                                                                           cGCAACTGTCTACAAGAACGGTCGGGGCTCCTTCCCGTCCGGCGGTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGccc                                                                   <  4:188848/147‑1 (MQ=255)
                                                                                cTGTATAAAAGAACTTTGTGGGCTCCTTACCGTCCGGGGTTACCGTCATCTCCA‑CCCAGGACGACGAGGGCTCGGTCGTCGGACTGACCGCCAGCGCCTTCTCCTCCCTGTCGATTGACCCGCCGCTGGTGCTGTTCTGCCCCAAc                                                               <  8:250227/146‑1 (MQ=37)
                                                                                   tCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTacac                                                          <  3:23318/148‑1 (MQ=255)
                                                                                   tCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTaca                                                           <  1:218335/147‑1 (MQ=255)
                                                                                     tACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGACTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       >  6:133598/1‑149 (MQ=255)
                                                                                      aCAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGcc                                                                    >  1:164780/1‑135 (MQ=255)
                                                                                      aCAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGcc                                                                    <  2:164780/135‑1 (MQ=255)
                                                                                          gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       <  4:97022/144‑1 (MQ=255)
                                                                                          gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       >  3:97022/1‑144 (MQ=255)
                                                                                          gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       >  5:323408/1‑144 (MQ=255)
                                                                                          gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       <  6:323408/144‑1 (MQ=255)
                                                                                               ttCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGc                                                                     <  4:79567/125‑1 (MQ=255)
                                                                                               ttCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGc                                                                     >  3:79567/1‑125 (MQ=255)
                                                                                               ttCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       >  1:122156/1‑139 (MQ=255)
                                                                                               ttCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTc                                                       <  2:122156/139‑1 (MQ=255)
                                                                                                 cTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTAcc                                            <  3:159659/148‑1 (MQ=255)
                                                                                                     gCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                       <  3:218487/149‑1 (MQ=255)
                                                                                                          ttCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGa                                   >  7:105711/1‑148 (MQ=255)
                                                                                                            cccGTCCGGCGTTACCGTCATCACCAGCCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTAc                                      >  1:167187/1‑144 (MQ=255)
                                                                                                            cccGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTAc                                      <  2:167187/143‑1 (MQ=255)
                                                                                                                           cGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCaa                                                                >  7:126921/1‑102 (MQ=255)
                                                                                                                           cGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCaa                                                                <  8:126921/102‑1 (MQ=255)
                                                                                                                            gTCATCACCA‑CGCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGc                <  1:175330/149‑1 (MQ=255)
                                                                                                                                tcaCCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCAtt            >  6:73518/1‑149 (MQ=255)
                                                                                                                                 cacca‑cCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCAtt            >  5:94577/1‑148 (MQ=255)
                                                                                                                                   cca‑ccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATc                                 >  6:109824/1‑125 (MQ=255)
                                                                                                                                   cca‑ccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATc                                 <  5:109824/125‑1 (MQ=255)
                                                                                                                                    ca‑ccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCAcc        >  6:32942/1‑149 (MQ=255)
                                                                                                                                     a‑ccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACct       <  5:32942/149‑1 (MQ=255)
                                                                                                                                           ggACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGt  >  1:214277/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCA‑CCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGT  >  NC_002947/4138435‑4138721

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