Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,680,195	C→T	11.1%	L96L (TTG→TTA)	cdd ←	cytidine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,680,195	0	C	T	11.1%	93.2 / 3.4	36	L96L (TTG→TTA)	cdd	cytidine deaminase
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base T (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GATTTATGTTCTGTCATTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGTTCAGCACATATCGATAAT  >  CP000730/1680061‑1680333                                                                                                                                       |                                                                                                                                           gATTTATGTTCTGTCATTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACt                                                                                                                                      >  1:209200/1‑141 (MQ=255)           ctttCATTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTTAAAACTTGACGACATg                                                                                                                            <  2:393599/138‑1 (MQ=255)            tGTCATTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGc                                                                                                                           <  1:256083/141‑1 (MQ=255)              tCATTTATTCTAAATCCTTTCCTGAAAATCCAAATTGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGcac                                                                                                                         <  1:203505/141‑1 (MQ=255)               cATTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGcacc                                                                                                                        <  2:176950/141‑1 (MQ=255)                aTTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTTAAAACTTGACGACATGcacca                                                                                                                       >  2:179560/1‑141 (MQ=255)                         aaaTCCTTTCCTGAAAATCCAAATGGTAGTAACTCGGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGtgat                                                                                                              >  1:102866/1‑141 (MQ=255)                                 tCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGttt                                                                                                      >  1:240365/1‑141 (MQ=255)                                 tCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGttt                                                                                                      >  2:92859/1‑141 (MQ=255)                                       aaaTCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGc                                                                                                <  1:162117/141‑1 (MQ=255)                                            cAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACTAGGCAATGACTGTTTATCTGCATCTa                                                                                           >  2:433125/1‑141 (MQ=255)                                              aaTGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACg                                                                                         >  2:83558/1‑141 (MQ=255)                                                   tagtaACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTAc                                                                                    <  2:452277/141‑1 (MQ=255)                                                      taACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGt                                                                                 >  1:260560/1‑141 (MQ=255)                                                                aCTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCa                                                                       <  2:274388/141‑1 (MQ=255)                                                                   gTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCaaaa                                                                    >  1:198797/1‑141 (MQ=255)                                                                       tcatAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCAc                                                                >  2:286965/1‑141 (MQ=255)                                                                           agtcATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCTGATAATTTCTTCGAATCTTGACGACGTGTACCACAAGGTGATGGCTGTTTATCTGCATCTGCTGATACAGTTATTGGTTCTAAATCACCTgg                                                            <  2:29952/138‑1 (MQ=255)                                                                           aaCCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTgg                                                            <  1:7422/141‑1 (MQ=255)                                                                               gcatCTACGTTATGATTTGTCATATACACCGGCATATCAACAGCACTTAATTGCTTCACAACTTGACGACGTGCACCACGAGGTGATGACGGTTAATCTGCATCTACGGTTACAGTTATTGCTCCAAAATCACCTGGTCtc                                                        <  2:249589/139‑2 (MQ=255)                                                                                tatCTCCTTTATGCTTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATCCACCACAACGTGATGACGGTTTAGCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTTGTCTGt                                                       <  2:116827/141‑1 (MQ=255)                                                                                    tcCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTTAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATcc                                                   <  1:280611/141‑1 (MQ=255)                                                                                       tttATGATTTGTCATATACACAGGCATATCATCATCACATAATCCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTg                                                >  1:341090/1‑141 (MQ=255)                                                                                        ttATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTga                                               >  1:46107/1‑141 (MQ=255)                                                                                              tttGTCATATACACAGGCATATCATCATCACATAATTCCTTTAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAAtt                                         >  1:43362/1‑141 (MQ=255)                                                                                               ttGTCATATACACAGGCATATCATCATCATATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTTAAAATCACCTGGTCTGTATCCTTGATAAATTg                                        >  2:326082/1‑141 (MQ=255)                                                                                                tGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGc                                       >  2:92002/1‑141 (MQ=255)                                                                                                    atatACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGAt                                   <  2:235560/141‑1 (MQ=255)                                                                                                          acaGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTAGCCTTGAGAAATTGCCGATACCaaa                             >  2:448868/1‑141 (MQ=255)                                                                                                              gCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAg                         >  1:205827/1‑141 (MQ=255)                                                                                                                   tcatcatcaCATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGt                    >  2:179073/1‑141 (MQ=255)                                                                                                                   tcatcatcaCATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGt                    >  1:384627/1‑141 (MQ=255)                                                                                                                          cacaTAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGTTCAGcac             >  1:342201/1‑141 (MQ=255)                                                                                                                            caTAATGCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGATAAATTGCCGATACCAAACTAGCTCGTTCAGCACat           <  2:38952/141‑1 (MQ=255)                                                                                                                             aTAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGTTCAGCACata          <  2:379675/141‑1 (MQ=255)                                                                                                                                     ttCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGTTCAGCACATATCGATAAt  <  1:175292/141‑1 (MQ=255)                                                                                                                                       |                                                                                                                                           GATTTATGTTCTGTCATTTATTCTAAATCCTTTCCTGAAAATCCAAATGGTAGTAACTCTGCGACTGTCATCATAACCATATCTCCTTTATGATTTGTCATATACACAGGCATATCATCATCACATAATTCCTTCAAAACTTGACGACATGCACCACAAGGTGATGACGGTTTATCTGCATCTACGGTTACAGTTATTGATTCAAAATCACCTGGTCTGTATCCTTGAGAAATTGCCGATACCAAACTAGCTCGTTCAGCACATATCGATAAT  >  CP000730/1680061‑1680333

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑