Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,033,374 C→G 22.0% I285M (ATC→ATG 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,3740CG22.0% 82.6 / 20.7 41I285M (ATC→ATGpgi‑Iglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (19/13);  new base G (6/3);  total (25/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 = 6.21e-01

GAAGCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGAC  >  NC_002947/2033228‑2033518
                                                                                                                                                  |                                                                                                                                                
gAAGCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGTGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGc                                                                                                                                                >  2:49029/1‑149 (MQ=21)
            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                          >  6:257722/3‑130 (MQ=11)
                                                                                                                                         agtGCCGATGGCCGGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGc                          <  5:257722/128‑1 (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     <  1:299702/147‑1 (MQ=32)
                                                                                                                                             ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg     <  5:156063/147‑1 (MQ=32)
                                                                                                                                              cGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGAc  <  1:270690/149‑1 (MQ=32)
                                                                                                                                               gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                            <  4:327449/122‑1 (MQ=32)
                                                                                                                                                  |                                                                                                                                                
GAAGCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGAC  >  NC_002947/2033228‑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: