Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,477 A→G 8.3% Q150R (CAG→CGG)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,4770AG 8.3% 132.3 / 2.6 48Q150R (CAG→CGG) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base G (2/2);  total (25/23)
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.66e-01

ATATCAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCT  >  NC_002947/571336‑571624
                                                                                                                                             |                                                                                                                                                   
atatCAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCg                                                                                                                                              >  1:23353/1‑149 (MQ=255)
 tatCAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCg                                                                                                                                              >  2:31312/1‑148 (MQ=255)
 tatCAACGGCACAATGCAGCTGCACGGAGGTCGTCCGCAGGGCCAGAAGCAGGGCGGCGACACGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCCCCCg                                                                                                                                              >  2:94401/1‑148 (MQ=255)
  atCAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGc                                                                                                                                             >  1:242635/1‑148 (MQ=255)
   tCAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCaa                                                                                                                                           >  6:2032/1‑149 (MQ=255)
    cAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCaa                                                                                                                                           >  2:57167/1‑148 (MQ=255)
     aaCGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACa                                                                                                                                         <  3:295080/149‑1 (MQ=255)
                               tcgtcCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAAc                                                                                                               <  2:23353/149‑1 (MQ=255)
                                 gtcCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACACCCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAg                                                                                                             >  5:96774/1‑149 (MQ=255)
                                                    gggCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCcggcgg                                                                                          <  2:304653/149‑1 (MQ=255)
                                                       cggcgACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCt                                                                                        >  8:79562/1‑148 (MQ=255)
                                                           gACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACtt                                                                                   <  1:31312/149‑1 (MQ=255)
                                                            aCCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTc                                                                                  <  1:57167/149‑1 (MQ=255)
                                                              ccGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACCGCCTGCGCCGCAGCCGGCGGCTGACTTCGa                                                                                >  2:187240/1‑149 (MQ=255)
                                                                    aaCCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCGGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCtt                                                                          >  2:318304/1‑149 (MQ=255)
                                                                    aaCCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCGGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCtt                                                                          >  4:213325/1‑149 (MQ=255)
                                                                    aaCCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCtt                                                                          <  3:44163/149‑1 (MQ=255)
                                                                      ccAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTc                                                                         >  2:321619/1‑148 (MQ=255)
                                                                      ccAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCACCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGACTGACTTCGACAGCTTc                                                                         >  7:13134/1‑148 (MQ=255)
                                                                        aaGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGAt                                                                      <  1:321619/149‑1 (MQ=255)
                                                                             ggcggcAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGa                                                                 >  2:202931/1‑149 (MQ=255)
                                                                               cggcAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGata                                                               <  4:313795/149‑1 (MQ=255)
                                                                                     ttATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGt                                                         <  2:3080/149‑1 (MQ=255)
                                                                                          gtggtggCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGa                                                    >  4:139449/1‑149 (MQ=255)
                                                                                             gtggCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATc                                                  <  2:242635/148‑1 (MQ=255)
                                                                                                 ccAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGc                                             >  3:81194/1‑149 (MQ=255)
                                                                                                   agcagcAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCGGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCgg                                           <  3:213325/149‑1 (MQ=255)
                                                                                                   agcagcAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCGGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCgg                                           <  1:318304/149‑1 (MQ=255)
                                                                                                   agcagcAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCgg                                           >  4:95748/1‑149 (MQ=255)
                                                                                                     cagcaACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTa                                         <  5:206867/149‑1 (MQ=255)
                                                                                                       gcaACAGTACACCCCGGCACCGCCCCGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCCGCGCCCCGCGCCGCAACAGCCAGCGCCGCAGCCGGCGGCTGAGTTCGACAGCGTCGCTGCCGCTATTCCGTTCTGATCCGAGCGGTAgg                                       >  6:239517/1‑149 (MQ=255)
                                                                                                       gcaACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAgg                                       <  7:79562/149‑1 (MQ=255)
                                                                                                                 aaCCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGc                              >  6:58036/1‑148 (MQ=255)
                                                                                                                 aaCCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGc                              >  5:216616/1‑148 (MQ=255)
                                                                                                                 aaCCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGc                              >  8:260598/1‑148 (MQ=255)
                                                                                                                 aaCCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCATCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGc                              >  8:1060/1‑148 (MQ=255)
                                                                                                                    cAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGa                          <  5:58036/149‑1 (MQ=255)
                                                                                                                     aGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGAc                         <  3:139449/149‑1 (MQ=255)
                                                                                                                      ggCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGAcc                        <  7:1060/149‑1 (MQ=255)
                                                                                                                       gCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGAc                         >  5:117984/1‑147 (MQ=255)
                                                                                                                       gCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCa                       <  7:80721/149‑1 (MQ=255)
                                                                                                                              cacgccaGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCaa                            <  1:301801/137‑1 (MQ=255)
                                                                                                                              cacgccaGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCaa                            >  2:301801/1‑137 (MQ=255)
                                                                                                                               acgccaGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCa               <  6:96774/149‑1 (MQ=255)
                                                                                                                                     gcagGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAgg             >  6:19603/1‑145 (MQ=255)
                                                                                                                                     gcagGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAgg             <  5:19603/145‑1 (MQ=255)
                                                                                                                                       agGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGtt       <  2:146650/149‑1 (MQ=255)
                                                                                                                                            cAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCt  <  6:216616/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ATATCAACGGCACCATGCAGCTGCTCGGCGGTCGTCCGCAGGGCCAGCAGCAGGGCGGCGACCCGTACAACCAAGGTGGCGGCAATTATGGTGGTGGCCAGCAGCAACAGTACAACCAGGCACCGCCACGCCAGCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCT  >  NC_002947/571336‑571624

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