Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,027,278 Δ1 bp 8.7% coding (588/993 nt) ocd → ornithine cyclodeaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,027,2780G. 8.7% 138.8 / 7.6 46coding (588/993 nt)ocdornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base . (2/2);  total (22/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 = 3.77e-01

AAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAT  >  NC_002947/5027135‑5027425
                                                                                                                                                |                                                                                                                                                    
aaaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAgcgc                                                                                                                                                 >  3:34557/1‑149 (MQ=255)
aaaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAgc                                                                                                                                                   <  8:195259/147‑1 (MQ=255)
aaaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAgc                                                                                                                                                   <  4:32772/147‑1 (MQ=255)
         cttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTg                                                                                                                                        >  8:65868/1‑148 (MQ=255)
               gtgCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCAGCGCTCGACATTGcttcat                                                                                                                                  <  8:270287/148‑6 (MQ=255)
                   gggATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGc                                                                                                                             <  6:211827/149‑1 (MQ=255)
                              cGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCaa                                                                                                                  <  7:65868/149‑1 (MQ=255)
                                   tCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCaa                                                                                                              <  7:296710/148‑1 (MQ=255)
                                       gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCAc                                                                                                                           <  6:22011/131‑1 (MQ=255)
                                       gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCAc                                                                                                                           >  5:22011/1‑131 (MQ=255)
                                       gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCAGCGCTCGACATTGct                                                                                                                                      >  1:161343/1‑119 (MQ=255)
                                       gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCAGCGCTCGACATTGct                                                                                                                                      <  2:161343/120‑2 (MQ=255)
                                                 ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAAt                                                                                                                 >  3:131857/1‑131 (MQ=255)
                                                 ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAAt                                                                                                                 <  4:131857/131‑1 (MQ=255)
                                                 ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGc                                                                                               <  6:115172/149‑1 (MQ=255)
                                                   aGTTCCGCAAAAGCGCAGCTACTGGCCTTGCAGCTATCGAGCCTGTACGGCCTGGAAGCCTCTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAg                                                                                             <  2:39961/149‑1 (MQ=255)
                                                        cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCAc                                                                                         >  1:253436/1‑149 (MQ=255)
                                                        cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACt                                                                                        >  8:118314/1‑149 (MQ=255)
                                                           aaaaGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGc                                                                                     >  2:235598/1‑149 (MQ=255)
                                                                       aCTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCggaca                                                                           >  1:115352/1‑144 (MQ=255)
                                                                           gCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCTAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                       <  4:34557/147‑1 (MQ=255)
                                                                              ttGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCAc                                                                  >  3:78621/1‑149 (MQ=255)
                                                                                 cAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACgg                                                                <  2:253436/148‑1 (MQ=255)
                                                                                 cAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGc                                                               >  5:236936/1‑149 (MQ=255)
                                                                                         gAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATggggg                                                       <  7:118314/149‑1 (MQ=255)
                                                                                           gCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCa                                                     <  2:115352/149‑1 (MQ=255)
                                                                                                  cGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACt                                               <  3:292268/148‑1 (MQ=255)
                                                                                                  cGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTg                                              <  6:236936/149‑1 (MQ=255)
                                                                                                    gCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                            <  1:235598/149‑1 (MQ=255)
                                                                                                     ccTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                            >  5:44158/1‑148 (MQ=255)
                                                                                                     ccTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                            <  6:44158/148‑1 (MQ=255)
                                                                                                       tGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCACGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAaca                                          >  2:225394/1‑148 (MQ=255)
                                                                                                         gAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAAcaca                                        >  5:45488/1‑148 (MQ=255)
                                                                                                                   ccGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGt                              <  4:102571/148‑1 (MQ=255)
                                                                                                                   ccGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAAGCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTCAACACAATAACGCGGtt                             >  3:230393/1‑149 (MQ=255)
                                                                                                                          agTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCCCCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGc                                                      >  8:255004/1‑118 (MQ=255)
                                                                                                                          agTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGc                                                      <  7:255004/118‑1 (MQ=255)
                                                                                                                              ggtcgATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                       >  6:110002/1‑97 (MQ=255)
                                                                                                                              ggtcgATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                       <  5:110002/97‑1 (MQ=255)
                                                                                                                                    tGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATag             <  5:140158/149‑1 (MQ=255)
                                                                                                                                      aaGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGCTTGCTCCTGAGGTAATag             >  5:216231/1‑147 (MQ=255)
                                                                                                                                       tgCAATGTCAGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCACATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag          >  2:210797/2‑149 (MQ=255)
                                                                                                                                        gCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag          >  8:290825/1‑148 (MQ=255)
                                                                                                                                          aaTGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGcagcag       <  7:290825/149‑1 (MQ=255)
                                                                                                                                            tGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTc     >  5:197722/1‑149 (MQ=255)
                                                                                                                                               cGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAt  >  4:168918/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCTGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAT  >  NC_002947/5027135‑5027425

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: