Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,238,748 A→G 6.9% V93A (GTG→GCG)  argF ← ornithine carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,238,7480AG 6.9% 166.2 / 2.6 58V93A (GTG→GCG) argFornithine carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (2/2);  total (29/29)
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.77e-01

CTGCATGTCGGCCAGCAGTTGGCACGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACTTTTCGAAGATCATGCC  >  NC_002947/1238606‑1238890
                                                                                                                                               |                                                                                                                                               
cTGCATGTCGGCCAGCAGTTGGCACGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGCCTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACAGCACCCATGCTCGACAGCACGATgg                                                                                                                                           >  8:301774/1‑149 (MQ=255)
      gTCGGCCAGCAGTTGGCACGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGt                                                                                                                                    >  1:209960/1‑149 (MQ=255)
         ggCCAGCAGTTGGCACGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCa                                                                                                                                  <  1:237279/148‑1 (MQ=255)
                 gTTGGCACGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGGCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGc                                                                                                                         <  6:99623/149‑1 (MQ=255)
                       aCGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGt                                                                                                                                    <  5:34404/132‑1 (MQ=255)
                       aCGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGt                                                                                                                                    >  6:34404/1‑132 (MQ=255)
                         gggTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCcgcg                                                                                                                 <  8:64406/149‑1 (MQ=255)
                          ggTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCgg                                                                                                                <  4:256831/149‑1 (MQ=255)
                             gCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATg                                                                                                                                            >  2:186188/1‑119 (MQ=255)
                             gCGATTCGGCGGACATGCCGTTGATCAAGGGGACACGCGACTTGGCGGCGAAATGGGTCAGGGTGCTGTGGGGGTGGGTCCGGATCATCAACCTCATCAACCATGCTCGACAGCACGATg                                                                                                                                            <  1:186188/120‑1 (MQ=255)
                                 ttCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACTTCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCt                                                                                                         <  7:301774/149‑1 (MQ=255)
                                    gtcgGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTTGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTgg                                                                                                       <  6:218653/148‑1 (MQ=255)
                                    gtcgGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGCTCACCAATTGGCTCGCCGCGGCCCAGCTg                                                                                                        <  1:208857/148‑1 (MQ=255)
                                    gtcgGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTgg                                                                                                       <  2:209960/148‑1 (MQ=255)
                                     tcgGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCGCCATCGTGCTGTcga                                                                                                                                 >  7:256094/1‑119 (MQ=255)
                                     tcgGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCGCCATCGTGCTGTcga                                                                                                                                 <  8:256094/121‑3 (MQ=255)
                                             gCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCggg                                                                                             >  5:300873/1‑149 (MQ=255)
                                                   gATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGAc                                                                                        <  6:300873/148‑1 (MQ=255)
                                                    aTCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACa                                                                                       <  5:17315/148‑1 (MQ=255)
                                                     tCACGGGCACACGCGACTTCGCGGCGAACTCGGTCAGGCTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATCCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGCCCCAGCTGGGTGTCCCGCGGCGACAg                                                                                      >  5:84351/1‑148 (MQ=255)
                                                      cACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAgg                                                                                     >  3:241983/1‑148 (MQ=255)
                                                       aCGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGa                                                                                    <  4:241983/148‑1 (MQ=255)
                                                        cGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAgg                                                                                     >  8:260482/1‑146 (MQ=255)
                                                               acGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCgggg                                                                                            <  4:197312/132‑1 (MQ=255)
                                                               acGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCgggg                                                                                            >  3:197312/1‑132 (MQ=255)
                                                                cgcgACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTcgccgc                                                                                                                  <  4:38660/109‑1 (MQ=255)
                                                                cgcgACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTcgccgc                                                                                                                  >  3:38660/1‑109 (MQ=255)
                                                                cgcgACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACGGGAAGAtggcctg                                                                          >  6:287660/1‑149 (MQ=255)
                                                                                aCTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGAt                                                           >  4:254730/1‑148 (MQ=255)
                                                                                    ggTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACCGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGc                                                      >  6:116067/1‑149 (MQ=255)
                                                                                     gTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGccc                                                                                                             >  7:206644/1‑93 (MQ=255)
                                                                                     gTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGccc                                                                                                             <  8:206644/93‑1 (MQ=255)
                                                                                      tCAGGGTGCTGTGGGCGTGGGTGCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGAGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGc                                                      <  6:84351/147‑1 (MQ=255)
                                                                                       cAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCgg                                                   >  5:191911/1‑149 (MQ=255)
                                                                                          ggTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCAt                                                        >  3:326027/1‑141 (MQ=255)
                                                                                          ggTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCAt                                                        <  4:326027/141‑1 (MQ=255)
                                                                                              cTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGa                                             <  7:260482/148‑1 (MQ=255)
                                                                                                 tgGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGa                                         >  7:252728/1‑149 (MQ=255)
                                                                                                   ggCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGac                                        <  3:263229/148‑1 (MQ=255)
                                                                                                              ggATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACacgggtacgg                            <  4:93102/149‑1 (MQ=255)
                                                                                                               gATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACacgggtacggg                           >  5:335913/1‑149 (MQ=255)
                                                                                                                atcatcACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTcgcg                                                                                               >  3:55243/1‑80 (MQ=255)
                                                                                                                   atcaCCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGa                       <  6:335913/149‑1 (MQ=255)
                                                                                                                    tcaCCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAg                      <  3:254730/149‑1 (MQ=255)
                                                                                                                           caTCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACttt                >  1:185720/1‑149 (MQ=255)
                                                                                                                            aTCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACtttt               <  8:252728/149‑1 (MQ=255)
                                                                                                                                 ccATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGccc                                                                                                             >  5:357587/1‑50 (MQ=255)
                                                                                                                                 ccATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGccc                                                                                                             <  6:357587/50‑1 (MQ=255)
                                                                                                                                  cATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACTTTTCg             >  1:114037/1‑145 (MQ=255)
                                                                                                                                  cATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACTTTTCg             >  3:162080/1‑145 (MQ=255)
                                                                                                                                       tCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACTTTTCGAAGATCATg    >  2:89343/1‑149 (MQ=255)
                                                                                                                                       tCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACGTTTCGATGGTCATg    >  8:130752/1‑149 (MQ=255)
                                                                                                                                        tgACAGCGCCATCGTGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACacgggt                                <  1:243625/119‑1 (MQ=255)
                                                                                                                                        tgACAGCGCCATCGTGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACacgggt                                >  2:243625/2‑120 (MQ=255)
                                                                                                                                         gACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAAGACTTTTCGAAGATCATGc   <  1:161093/148‑1 (MQ=255)
                                                                                                                                          aCAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACTTTTCGAAGATCATGcc  >  6:118385/1‑148 (MQ=255)
                                                                                                                                             gCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATc                                                          <  6:22206/89‑1 (MQ=255)
                                                                                                                                             gCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATc                                                          >  5:22206/1‑89 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CTGCATGTCGGCCAGCAGTTGGCACGGGTGCGATTCGTCGGACAGGCCGTTGATCACGGGCACACGCGACTTGGCGGCGAACTCGGTCAGGGTGCTGTGGGCGTGGGTCCGGATCATCACCACATCAACCATGCTCGACAGCACGATGGCGCTGTCACCAATTGGCTCGCCGCGGCCCAGCTGGGTGTCGCGGGGCGACAGGAAGATGGCCTGGCCGCCGAGCTGGATCATGCCGGCCTCGAACGACACACGGGTACGGGTCGAGGACTTTTCGAAGATCATGCC  >  NC_002947/1238606‑1238890

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: