Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,372,120 T→C 10.8% intergenic (‑48/‑71) csoR ← / → PP_2970 copper utilization transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,372,1200TC10.8% 91.7 / 3.2 37intergenic (‑48/‑71)csoR/PP_2970copper utilization transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/15);  new base C (2/2);  total (20/17)
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.00e-01

TCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGC  >  NC_002947/3371973‑3372263
                                                                                                                                                   |                                                                                                                                               
tCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTg                                                                                                                                                <  7:99909/149‑1 (MQ=255)
       gACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGa                                                                                                                                         >  4:23430/1‑149 (MQ=255)
       gACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGa                                                                                                                                         >  3:45265/1‑149 (MQ=255)
           aTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATg                                                                                                                                     <  1:43665/149‑1 (MQ=255)
              ctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGt                                                                                                                                   <  2:33170/148‑1 (MQ=255)
               tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGataggat                                                                                                                                   >  1:263746/1‑143 (MQ=255)
               tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGataggat                                                                                                                                   <  2:263746/147‑5 (MQ=255)
                 cgcAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCATATCCCCCCCAGGGGGataggatgc                                                                                                                                 <  3:218111/147‑7 (MQ=255)
                          cccTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTAcc                                                                                                                      <  6:118434/149‑1 (MQ=255)
                                 cACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCtt                                                                                                                  >  3:318186/1‑146 (MQ=255)
                                                  atgatATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAg                                                                                              >  2:73957/1‑149 (MQ=255)
                                                      tatCGCTATGACTTTGGTGCTGATGGCAGTGTTCGTGATCGCTCATGCCTTGAGTCTCATTCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCa                                                                                           >  1:304173/1‑148 (MQ=255)
                                                           cTATGACTTTGGTTCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                     <  3:107844/149‑1 (MQ=255)
                                                           cTATGACTTTGGTGCTGATGGCCGTGTTCGTGAACGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACAAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                     >  7:66138/1‑149 (MQ=255)
                                                                              ggCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTa                                                                                                                        <  2:176012/95‑1 (MQ=255)
                                                                              ggCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTa                                                                                                                        >  1:176012/1‑95 (MQ=255)
                                                                               gCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAacac                                                                              >  2:98285/1‑136 (MQ=255)
                                                                               gCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAacac                                                                              <  1:98285/136‑1 (MQ=255)
                                                                                  gtgtTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAAcc                                                              >  8:168624/1‑149 (MQ=255)
                                                                                         tGATCGCTCATGCCTTGAGTATCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAg                                                       <  1:73957/149‑1 (MQ=255)
                                                                                                  aTGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGt                                              >  8:12206/1‑149 (MQ=255)
                                                                                                    gCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTc                                                                     <  5:141183/124‑1 (MQ=255)
                                                                                                    gCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTc                                                                     >  6:141183/1‑124 (MQ=255)
                                                                                                          aGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCACAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGAttttt                                       <  4:45265/148‑1 (MQ=255)
                                                                                                            tctcATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGcc                                    >  7:118405/1‑149 (MQ=255)
                                                                                                            tctcATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGcc                                    >  5:157972/1‑149 (MQ=255)
                                                                                                             ctcATGCATCGATGAATTACGCTCGCATCCTATCCCCTGGGGGGGATATGGTCAAACCAAATACCCCTACAGAACAACCATCCGCTTAGGCATCAAAACaaa                                                                                 >  6:282253/1‑102 (MQ=255)
                                                                                                                aTGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGa                                            <  6:123446/137‑1 (MQ=255)
                                                                                                                aTGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGa                                            >  5:123446/1‑137 (MQ=255)
                                                                                                                  gCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACAGGACATGATCAGCAACCCCATCACCAATAAGGTGATTTTTGCCCTTGc                               >  4:293752/1‑148 (MQ=255)
                                                                                                                     tCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTg                           <  6:157972/149‑1 (MQ=255)
                                                                                                                             ttACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCtt                   <  7:12206/149‑1 (MQ=255)
                                                                                                                                     gCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCACCCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCt           >  4:39165/1‑149 (MQ=255)
                                                                                                                                       accatATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTgg         >  5:307895/5‑149 (MQ=255)
                                                                                                                                           tATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCg     >  5:283986/1‑149 (MQ=255)
                                                                                                                                             tCCCCCTGGGGGGTATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCTCTTATGCATCAAAATAAACACGACATGCTCAGCAATCTCATTAGCACTAAGGTGATTTTAGCTCTTGCTCTGGTCATTTTGATAGTCTGGGCCgtg   <  8:48282/149‑3 (MQ=255)
                                                                                                                                              cccccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGc  <  7:168624/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGC  >  NC_002947/3371973‑3372263

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: