Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,362,215 C→G 11.6% V49L (GTG→CTG)  tpiA ← triose phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,362,2150CG11.6% 113.8 / 7.6 43V49L (GTG→CTG) tpiAtriose phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (21/17);  new base G (2/3);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

GCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAC  >  NC_002947/5362078‑5362355
                                                                                                                                         |                                                                                                                                            
gCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatgg                                                                                                                                   <  1:223003/149‑2 (MQ=255)
        acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAc                                                                                                                                              >  1:78574/1‑130 (MQ=255)
        acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAc                                                                                                                                              <  2:78574/130‑1 (MQ=255)
        acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatggc                                                                                                                                  <  1:177705/142‑3 (MQ=255)
        acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatggc                                                                                                                                  >  2:177705/1‑140 (MQ=255)
        acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                     >  4:30011/1‑139 (MQ=255)
        acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                     <  3:30011/139‑1 (MQ=255)
            cAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTg                                                                                                                       >  4:69108/1‑149 (MQ=255)
                  gCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg                                                                                                                      >  4:325266/1‑144 (MQ=255)
                  gCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg                                                                                                                      <  3:325266/144‑1 (MQ=255)
                           ttCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                     <  8:211574/120‑1 (MQ=255)
                           ttCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                     >  7:211574/1‑120 (MQ=255)
                             cAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACtt                                                                                                        <  4:168926/147‑1 (MQ=255)
                                        tCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTc                                                                                             <  7:147533/147‑1 (MQ=255)
                                        tCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTc                                                                                             >  8:147533/1‑147 (MQ=255)
                                              gCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGAGCCGACAGTAATTTCTTTACCTGCAAGGCCAACAATCACCTGATTGATGAACAAGGCCGGTCGAAATACgg                                                                                                              >  2:154696/1‑124 (MQ=255)
                                               cAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTGCAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                       >  6:25095/1‑130 (MQ=255)
                                               cAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTGCAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                       <  5:25095/130‑1 (MQ=255)
                                                aaCTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCga                                                                                   <  6:47856/149‑1 (MQ=255)
                                                  cTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCgaga                                                                                 <  7:113503/149‑1 (MQ=255)
                                                       ccGGTCAGCGCACCCTGTTCGGGTTGTAAAGCAGAATTCTGTGCGCCGTCAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACt                                                                            >  6:75236/1‑149 (MQ=255)
                                                                   cccTGTTCGGGTTGTAAAGCAGAATTCTGTGCGCCGTCAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                <  5:75236/149‑1 (MQ=255)
                                                                    ccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCaa                                                                         >  2:46665/1‑139 (MQ=255)
                                                                    ccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCaa                                                                         <  1:46665/139‑1 (MQ=255)
                                                                    ccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                >  2:282199/1‑148 (MQ=255)
                                                                            tgTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCTGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                  <  6:249410/137‑1 (MQ=255)
                                                                            ggTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                  >  5:249410/1‑138 (MQ=255)
                                                                                    agcagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAg                                               >  5:1407/1‑149 (MQ=255)
                                                                                      cagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                           >  3:325055/1‑119 (MQ=255)
                                                                                      cagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                           <  4:325055/119‑1 (MQ=255)
                                                                                           ttCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCggg                                         >  2:113628/1‑148 (MQ=255)
                                                                                           ttCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCggg                                         >  7:20292/1‑148 (MQ=255)
                                                                                                tgCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGt                                   >  5:213485/1‑149 (MQ=255)
                                                                                                         aGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTcc                          <  8:20292/149‑1 (MQ=255)
                                                                                                          gTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa                         <  3:69108/149‑1 (MQ=255)
                                                                                                                tCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTAc                    >  3:63983/1‑148 (MQ=255)
                                                                                                                      ccTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTa              <  6:1407/148‑1 (MQ=255)
                                                                                                                       cTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTAc             >  1:34346/1‑148 (MQ=255)
                                                                                                                         gCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCa           >  7:179375/1‑148 (MQ=255)
                                                                                                                          cttgttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTc                           <  7:126703/125‑1 (MQ=255)
                                                                                                                          cttgttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTc                           >  8:126703/7‑131 (MQ=255)
                                                                                                                            aGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAgg       >  6:77060/1‑149 (MQ=255)
                                                                                                                             ggCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg                                          <  4:91083/113‑1 (MQ=255)
                                                                                                                             ggCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg                                          >  3:91083/1‑113 (MQ=255)
                                                                                                                                 atcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAc  <  6:213485/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
GCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAC  >  NC_002947/5362078‑5362355

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: