Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,712 C→T 17.2% R162R (CGC→CGT 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,7120CT17.2% 74.0 / 7.1 29R162R (CGC→CGTyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (16/8);  new base T (2/3);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTC  >  NC_002947/1486581‑1486835
                                                                                                                                   |                                                                                                                           
cgcAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgc                                                                                                                         >  1:470735/1‑135 (MQ=255)
     gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGc                                                                                                                    >  2:33611/1‑136 (MQ=255)
               cgcAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCa                                                                                                          >  2:54641/1‑136 (MQ=255)
                gcAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGaa                                                                                                               >  2:554320/1‑130 (MQ=255)
                gcAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGaa                                                                                                               <  1:554320/130‑1 (MQ=255)
                                  ccATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTcgacga                                                                                       >  2:276526/1‑136 (MQ=255)
                                      gCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCg                                                                                   >  2:130730/1‑136 (MQ=9)
                                              tcaaGCGGGTCCGTATTCCCGATCAGGCCGATAAATATCCAAGCCAGTTGTCCGGCGGGCAGCAGCAGCGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGc                                                                           <  1:543344/133‑1 (MQ=12)
                                               ggAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGcc                                                                          >  1:188473/1‑135 (MQ=255)
                                               ggAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCa                                                                          >  2:263783/1‑136 (MQ=255)
                                               ggAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCAAAGATCATGCTGTTCGACGAGCCGACCTCGGcc                                                                          >  1:26171/1‑135 (MQ=255)
                                               ggAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGACCAAGATCATGCTGTTCGACGAGCCGACCTCggaa                                                                          >  1:23431/1‑134 (MQ=255)
                                                    gtgtGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGc                                                                           <  2:477427/130‑1 (MQ=255)
                                                    gtgtGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGc                                                                           >  1:477427/1‑130 (MQ=255)
                                                    gtgtGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGATGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGa                                                                     >  2:268660/1‑136 (MQ=255)
                                                       tgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGAccc                                                                  >  1:108302/1‑136 (MQ=255)
                                                        gCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCg                                                                 <  1:413334/136‑1 (MQ=255)
                                                        gCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCg                                                                 <  1:437005/136‑1 (MQ=255)
                                                                    gCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAg                                                     <  1:436779/136‑1 (MQ=255)
                                                                               aGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGa                                          <  2:108302/136‑1 (MQ=255)
                                                                               aGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGa                                          <  1:536520/136‑1 (MQ=255)
                                                                                         ccAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAtggtg                                >  1:153611/1‑136 (MQ=255)
                                                                                           aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTggg                              >  1:187054/1‑136 (MQ=255)
                                                                                           aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTggg                              >  2:395842/1‑136 (MQ=255)
                                                                                           aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTggg                              >  2:529073/1‑136 (MQ=255)
                                                                                           aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTggg                              >  1:148820/1‑136 (MQ=255)
                                                                                                  gggcggcCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGcc                       <  2:32819/135‑1 (MQ=255)
                                                                                                              gcagcGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATg           <  1:274147/136‑1 (MQ=17)
                                                                                                                   gtgtCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGAt      <  2:467673/136‑1 (MQ=255)
                                                                                                                       ggCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGatc  >  1:173244/1‑136 (MQ=12)
                                                                                                                                   |                                                                                                                           
CGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTC  >  NC_002947/1486581‑1486835

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: