Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,198,908	A→G	Q40R (CAG→CGG)	nfi →	endonuclease V

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,198,908	0	A	G	97.9%	155.4 / ‑5.0	47	Q40R (CAG→CGG)	nfi	endonuclease V
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (18/28);  total (18/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.96e-01

TCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGA  >  NC_000913/4198771‑4199034                                                                                                                                          |                                                                                                                               tCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCggg                                                                                                                              >  2:282112/1‑140 (MQ=255)  cACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCgggg                                                                                                                             <  1:408499/140‑1 (MQ=255)    ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCg                                                                                                                           <  1:93886/140‑1 (MQ=255)    ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCg                                                                                                                           >  2:475593/1‑140 (MQ=255)    ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCg                                                                                                                           >  1:570596/1‑140 (MQ=255)                 aTTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGAcgcg                                                                                                              <  1:333473/140‑1 (MQ=255)                    aTGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAgc                                                                                                           >  1:31267/1‑140 (MQ=255)                      ggATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAgccg                                                                                                         >  2:759438/1‑138 (MQ=255)                            cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATgg                                                                                                   >  2:420645/1‑140 (MQ=255)                            cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATgg                                                                                                   >  2:289159/1‑140 (MQ=255)                              cgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtg                                                                                                 >  1:757369/1‑140 (MQ=255)                              cgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtg                                                                                                 <  2:31267/140‑1 (MQ=255)                              cgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtg                                                                                                 <  1:475593/140‑1 (MQ=255)                                tCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtgct                                                                                               <  2:716702/140‑1 (MQ=255)                                 cATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                                              <  1:433126/140‑1 (MQ=255)                                 cATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                                              <  2:464368/140‑1 (MQ=255)                                  aTTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctgc                                                                                             >  1:664006/1‑140 (MQ=255)                                      cgcgCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGt                                                                                                  <  2:793081/131‑1 (MQ=255)                                      cgcgCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGt                                                                                                  >  1:793081/1‑131 (MQ=255)                                            caacaaATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATcc                                                                                   <  1:282112/140‑1 (MQ=255)                                            caacaaATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATcc                                                                                   <  1:774818/140‑1 (MQ=255)                                            caacaaATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATcc                                                                                   <  2:570596/140‑1 (MQ=255)                                                  aTCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCt                                                                             <  1:289159/140‑1 (MQ=255)                                                  aTCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCt                                                                             <  1:420645/140‑1 (MQ=255)                                                      aaCTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAg                                                                         <  1:242199/140‑1 (MQ=255)                                                           gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGt                                                                    <  1:347611/140‑1 (MQ=255)                                                            cttcttCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTc                                                                   <  2:33505/140‑1 (MQ=255)                                                                ttctGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                               <  2:757369/140‑1 (MQ=255)                                                                ttctGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                               >  2:740228/1‑140 (MQ=255)                                                                     tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATaaa                                                          >  1:770183/1‑140 (MQ=255)                                                                          cgcgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGc                                                     >  1:680338/1‑140 (MQ=255)                                                                            cgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTTACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGccc                                                   >  2:568977/1‑140 (MQ=255)                                                                            cgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGccc                                                   <  2:522101/140‑1 (MQ=255)                                                                              aGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGc                                                 <  1:603186/140‑1 (MQ=255)                                                                                   cGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGc                                            <  2:236357/140‑1 (MQ=255)                                                                                    gACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGc                                                 >  1:579579/1‑134 (MQ=255)                                                                                    gACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGc                                                 <  2:579579/134‑1 (MQ=255)                                                                                     aCTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGcca                                          <  2:746341/140‑1 (MQ=255)                                                                                        cGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccacca                                       <  1:216745/140‑1 (MQ=255)                                                                                        cGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccacca                                       <  2:581602/140‑1 (MQ=255)                                                                                          aTAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccacc                                     <  1:740228/140‑1 (MQ=255)                                                                                          aTAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccacc                                     >  2:49277/1‑140 (MQ=255)                                                                                           tAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccacca                                    >  1:382091/1‑140 (MQ=255)                                                                                            aaaGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCAt                                   <  1:765295/140‑1 (MQ=255)                                                                                                           cTGATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAgg                    <  2:770183/140‑1 (MQ=255)                                                                                                             gATCGCCGGAGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGtt                  >  1:545468/1‑140 (MQ=255)                                                                                                                      aGCCGATGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTttcct         <  1:49277/140‑1 (MQ=255)                                                                                                                           aTGTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCcgc    >  2:487324/1‑140 (MQ=255)                                                                                                                             gTCGGGTTTGAGCGGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGa  <  2:664006/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                               TCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGA  >  NC_000913/4198771‑4199034

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