Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,968,703 C→A 10.8% M53I (ATG→ATT)  gph ← phosphoglycolate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,968,7030CA10.8% 112.6 / 2.1 38M53I (ATG→ATT) gphphosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base A (2/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTCTCGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTC  >  NC_002947/1968557‑1968851
                                                                                                                                                  |                                                                                                                                                    
ccTTTCTCGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATg                                                                                                                                                   <  1:308830/150‑1 (MQ=255)
 cTTTCTCGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACGAtcc                                                                                                                                                  >  2:310049/1‑148 (MQ=255)
 cTTTCTCGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATgg                                                                                                                                                  >  1:380533/1‑150 (MQ=255)
            tCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCa                                                                                                                                       >  1:322766/1‑150 (MQ=255)
                 cagcagCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatggtt                                                                                                                                  <  1:310049/150‑6 (MQ=255)
                        tCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACg                                                                                                                           >  1:230662/1‑150 (MQ=255)
                            ccATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCATCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCAc                                                                                                                       <  2:372737/150‑1 (MQ=255)
                                gCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGAt                                                                                                                   <  2:380533/150‑1 (MQ=255)
                                 cccTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGCGATCACGCCACGGATc                                                                                                                  >  1:81483/1‑150 (MQ=255)
                                   cTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACCACAACCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAg                                                                                                                <  2:323702/150‑1 (MQ=255)
                                         cAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTgtc                                                                                                          <  1:289447/150‑1 (MQ=255)
                                                  gtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACAATGGCGCGGGCGCCACCCGAGCTCACGCCACGGATCAGGTTGTCGTCAAcggc                                                                                                 <  1:142607/150‑1 (MQ=255)
                                                  gtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACAATGGCGCGGGCGCCACCCGAGATTACCCCACGGATCAGGTTGTCGTCAAcggc                                                                                                 <  1:142654/150‑1 (MQ=255)
                                                  gtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAAAATTGAAGCGTTCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggc                                                                                                 >  2:142607/1‑150 (MQ=255)
                                                  gtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAAAAATGGCACGTGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggc                                                                                                 >  2:142654/1‑150 (MQ=255)
                                                   tgCACCGCGCAGTCGCGCGGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCTCAGGATGAGGGTGTAGTCAAcgggc                                                                                                >  2:11878/1‑148 (MQ=255)
                                                        cgcgCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCAccc                                                                                                                                    <  1:501110/109‑1 (MQ=255)
                                                        cgcgCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCAccc                                                                                                                                    >  2:501110/1‑109 (MQ=255)
                                                                 gcgcTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTc                                                                                  >  2:289804/1‑150 (MQ=255)
                                                                        tAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGc                                                                           <  2:322766/150‑1 (MQ=255)
                                                                              tCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGcc                                                                     >  2:93100/1‑150 (MQ=255)
                                                                                         cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACTCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAgg                                                                                                               <  1:317127/97‑1 (MQ=255)
                                                                                         cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAgg                                                                                                               >  2:317127/1‑97 (MQ=255)
                                                                                               cAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGaa                                                    <  1:93100/150‑1 (MQ=255)
                                                                                                    ccTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCg                                               <  1:15998/150‑1 (MQ=255)
                                                                                                     cTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCgg                                              >  2:104506/1‑150 (MQ=255)
                                                                                                        gAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGAt                                                       >  2:313907/1‑138 (MQ=255)
                                                                                                        gAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGAt                                                       <  1:313907/138‑1 (MQ=255)
                                                                                                             cGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCgtgt                                      <  2:35692/150‑1 (MQ=255)
                                                                                                                     tCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAg                              <  2:230662/150‑1 (MQ=255)
                                                                                                                     tCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAg                              <  1:509146/150‑1 (MQ=255)
                                                                                                                     tCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAg                              >  2:509146/1‑150 (MQ=255)
                                                                                                                              aTGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTcc                     >  1:463471/1‑150 (MQ=255)
                                                                                                                                   gAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTc                >  2:249582/1‑150 (MQ=255)
                                                                                                                                     aCGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGa              >  2:315173/1‑150 (MQ=255)
                                                                                                                                        cAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCGTGGCACATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAaga           <  1:11878/150‑1 (MQ=255)
                                                                                                                                        cAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAaga           <  1:104506/150‑1 (MQ=255)
                                                                                                                                           tGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTa        <  1:289804/150‑1 (MQ=255)
                                                                                                                                                 ccATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTc  <  2:81483/150‑1 (MQ=255)
                                                                                                                                                 ccATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCACGTTCGGCGAGCATGGCCTGGCAAATGCCGATCAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTTCATGTCGAAGAGTActcctc  <  2:176639/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCTTTCTCGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTC  >  NC_002947/1968557‑1968851

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