Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,033,379 T→G 12.5% L287R (CTG→CGG)  pgi‑I → glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,033,3790TG12.5% 97.6 / 7.4 40L287R (CTG→CGG) pgi‑Iglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (21/14);  new base G (3/2);  total (24/16)
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.38e-01

TACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGAC  >  NC_002947/2033240‑2033518
                                                                                                                                           |                                                                                                                                           
tACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTCGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCgg                                                                                                                                    >  1:186923/1‑149 (MQ=21)
                ccGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAg                                                                                                                    >  2:308364/1‑149 (MQ=34)
                     tCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAAct                                                                                                               >  3:42493/1‑149 (MQ=34)
                                   cGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGcc                                                                                                  >  1:194045/1‑148 (MQ=37)
                                        ctgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGCGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTacac                                                                                             >  6:156063/1‑148 (MQ=37)
                                         tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCa                                                                                           >  1:335168/1‑149 (MQ=37)
                                                ttCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGAc                                                                                      >  4:366553/1‑147 (MQ=37)
                                                ttCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGAc                                                                                      >  6:192625/1‑147 (MQ=37)
                                                  cGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGTGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAg                                                                                   <  1:49029/148‑1 (MQ=37)
                                                         cgcgAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTc                                                                             >  2:270690/1‑147 (MQ=37)
                                                                 aaaCATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCAc                                                                   >  5:296896/1‑149 (MQ=37)
                                                                   aCATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGcc                                                                                                                             >  5:158653/1‑89 (MQ=38)
                                                                              aTGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                       >  7:10467/1‑148 (MQ=34)
                                                                                  gggACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATg                                                  >  6:174443/1‑149 (MQ=34)
                                                                                                 ggCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGcc                                   <  8:10467/149‑1 (MQ=32)
                                                                                                    gcTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctg                                >  2:299702/1‑149 (MQ=32)
                                                                                                            cTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgt                        <  6:296896/149‑1 (MQ=32)
                                                                                                              gtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                      >  1:145271/1‑149 (MQ=32)
                                                                                                               tgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                      <  5:192625/148‑1 (MQ=32)
                                                                                                                ggTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctgc                               >  3:277615/1‑138 (MQ=32)
                                                                                                                ggTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctgc                               <  4:277615/138‑1 (MQ=32)
                                                                                                                ccggggCGATCGGCAGGCCGATGGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                     <  8:315535/144‑1 (MQ=11)
                                                                                                                ccagggCGATCGGCAGGCCGATGGCCCTGGCCATCGGCACGGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGGGCTACTGGCCCTGCTGGGCGTGTgg                    >  7:315535/5‑149 (MQ=11)
                                                                                                                   gggCGATCGGCAGGCCGATGGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACAGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTAt                 >  2:229337/2‑149 (MQ=11)
                                                                                                                    ggCGATCGGCAGGCCGATGGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                >  7:200565/1‑149 (MQ=255)
                                                                                                                    ggCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                >  8:272589/1‑149 (MQ=32)
                                                                                                                    ggCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                <  3:366553/149‑1 (MQ=32)
                                                                                                                       cATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCTGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTg                                      >  5:225136/1‑124 (MQ=17)
                                                                                                                       cATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCTGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTg                                      <  6:225136/124‑1 (MQ=17)
                                                                                                                         tCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAAc           <  2:194045/149‑1 (MQ=32)
                                                                                                                             agtGCCGATGGCCGGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGc                          <  5:257722/128‑1 (MQ=11)
                                                                                                                             agtGCCGATGGCCGGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGc                          >  6:257722/3‑130 (MQ=11)
                                                                                                                              gtGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                            <  6:238041/126‑1 (MQ=11)
                                                                                                                              gtGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                            >  5:238041/2‑127 (MQ=11)
                                                                                                                               tGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg     >  4:126164/1‑149 (MQ=255)
                                                                                                                                 ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCATACCGCACCGTTCGAAAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg     <  1:56589/147‑1 (MQ=17)
                                                                                                                                 ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg     <  5:156063/147‑1 (MQ=32)
                                                                                                                                 ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg     <  1:299702/147‑1 (MQ=32)
                                                                                                                                  cGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGAc  <  1:270690/149‑1 (MQ=32)
                                                                                                                                   gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                            <  4:327449/122‑1 (MQ=32)
                                                                                                                                           |                                                                                                                                           
TACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGAC  >  NC_002947/2033240‑2033518

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