Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,513,964 C→T 8.0% A142A (GCG→GCA)  gph ← phosphoglycolate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,513,9640CT 8.0% 146.6 / 2.9 50A142A (GCG→GCA) gphphosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base T (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

GACAAACAGCATCTGTTGTGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGA  >  NC_000913/3513820‑3514111
                                                                                                                                                |                                                                                                                                                   
gACAAACAGCATCTGTTGTGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGcgac                                                                                                                                                >  2:87617/1‑150 (MQ=255)
  cAAACAGCATCTGTTGTGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGcgacga                                                                                                                                              >  1:164565/1‑150 (MQ=255)
      cAGCATCTGTTGTGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACgg                                                                                                                                          >  2:1056303/1‑150 (MQ=255)
          aTCTGTTGTGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGGCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTc                                                                                                                                      <  1:1240674/150‑1 (MQ=255)
            cTGTTGTGGGGCAATTCCCATCCGCTCAGCCA‑CAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGt                                                                                                                                   >  2:1236458/1‑150 (MQ=255)
             tgttgtGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGt                                                                                                                                   <  2:1181837/150‑1 (MQ=255)
              gttgtGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGtt                                                                                                                                  <  1:1311474/150‑1 (MQ=255)
                tgtgGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTttg                                                                                                                                >  1:794438/1‑150 (MQ=255)
                   ggggCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTttgttg                                                                                                                             <  1:1056303/150‑1 (MQ=255)
                           tCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGCCCAg                                                                                                                     >  1:943803/1‑150 (MQ=255)
                                  cGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGc                                                                                                              >  2:715805/1‑150 (MQ=255)
                                   gCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCg                                                                                                             >  1:230406/1‑150 (MQ=255)
                                     tCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGc                                                                                                           <  2:304613/150‑1 (MQ=255)
                                      cAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGTGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCa                                                                                                          >  1:36698/1‑150 (MQ=255)
                                      cAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGTGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCa                                                                                                          >  1:584503/1‑150 (MQ=255)
                                        gCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGTTATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAgg                                                                                                        >  2:824956/1‑150 (MQ=255)
                                        gCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAgg                                                                                                        <  2:252249/150‑1 (MQ=255)
                                         ccaccaGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGc                                                                                                       <  2:1385613/150‑1 (MQ=255)
                                            ccaGTAACAGCGGGTCCGGATGCGGGTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCtt                                                                                                    <  1:9666/150‑1 (MQ=255)
                                               gTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAg                                                                                                 >  1:66336/1‑150 (MQ=255)
                                               gTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAg                                                                                                 <  1:258594/150‑1 (MQ=255)
                                                tAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCAGTCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGTGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGc                                                                                                <  2:36698/150‑1 (MQ=255)
                                                tAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGTGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGc                                                                                                <  2:584503/150‑1 (MQ=255)
                                                          gccGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAAcgcg                                                                                      <  2:339063/149‑1 (MQ=255)
                                                            cGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGcc                                                                                    >  2:565892/1‑150 (MQ=255)
                                                              gATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCa                                                                                  <  1:565892/150‑1 (MQ=255)
                                                                      tttttGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCg                                                                          >  2:304913/1‑150 (MQ=255)
                                                                           gttttgCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAAc                                                                     <  2:100748/150‑1 (MQ=255)
                                                                           gttttgCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAAc                                                                     <  2:95658/150‑1 (MQ=255)
                                                                                 cacaTCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCgg                                                               <  1:1312394/150‑1 (MQ=255)
                                                                                   catcatCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGa                                                             >  2:466428/1‑150 (MQ=255)
                                                                                   catcatCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGa                                                             >  2:243472/1‑150 (MQ=255)
                                                                                     tcatcaCCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAAc                                                           >  1:963573/1‑150 (MQ=255)
                                                                                            caccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACg                                                    <  1:316206/150‑1 (MQ=255)
                                                                                             accaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGt                                                   >  1:7732/1‑150 (MQ=255)
                                                                                               caATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTcc                                                 <  1:87617/150‑1 (MQ=255)
                                                                                                aaTCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAACCGTccc                                                >  2:177383/1‑150 (MQ=255)
                                                                                                 aTCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTcccc                                               <  2:230406/150‑1 (MQ=255)
                                                                                                   caacaagCTGAAGTATTTGGATATATCTAAGGCTTCGATCAGCGGCGAGAACAACGGAGACGTTTTATTAATGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCctc                                             >  2:868407/7‑150 (MQ=255)
                                                                                                     ccacGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCtt                                           >  1:103838/1‑150 (MQ=255)
                                                                                                        cGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCgg                                        <  2:164565/150‑1 (MQ=255)
                                                                                                        cGCTGAACAATTTGGCGATATCTAAGACTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTTACCAGGCCTAGCGCCAGGCCTTTACCCTGCAACCCGCCCAACCTATCGGCAACGTGCGCCAACAAAAACGTCCCCTCTTctg                                        >  2:331055/1‑148 (MQ=255)
                                                                                                          cTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCa                                      <  1:715805/150‑1 (MQ=255)
                                                                                                              aGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCt                                  >  1:896704/1‑150 (MQ=255)
                                                                                                                 aTTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGc                               >  1:1113092/1‑150 (MQ=255)
                                                                                                                 aTTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGc                               <  2:117/150‑1 (MQ=255)
                                                                                                                       cGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCAtagt                         >  1:165105/1‑150 (MQ=255)
                                                                                                                          tatCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGc                      <  1:1008798/150‑1 (MQ=255)
                                                                                                                              tAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATc                  <  2:162269/150‑1 (MQ=255)
                                                                                                                               aaGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCg                 >  2:166813/1‑150 (MQ=255)
                                                                                                                                              ggcgCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGa  >  1:473108/1‑150 (MQ=255)
                                                                                                                                              ggcgCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGa  >  1:212045/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GACAAACAGCATCTGTTGTGGGGCAATTCCCATCCGCTCAGCCACCAGTAACAGCGGGTCCGGATGCGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGA  >  NC_000913/3513820‑3514111

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: ‑