Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,261,988 C→T 16.8% R227R (CGC→CGT yeiR → Zn‑stimulated GTPase involved in zinc homeostasis; mutants are cadmium and EDTA sensitive; Zn(2+) binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,261,9880CT16.8% 44.1 / 4.8 24R227R (CGC→CGTyeiRZn‑stimulated GTPase involved in zinc homeostasis; mutants are cadmium and EDTA sensitive; Zn(2+) binding protein
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base T (2/2);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

TCACAGTGAGCATGGGAAAGTTGACGGTCATCTTCTGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGCGTGCTGCGTATTCCCGAAG  >  CP009273/2261844‑2262136
                                                                                                                                                |                                                                                                                                                    
tCACGGTGAGCTTGGGAAAGTTGACGGTCATCTTCTGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGc                                                                                                                                                  <  2:348167/149‑1 (MQ=255)
          cATGCGAAAGTTGACGGTCATCTTCTGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCg                                                                                                                                              <  1:144763/143‑1 (MQ=255)
                   ggtGACGGTCATCTTCTGGATTTGCGGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGaa                                                                                                                                     <  1:326589/141‑1 (MQ=255)
                      gACGGTCATCTTCTGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAg                                                                                                                                  <  1:202814/143‑1 (MQ=255)
                                ttctGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCCAGGATATCAg                                                                                                                  >  2:129202/1‑149 (MQ=255)
                                 tctGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAgg                                                                                                                 <  2:340653/149‑1 (MQ=255)
                                            cgcgTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTgg                                                                                                      >  2:89188/1‑149 (MQ=255)
                                                        ctccggatttGCCCGCCAGCGCCGCGCTTTCTCATCAGCATTGCGTGAAAAAAGGGTTAGCAGCGTTACGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTc                                                                                                <  1:67372/135‑1 (MQ=255)
                                                                    ccGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGTTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACg                                                                                    >  1:288163/1‑143 (MQ=255)
                                                                    ccGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGTTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACg                                                                                    >  1:288121/1‑143 (MQ=255)
                                                                          gcgcCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTc                                                                              <  1:377499/143‑1 (MQ=255)
                                                                                 ggaTTTTCATCGGCATGTCGTGAACAAAGGTTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGTAGTCTGACCAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCa                                                                       <  1:323582/141‑1 (MQ=255)
                                                                                              cATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGgc                                                    >  2:126528/1‑149 (MQ=255)
                                                                                                     tGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGGTATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACtt                                             >  2:203325/1‑149 (MQ=255)
                                                                                                         aaaaGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCAc                                         <  2:370496/149‑1 (MQ=255)
                                                                                                         aaaaGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCAc                                         <  2:305153/149‑1 (MQ=255)
                                                                                                            aGGGTTAGCGGCGTTAAGCCTGCCGGTGCATCAACGCTGGCGTCGCAGTCTGAGCAGCGGGCAAGGATATCAGGCCTGCGGCTGGAGATTCGAGGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTg                                            <  1:56546/143‑1 (MQ=255)
                                                                                                            aGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTg                                            <  1:120392/143‑1 (MQ=255)
                                                                                                              ggTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCa                                          >  1:368327/1‑143 (MQ=255)
                                                                                                                             gCCTGCCAGAGCATCAACGTTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGc                     <  2:323856/149‑1 (MQ=255)
                                                                                                                             gCCTGCCAGAGCATCAACGTTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGc                     <  2:323812/149‑1 (MQ=255)
                                                                                                                              ccTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGCg                    >  2:353647/1‑149 (MQ=255)
                                                                                                                                 gCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGACTGGCCGCGACTTGCACCGGTGGAACGCGGCAAAGGCGtgc                 >  2:405008/1‑149 (MQ=255)
                                                                                                                                          tCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGCGtgctgc              <  1:111456/143‑1 (MQ=255)
                                                                                                                                                cTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGCCGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGGCTTGCACCGGTGGAACGCGTCAAAGGCGTGCTGCGTATTCCCGAAg  >  2:13263/1‑149 (MQ=255)
                                                                                                                                                cTGGCGTCGCAGTCGGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGCGTGCTGCGTATTCCCGAAg  <  2:7930/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TCACAGTGAGCATGGGAAAGTTGACGGTCATCTTCTGGATTTGCCGCGTCGCAATTTAGCCGAGTTGCCCGCCAGCGCCGCGCATTCTCATCAGCATGTCGTGAAAAAAGGGTTAGCAGCGTTAAGCCTGCCAGAGCATCAACGCTGGCGTCGCAGTCTGAACAGCGGGCAAGGATATCAGGCCTGCGGCTGGATATTCGACGCTGATACGGTATTCGACACCATTGGCATTCTGGAATGGGCGCGACTTGCACCGGTGGAACGCGTCAAAGGCGTGCTGCGTATTCCCGAAG  >  CP009273/2261844‑2262136

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: