Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,674 T→C 14.3% L150L (TTG→CTG)  yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,6740TC14.3% 65.6 / 4.7 28L150L (TTG→CTG) yhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (12/12);  new base C (2/2);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTG  >  NC_002947/1486555‑1486805
                                                                                                                       |                                                                                                                                   
aCTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCcagca                                                                                                                     <  1:215121/136‑1 (MQ=255)
              ccGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGcga                                                                                                       >  2:350972/1‑134 (MQ=255)
              ccGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGCCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGcga                                                                                                       >  2:294619/1‑134 (MQ=255)
                               gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGc                                                                                      <  2:244470/136‑1 (MQ=255)
                               gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGc                                                                                      <  1:221703/136‑1 (MQ=255)
                                            aaGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGa                                                                         <  2:253668/136‑1 (MQ=255)
                                                     gAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGt                                                                >  2:55615/1‑136 (MQ=255)
                                                      aaaTTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGtt                                                               >  2:496324/1‑136 (MQ=255)
                                                      aaaTTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGtt                                                               >  2:356319/1‑136 (MQ=255)
                                                                              gtgtGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGa                                       >  1:27457/1‑136 (MQ=255)
                                                                                 tgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGAccc                                    >  2:362665/1‑136 (MQ=255)
                                                                                 tgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGAccc                                    >  1:231284/1‑136 (MQ=255)
                                                                                  gCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCg                                   <  1:356319/136‑1 (MQ=255)
                                                                                   cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCgg                                  >  1:138310/1‑136 (MQ=255)
                                                                                   cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCgg                                  >  1:123704/1‑136 (MQ=255)
                                                                                    gtaTTCCCGATCAGGCCGATAAATATCCAAGCCAGTTGTCCGGCGGGCAGCAGCAGCGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGa                                 >  1:264948/3‑136 (MQ=12)
                                                                                    ggATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGTGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGa                                 >  2:486887/1‑136 (MQ=255)
                                                                                            gaacAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCa                         <  2:287990/133‑1 (MQ=9)
                                                                                            caacAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCa                         <  1:416217/133‑1 (MQ=9)
                                                                                                   cacacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGt                  <  1:496324/136‑1 (MQ=255)
                                                                                                   cacacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGt                  <  1:263422/136‑1 (MQ=255)
                                                                                                   cacacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGt                  <  1:262922/136‑1 (MQ=255)
                                                                                                     cacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAg                       <  2:1012/129‑1 (MQ=255)
                                                                                                     cacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAg                       >  1:1012/1‑129 (MQ=255)
                                                                                                        aaGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCg             <  1:55615/136‑1 (MQ=255)
                                                                                                               cTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAt      <  2:435037/136‑1 (MQ=255)
                                                                                                                 agCCAGTTGTCCGGCGGGCAGCAGCAGCGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGAtgg    <  2:399601/135‑1 (MQ=12)
                                                                                                                   ccAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAtggtg  >  1:120438/1‑136 (MQ=255)
                                                                                                                       |                                                                                                                                   
ACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTG  >  NC_002947/1486555‑1486805

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