Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,218 G→C 12.0% V63L (GTT→CTT)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,2180GC12.0% 160.6 / 11.8 56V63L (GTT→CTT) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base C (3/4);  total (28/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.22e-01

CCTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTCG  >  NC_002947/571069‑571365
                                                                                                                                                     |                                                                                                                                                    
ccTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAgg                                                                                                                                                      <  1:604418/150‑1 (MQ=255)
        cAGGATTCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGaa                                                                                                                                             >  2:434088/1‑150 (MQ=255)
            aTCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCg                                                                                                                                         >  2:1009510/1‑150 (MQ=255)
               ccGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGtt                                                                                                                                                   <  2:327966/137‑1 (MQ=255)
               ccGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGtt                                                                                                                                                   >  1:327966/1‑137 (MQ=255)
               ccGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTg                                                                                                                                      >  2:300872/1‑150 (MQ=255)
                cGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                                     >  1:420118/1‑150 (MQ=255)
                cGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                                     >  1:305555/1‑150 (MQ=255)
                cGAAGTCCGCTACCTGCCCAACGGTAACGCCGTCACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGCGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTgg                                                                                                                                     >  1:819885/1‑150 (MQ=255)
                     tCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAAt                                                                                                                                <  2:916728/150‑1 (MQ=255)
                      ccGCTACCTGCCCAACGGTACCGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCCTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATa                                                                                                                               <  1:893770/150‑1 (MQ=255)
                       cGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATAc                                                                                                                              <  1:924882/150‑1 (MQ=255)
                            ccTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTgcg                                                                                                                         >  2:1005425/1‑150 (MQ=255)
                                 ccAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAgg                                                                                                                    <  2:183/150‑1 (MQ=255)
                                   aaCGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGt                                                                                                                  >  1:126881/1‑150 (MQ=255)
                                      ggTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCAAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGaa                                                                                                                                 <  2:499437/132‑1 (MQ=255)
                                      ggTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCAAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGaa                                                                                                                                 >  1:499437/1‑132 (MQ=255)
                                           cGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCGATCATTCGAGAAGGGCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGcncgc                                                                                                          <  2:705045/150‑5 (MQ=255)
                                             ccGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaacagc                                                                                                                                         <  1:304355/117‑5 (MQ=255)
                                             ccGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAACCTTGCCGaacagc                                                                                                                                         >  2:304355/1‑113 (MQ=255)
                                                 gACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGAAACGTGAGTCGCTGTTCGGAAGCCTGTGCCGTTTCTCGTGCGAATACCTGCGCAAGGGTTCGCAGTGCTACAt                                                                                                    <  1:300872/150‑1 (MQ=255)
                                                   ccAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCg                                                                                                  <  1:439646/150‑1 (MQ=255)
                                                    cAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATAACTGCGAAACGGTTCGCAAGGCTACATTGa                                                                                                 >  1:19893/1‑150 (MQ=255)
                                                       ccTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGGGCTACATCGAAgg                                                                                              >  2:112801/1‑150 (MQ=255)
                                                               tGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCCCAGTGTGTAGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGc                                                                                      >  2:678471/1‑150 (MQ=255)
                                                                   gaccaGCTAGCTGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGAc                                                                                  >  2:856649/2‑150 (MQ=255)
                                                                         cGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGa                                                                            <  1:1005425/150‑1 (MQ=255)
                                                                               gTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGa                                                                      >  2:195924/1‑150 (MQ=255)
                                                                                  gACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTgg                                                                        <  1:1109345/145‑1 (MQ=255)
                                                                                  gACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTgg                                                                        >  2:1109345/1‑145 (MQ=255)
                                                                                       aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTgg                                                                        <  1:1081810/140‑1 (MQ=255)
                                                                                       aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTgg                                                                        >  2:1081810/1‑140 (MQ=255)
                                                                                       aCAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACg                                                              >  1:1057602/1‑150 (MQ=255)
                                                                                         aaGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGc                                                            <  1:856649/150‑1 (MQ=255)
                                                                                           gCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTACAGACCCGCGAGTGGGAAAAAGACGGc                                                            >  1:398587/1‑148 (MQ=255)
                                                                                           gCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTACAGACCCGCGAGTGGGAAAAAGACGGc                                                            <  2:398587/148‑1 (MQ=255)
                                                                                                  ggCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAgcgc                                                   <  1:434088/150‑1 (MQ=255)
                                                                                                    ccAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTa                                                 >  1:1005052/1‑150 (MQ=255)
                                                                                                        aaGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTAcacc                                             <  2:305555/150‑1 (MQ=255)
                                                                                                            tCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATAc                                                                                                                              <  2:605859/65‑1 (MQ=255)
                                                                                                            tCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATAc                                                                                                                              >  1:605859/1‑65 (MQ=255)
                                                                                                                       gAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCAtcgtc                              <  1:848909/150‑1 (MQ=255)
                                                                                                                          tGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAATTTTTCTTCAACAATCGCTACACCACCGAGATCATCGTCGat                           <  2:158111/150‑1 (MQ=255)
                                                                                                                           ggCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGata                          <  1:1009510/150‑1 (MQ=255)
                                                                                                                               ccGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGc                                                                                         <  2:195683/83‑1 (MQ=255)
                                                                                                                               ccGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGc                                                                                         >  1:195683/1‑83 (MQ=255)
                                                                                                                                cGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAAc                     <  1:195924/150‑1 (MQ=255)
                                                                                                                                  tgtgtCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACgg                   >  2:1044293/1‑150 (MQ=255)
                                                                                                                                        gCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAACAGACGGCATCAAGCACGACACCACCGAGATCATCGTCGATATCAACGGCCCCag             >  1:505113/1‑149 (MQ=255)
                                                                                                                                        gCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCAt             >  1:382557/1‑150 (MQ=255)
                                                                                                                                         gatttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGa                                                                >  1:171567/4‑98 (MQ=255)
                                                                                                                                         gatttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGa                                                                <  2:171567/95‑1 (MQ=255)
                                                                                                                                           tttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGa                                                                >  1:454232/2‑96 (MQ=255)
                                                                                                                                           tttCGGCAACCTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGa                                                                <  2:454232/95‑1 (MQ=255)
                                                                                                                                           gTTCGGCAAGGTTGCCGAAATCGCTGGAGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAATTATTTGGCAACGAGCGCTCCACCACCGAGATCATCGTCGATATCAACGGCACCATGCa          <  2:420118/150‑1 (MQ=255)
                                                                                                                                                  aaGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTc   <  2:1005052/150‑1 (MQ=255)
                                                                                                                                                   aGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTcg  <  2:1057602/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGTGGCCAGGATCCCGAAGTCCGCTACCTGCCCAACGGTAACGCCGTGACCAACCTGAGCCTGGCCACCAGCGAGCAGTGGACTGACAAGCAGTCGGGCCAGAAGGTCGAGCGTACCGAGTGGCACCGTGTGTCGCTGTTCGGCAAGGTTGCCGAAATCGCTGGCGAATACCTGCGCAAGGGTTCGCAGTGCTACATCGAAGGCAAGCTGCAGACCCGCGAGTGGGAAAAAGACGGCATCAAGCGCTACACCACCGAGATCATCGTCGATATCAACGGCACCATGCAGCTGCTCG  >  NC_002947/571069‑571365

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: