Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,027,284:1	+C	9.1%	coding (594/993 nt)	ocd →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,027,284	1	.	C	9.1%	160.2 / 10.8	44	coding (594/993 nt)	ocd	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base . (19/21);  new base C (2/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACG  >  NC_002947/5027144‑5027432                                                                                                                                               |                                                                                                                                                     cttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTg                                                                                                                                               >  8:65868/1‑148 (MQ=255)       gtgCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGcttcat                                                                                                                                         <  8:270287/148‑6 (MQ=255)           gggATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGc                                                                                                                                    <  6:211827/149‑1 (MQ=255)                      cGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCaa                                                                                                                         <  7:65868/149‑1 (MQ=255)                           tCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCaa                                                                                                                     <  7:296710/148‑1 (MQ=255)                               gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCAc                                                                                                                                  <  6:22011/131‑1 (MQ=255)                               gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCAc                                                                                                                                  >  5:22011/1‑131 (MQ=255)                               gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                             >  1:161343/1‑119 (MQ=255)                               gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                             <  2:161343/120‑2 (MQ=255)                                         ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAAt                                                                                                                        >  3:131857/1‑131 (MQ=255)                                         ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAAt                                                                                                                        <  4:131857/131‑1 (MQ=255)                                         ggAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGc                                                                                                      <  6:115172/149‑1 (MQ=255)                                           aGTTCCGCAAAAGCGCAGCTACTGGCCTTGCAGCTATCGAGCCTGTACGGCCTGGAAGCCTCTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAg                                                                                                    <  2:39961/149‑1 (MQ=255)                                                cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCAc                                                                                                >  1:253436/1‑149 (MQ=255)                                                cGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACt                                                                                               >  8:118314/1‑149 (MQ=255)                                                   aaaaGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGc                                                                                            >  2:235598/1‑149 (MQ=255)                                                               aCTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCggaca                                                                                  >  1:115352/1‑144 (MQ=255)                                                                   gCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCTAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                              <  4:34557/147‑1 (MQ=255)                                                                      ttGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCAc                                                                         >  3:78621/1‑149 (MQ=255)                                                                         cAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACgg                                                                       <  2:253436/148‑1 (MQ=255)                                                                         cAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGc                                                                      >  5:236936/1‑149 (MQ=255)                                                                                 gAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATggggg                                                              <  7:118314/149‑1 (MQ=255)                                                                                   gCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCa                                                            <  2:115352/149‑1 (MQ=255)                                                                                          cGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACt                                                      <  3:292268/148‑1 (MQ=255)                                                                                          cGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTg                                                     <  6:236936/149‑1 (MQ=255)                                                                                            gCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                                   <  1:235598/149‑1 (MQ=255)                                                                                             ccTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                                   <  6:44158/148‑1 (MQ=255)                                                                                             ccTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                                   >  5:44158/1‑148 (MQ=255)                                                                                               tGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCACGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAaca                                                 >  2:225394/1‑148 (MQ=255)                                                                                                 gAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAAcaca                                               >  5:45488/1‑148 (MQ=255)                                                                                                           cc‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGt                                     <  4:102571/148‑1 (MQ=255)                                                                                                           cc‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAAGCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTCAACACAATAACGCGGtt                                    >  3:230393/1‑149 (MQ=255)                                                                                                                  agTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCCCCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGc                                                             >  8:255004/1‑118 (MQ=255)                                                                                                                  agTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGc                                                             <  7:255004/118‑1 (MQ=255)                                                                                                                      ggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                              >  6:110002/1‑97 (MQ=255)                                                                                                                      ggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                              <  5:110002/97‑1 (MQ=255)                                                                                                                            tGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATag                    <  5:140158/149‑1 (MQ=255)                                                                                                                              aaGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGCTTGCTCCTGAGGTAATag                    >  5:216231/1‑147 (MQ=255)                                                                                                                               tgCAATGTC‑AGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCACATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag                 >  2:210797/2‑149 (MQ=255)                                                                                                                                gCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcag                 >  8:290825/1‑148 (MQ=255)                                                                                                                                  aaTGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGcagcag              <  7:290825/149‑1 (MQ=255)                                                                                                                                    tGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTc            >  5:197722/1‑149 (MQ=255)                                                                                                                                       cGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGAt         >  4:168918/1‑149 (MQ=255)                                                                                                                                         aGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCa       <  6:216231/149‑1 (MQ=255)                                                                                                                                              t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcagc                <  4:192123/135‑1 (MQ=255)                                                                                                                                              t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcagc                >  3:192123/1‑135 (MQ=255)                                                                                                                                              t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACg  >  4:121582/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     CTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCC‑GCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACG  >  NC_002947/5027144‑5027432

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