Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,655 A→C 25.1% A143A (GCA→GCC 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,6550AC25.1% 38.0 / 11.6 20A143A (GCA→GCCyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base C (3/2);  total (13/7)
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.67e-01

CAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGT  >  NC_002947/1486523‑1486789
                                                                                                                                    |                                                                                                                                      
cAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCATGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGcacac                                                                                                                                     >  2:267491/1‑136 (MQ=255)
cAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGcacac                                                                                                                                     >  1:146578/1‑136 (MQ=255)
cAACCTGTTCCCGCAACTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGAAAGTCCGAGGAAATTTCCATTCAATACCTGGAGAGAGTGCGGAGTCAGGAGCAGGcacac                                                                                                                                     >  1:178066/1‑136 (MQ=255)
         cccGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCt                                                                                                                            >  1:253668/1‑136 (MQ=255)
                                aCTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCcagca                                                                                                     <  1:215121/136‑1 (MQ=255)
                                              ccGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGcga                                                                                       >  2:350972/1‑134 (MQ=255)
                                              ccGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGCCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGcga                                                                                       >  2:294619/1‑134 (MQ=255)
                                                               gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGc                                                                      <  1:221703/136‑1 (MQ=255)
                                                               gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGc                                                                      <  2:244470/136‑1 (MQ=255)
                                                                            aaGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGa                                                         <  2:253668/136‑1 (MQ=255)
                                                                                     gAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGt                                                >  2:55615/1‑136 (MQ=255)
                                                                                      aaaTTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGtt                                               >  2:356319/1‑136 (MQ=255)
                                                                                      aaaTTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGtt                                               >  2:496324/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)
                                                                                                                                    |                                                                                                                                      
CAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGT  >  NC_002947/1486523‑1486789

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