Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,769 C→T 15.1% D181D (GAC→GAT 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,7690CT15.1% 85.0 / 6.2 33D181D (GAC→GATyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base T (2/3);  total (14/19)
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.82e-01

TGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATG  >  NC_002947/1486636‑1486892
                                                                                                                                     |                                                                                                                           
tgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGAccc                                                                                                                           >  2:362665/1‑136 (MQ=255)
tgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGAccc                                                                                                                           >  1:231284/1‑136 (MQ=255)
 gCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCg                                                                                                                          <  1:356319/136‑1 (MQ=255)
  cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCgg                                                                                                                         >  1:123704/1‑136 (MQ=255)
  cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCgg                                                                                                                         >  1:138310/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:262922/136‑1 (MQ=255)
                  cacacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGt                                                                                                         <  1:263422/136‑1 (MQ=255)
                  cacacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGt                                                                                                         <  1:496324/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)
                                        gTCCGGCGGGCAGCAGAAGCGTGTCGCCATAGCGCGTGCCTGGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGGCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGAtggtgtctctc                                                                                   <  2:29903/136‑7 (MQ=255)
                                           aggcggcCAGCAGCAGCGTGGGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGcc                                                                                <  2:198124/135‑1 (MQ=17)
                                                  cagcagcagCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATg                                                                         <  1:267491/136‑1 (MQ=255)
                                                  cagcagcagCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATg                                                                         <  2:250628/136‑1 (MQ=255)
                                                   agcagcagcGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATgg                                                                        >  1:452903/1‑136 (MQ=255)
                                                                    aTTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCGGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGGGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCg                                                       <  2:345550/136‑1 (MQ=255)
                                                                      tGCCCGCGCGTTGTGCATGAAGCCGAAGATCATGCTGTTCGACGAGCAGACCTCGGCACTGGTCCGGGTAATGTTCTAGGAAGTGCTAGATACCATGGTGTGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTg                                                     <  2:178066/136‑1 (MQ=38)
                                                                        cccGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGAc                                                   >  1:437068/1‑136 (MQ=255)
                                                                        cccGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGACGATGGCATGACCATGCTCTGCGTGAc                                                   >  1:205139/1‑136 (MQ=255)
                                                                               gctttgCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGat                                            >  1:352492/5‑135 (MQ=12)
                                                                                  gtgCATGAAGNCCAAGATCATGCTGTTCGACGAGCCGACCTCGGNACTGGANCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCNTGCTCTGCGTGACCCACGAAATg                                         <  2:171157/136‑1 (MQ=255)
                                                                                  gtgCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATg                                         <  2:452903/136‑1 (MQ=255)
                                                                                      aTGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCt                                     >  2:118542/1‑136 (MQ=255)
                                                                                                aGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTAc                           >  2:26675/1‑136 (MQ=255)
                                                                                                  gtCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGtcagg                         >  2:466503/2‑132 (MQ=17)
                                                                                                         tGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGaa                  >  1:483496/1‑136 (MQ=255)
                                                                                                                         cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATg  <  1:145190/136‑1 (MQ=255)
                                                                                                                         cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATg  <  2:483496/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                           
TGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATG  >  NC_002947/1486636‑1486892

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: