Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,217 G→C 8.9% K62N (AAG→AAC)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,2170GC 8.9% 164.6 / 5.7 54K62N (AAG→AAC) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base C (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCA  >  NC_002947/571070‑571357
                                                                                                                                                   |                                                                                                                                            
cTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGt                                                                                                                                            >  1:864874/1‑150 (MQ=255)
   tgGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGc                                                                                                                                         >  2:552555/1‑150 (MQ=255)
     gCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCg                                                                                                                                       >  1:797562/1‑150 (MQ=255)
       cAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGaa                                                                                                                                     >  2:629302/1‑150 (MQ=255)
         ggATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaaca                                                                                                                                   <  2:886466/150‑3 (MQ=255)
             cccGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCt                                                                                                                               >  2:880732/1‑150 (MQ=255)
               cGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                             <  1:110288/150‑1 (MQ=255)
               cGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                             >  2:1178039/1‑150 (MQ=255)
                  aGTCCGCTCCCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGAGGGCCGTCTACGGTCCAGCAGTCGGGCCAGAAGGTCGAGAGTACCGAGTGGCACCGGTCGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGa                                                                                                                          <  1:1178039/150‑1 (MQ=255)
                   gTCCGCTACCTACCCATCTGTATCGCCGTGACCAACCTGATCCTGGCTACCAGCGAGCAGTGAACTCACAACCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCTGCAAGTTTGCCGAAATCGCTGGCGaa                                                                                                                         <  2:468737/150‑1 (MQ=255)
                      cGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATAc                                                                                                                      >  1:461429/1‑150 (MQ=255)
                               cccAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGaa                                                                                                                         <  2:1040598/138‑1 (MQ=255)
                               cccAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGaa                                                                                                                         >  1:1040598/1‑138 (MQ=255)
                                      gTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGc                                                                                                      >  2:272261/1‑150 (MQ=255)
                                         aCGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGt                                                                                                   >  2:706122/1‑150 (MQ=255)
                                         aCGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGt                                                                                                   <  2:54023/150‑1 (MQ=255)
                                            ccGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAAAAGGTCGAGTGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCt                                                                                                <  2:342758/150‑1 (MQ=255)
                                              gTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTAc                                                                                              <  1:1047692/150‑1 (MQ=255)
                                                 aCCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTg                                                                                                  >  1:499987/1‑143 (MQ=255)
                                                 aCCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTg                                                                                                  <  2:499987/143‑1 (MQ=255)
                                                 aCCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATc                                                                                           >  2:1021385/1‑150 (MQ=255)
                                                     aCCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGcc                                                                                                                                        >  1:559388/1‑101 (MQ=255)
                                                     aCCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGcc                                                                                                                                        <  2:559388/101‑1 (MQ=255)
                                                       cTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAACGCCTGAGCAAGTGTTCGCAGTGCTACATCGAAGGc                                                                                     <  1:285929/150‑1 (MQ=255)
                                                         gAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCaa                                                                                   <  2:398206/150‑1 (MQ=255)
                                                          aGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTTGCAGTGCTACATCGAAGGCAAg                                                                                  <  2:797562/150‑1 (MQ=255)
                                                               ggCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGAAGTGCTACATCCAAGGCAAGCGGCa                                                                             >  2:672266/1‑150 (MQ=255)
                                                                    ccaGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGAccc                                                                        >  2:382985/1‑150 (MQ=255)
                                                                      aGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCcgc                                                                      <  1:552555/150‑1 (MQ=255)
                                                                            cAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTgg                                                                <  2:549678/150‑1 (MQ=255)
                                                                                   cTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAg                                                         >  2:382611/1‑150 (MQ=255)
                                                                                    tGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGa                                                        <  1:706122/150‑1 (MQ=255)
                                                                                     gACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGAc                                                       <  1:382985/150‑1 (MQ=255)
                                                                                      aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGc                                                                              >  1:307554/1‑126 (MQ=255)
                                                                                      aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGc                                                                              <  2:307554/126‑1 (MQ=255)
                                                                                      aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACg                                                      >  1:1026958/1‑150 (MQ=255)
                                                                                      aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGGGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACg                                                      >  2:239811/1‑150 (MQ=255)
                                                                                       cAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCTCGATTGGGAAAAAGACgg                                                     <  2:864874/150‑1 (MQ=255)
                                                                                       cAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACgg                                                     <  1:382611/150‑1 (MQ=255)
                                                                                        aaGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGc                                                    >  2:657707/1‑150 (MQ=255)
                                                                                                      gAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTacac                                      >  1:423973/1‑150 (MQ=255)
                                                                                                         ggTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTAcaccac                                   <  1:1021385/150‑1 (MQ=255)
                                                                                                          gTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGCGAAAAAGACGGCATCAAGCGCTACaccacc                                  <  1:629302/150‑1 (MQ=255)
                                                                                                              aGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGTCATCAAGCGCTACACCACCgaga                              <  2:461429/150‑1 (MQ=255)
                                                                                                              aGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCgaga                              <  1:672266/150‑1 (MQ=255)
                                                                                                                  tACCGAGTGGCACCGTGTGTCGCTGTTCGGACAAGTTGCCAAAATCCCTGTCGAATACCTGCGCAAGGGGTCGCAGCGCCCCACACCTCCCGACCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGatcat                          >  1:574986/1‑150 (MQ=255)
                                                                                                                                    ggcGCTGCTCGGCAACGTTGCCGAAATCGCGGGCGAATACCTGAGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAAACCCGCGAGTGAGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCAc        <  2:423973/148‑1 (MQ=255)
                                                                                                                                    cagcgatttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATAATCGTCGGTAACAATGGAAc        >  2:670834/8‑150 (MQ=255)
                                                                                                                                     tCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCAcc       >  2:115889/1‑150 (MQ=255)
                                                                                                                                      gcgatttCGGCAACCTTGCCGAACTCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCa      >  1:1173933/6‑150 (MQ=255)
                                                                                                                                       gCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCAt     <  1:638219/150‑1 (MQ=255)
                                                                                                                                        cTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATg    <  1:115889/150‑1 (MQ=255)
                                                                                                                                         tGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGc   <  1:880732/150‑1 (MQ=255)
                                                                                                                                          tttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGAGATCAACGGCACCATGCa  >  1:244702/2‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCA  >  NC_002947/571070‑571357

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