Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,974,790 G→A 7.7% R267H (CGC→CAC)  serC → 3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,974,7900GA 7.7% 150.9 / 2.0 52R267H (CGC→CAC) serC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base A (2/2);  total (27/25)
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.97e-01

CAGCTGCCCGACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCTCGCTG  >  NC_002947/1974647‑1974933
                                                                                                                                                |                                                                                                                                               
cagcTGCCCGACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAAcc                                                                                                                                            <  4:142767/149‑1 (MQ=255)
  gcTGCCCGACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGt                                                                                                                                          >  5:168555/1‑149 (MQ=255)
  gcTGCCCGACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGt                                                                                                                                          <  4:132390/149‑1 (MQ=255)
      cccGACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGc                                                                                                                                        <  6:285386/147‑1 (MQ=255)
         gACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCaag                                                                                                                                    >  8:5626/1‑148 (MQ=255)
            cATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCaagaaag                                                                                                                                <  1:245542/149‑1 (MQ=255)
             aTGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCaagaaaga                                                                                                                               <  7:5626/149‑1 (MQ=255)
             aTGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCaagaaaga                                                                                                                               <  6:168555/149‑1 (MQ=255)
                cTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGCAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGAc                                                                                                                              <  1:80263/148‑1 (MQ=255)
                cTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGAcc                                                                                                                             <  7:16118/148‑1 (MQ=255)
                   gACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGc                                                                                                                           <  2:118251/147‑1 (MQ=255)
                           ggTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCt                                                                                                                 <  2:99262/149‑1 (MQ=255)
                           ggTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCt                                                                                                                 <  8:29786/149‑1 (MQ=255)
                                  gCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCACAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGa                                                                                                          >  8:64817/1‑149 (MQ=255)
                                  gCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCACAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGa                                                                                                          >  8:69022/1‑149 (MQ=255)
                                            gCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCACAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAg                                                                                                <  7:64817/149‑1 (MQ=255)
                                            gCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCACAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAg                                                                                                <  7:69022/149‑1 (MQ=255)
                                             cTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGt                                                                                               >  8:282079/1‑149 (MQ=255)
                                                gATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTc                                                                                             >  6:106994/1‑148 (MQ=255)
                                                gATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCt                                                                                            <  7:282079/149‑1 (MQ=255)
                                                         cacCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGCCTGAAAGAGCAGGGCGCTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAACAAAGCCCGCCTTTACGGCTTCATCGACCCCAGCGAGTTCAACACCAAcc                                                                                   >  5:56180/1‑149 (MQ=255)
                                                            cccGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGa                                                                                >  5:242392/1‑149 (MQ=255)
                                                                           cTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATGGTCCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGccc                                                                 >  1:102665/1‑149 (MQ=255)
                                                                           cTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCa                                                                             <  1:237132/137‑1 (MQ=255)
                                                                           cTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCa                                                                             >  2:237132/1‑137 (MQ=255)
                                                                              gTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGt                                                               <  5:106994/148‑1 (MQ=255)
                                                                                 ccTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATTGAGCAGCGCAACAGTGCCAAGAAAGATTTTTTTTTTTTTTTTATGGAGAGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGt                                                           <  6:89987/149‑1 (MQ=255)
                                                                                           cTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGt                                                 >  6:272815/1‑149 (MQ=255)
                                                                                            tGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTg                                                >  1:106788/1‑149 (MQ=255)
                                                                                            tGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTg                                                >  8:41878/1‑149 (MQ=255)
                                                                                                cTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGt                                            >  5:153654/1‑149 (MQ=255)
                                                                                                      gTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGcc                                      >  1:140493/1‑149 (MQ=255)
                                                                                                      gTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCATCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGcc                                      >  4:217677/1‑149 (MQ=255)
                                                                                                         gCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTgg                                   >  4:216124/1‑149 (MQ=255)
                                                                                                          cTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGc                                  >  2:220454/1‑149 (MQ=255)
                                                                                                              aaaGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACGGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGa                               >  5:169723/1‑149 (MQ=255)
                                                                                                              aaaGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTATACACAAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCAGTTCCGCCTGGCTg                                >  5:53974/1‑148 (MQ=255)
                                                                                                                  agCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCATTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAg                           >  1:235318/1‑149 (MQ=255)
                                                                                                                    cAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAgcg                         >  1:220072/1‑149 (MQ=255)
                                                                                                                      gggCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACATTTTTTAATTTTTTATGGAGAGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGcc                       <  6:153654/149‑1 (MQ=255)
                                                                                                                        gCGGTGTCCAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGATTGTTTTTTTTGTTTCATGGAGAGCAGCGAGATCTACACCTACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTg                     <  4:41522/149‑1 (MQ=255)
                                                                                                                         cGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTgg                    <  1:220454/149‑1 (MQ=255)
                                                                                                                                gAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGACCGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACACGGc             >  5:67339/1‑149 (MQ=255)
                                                                                                                                gAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCATCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAgg              <  3:217677/148‑1 (MQ=255)
                                                                                                                                gAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACAGTTTTTATTTTTTGATGGAGAGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAgg              <  2:106788/148‑1 (MQ=255)
                                                                                                                                 aGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGcc            <  7:41878/149‑1 (MQ=255)
                                                                                                                                     cATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCAGCCTGGCTGACGAGCGCCCGGACAAGGCCTTcc        >  7:116318/1‑149 (MQ=255)
                                                                                                                                      aTGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCt       >  3:270133/1‑149 (MQ=255)
                                                                                                                                         gAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCTCGc    >  3:107353/1‑149 (MQ=255)
                                                                                                                                          agcagcGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACTGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCTCGCt   <  6:169723/149‑1 (MQ=255)
                                                                                                                                          agcagcGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCTCGCt   <  5:272815/149‑1 (MQ=255)
                                                                                                                                           gcagcGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCTCGCTg  >  6:150866/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CAGCTGCCCGACCATGCTCGACTACAAGGTTTCGGCTGACAACGGCTCGATGTACAACACCCCGGCCACCTACTCCTGGTACCTCTCTGGCCTGGTCTTCGAGTGGCTGAAAGAGCAGGGCGGTGTCGAGGCCATGGAGCAGCGCAACCGTGCCAAGAAAGACCGCCTTTACGGCTTCATCGACCGCAGCGAGTTCTACACCAACCCGATCAGCGTCAACGCCCGTTCGTGGATGAACGTGCCGTTCCGCCTGGCTGACGAGCGCCTGGACAAGGCCTTCCTCGCTG  >  NC_002947/1974647‑1974933

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