Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,606,071	T→C	Q202R (CAG→CGG)	SAUSA300_RS07935 ←	ribonuclease Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,606,071	0	T	C	93.9%	74.3 / ‑3.2	33	Q202R (CAG→CGG)	SAUSA300_RS07935	ribonuclease Z
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (4/27);  total (4/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 = 2.33e-02

TAGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCA  >  NC_007793/1605975‑1606169                                                                                                 |                                                                                                   taGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGtt                                                                                                 <  1:1103621/100‑1 (MQ=255) taGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGtt                                                                                                 <  1:5299996/100‑1 (MQ=255) taGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGtt                                                                                                 <  1:3108766/100‑1 (MQ=255) taGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGta                                                                                                 <  1:2194809/100‑2 (MQ=255)  aGGTCGCTTCATGTACCATCCCATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGttg                                                                                                <  1:2671825/100‑2 (MQ=255)  aGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGttt                                                                                                <  1:3985342/100‑1 (MQ=255)  aGGTCGCCTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGttt                                                                                                <  1:5479023/100‑1 (MQ=255)   ggTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTg                                                                                               <  1:4873131/100‑1 (MQ=255)    gTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTgg                                                                                              <  1:4287569/100‑1 (MQ=255)    gTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTgg                                                                                              <  1:3532817/100‑1 (MQ=255)    gTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTgg                                                                                              <  1:2909451/100‑1 (MQ=255)    gTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTgg                                                                                              <  1:1900180/100‑1 (MQ=255)       gCTTCATGTGCCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACt                                                                                           <  1:1180280/100‑1 (MQ=255)       gCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACt                                                                                           <  1:2275748/100‑1 (MQ=255)       gCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACt                                                                                           <  1:3481293/100‑1 (MQ=255)       gCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACt                                                                                           <  1:4647286/100‑1 (MQ=255)        cTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACtg                                                                                          <  1:630691/100‑2 (MQ=255)          tCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCAc                                                                                        <  1:4675460/100‑1 (MQ=255)             tGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACctc                                                                                     <  1:2048396/100‑1 (MQ=255)              gTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACctct                                                                                    <  1:3693965/100‑1 (MQ=255)              gTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACctct                                                                                    <  1:4376747/100‑1 (MQ=255)               tACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACGTGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTa                                                                                   <  1:4046612/100‑1 (MQ=255)               tACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTa                                                                                   <  1:882017/100‑1 (MQ=255)                    acacatcaGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTAAAATc                                                                              >  1:4387889/2‑100 (MQ=255)                       catcaGCATCTCTAGAAATCACACGTTCATTAGAACGTGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTAAAATCttt                                                                           >  1:4120231/1‑100 (MQ=255)                        atcaGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTAAAATCtttt                                                                          <  1:47318/100‑1 (MQ=255)                          caGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTAAAATCTTttga                                                                        >  1:5310372/1‑100 (MQ=255)                                                  cATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACCCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGtaa                                                <  1:5238632/100‑3 (MQ=255)                                                           ttGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCCGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATc                                       >  1:4283610/2‑100 (MQ=255)                                                                                           ctcccccGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCtt       <  1:408341/97‑2 (MQ=255)                                                                                            gACCCCGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTgg      <  1:4033777/100‑1 (MQ=255)                                                                                            gACCCCGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTgg      <  1:3304179/100‑1 (MQ=255)                                                                                              ccccGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGtt    <  1:5790003/100‑1 (MQ=255)                                                                                               tcTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTc   <  1:2039202/99‑1 (MQ=255)                                                                                                cTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCa  <  1:3370442/100‑1 (MQ=255)                                                                                                 |                                                                                                   TAGGTCGCTTCATGTACCATCACATCAGCATCTCTAGAAATCACACGTTCATTAGAACATGGTTTTGTATCACCAAAAATTGCTACAACTGGACCCTGTTTGGACTCACCTCTAAAATCTTTTGATTGATAAACTTGACCATTATGTTCAAATGTATCATGAGATTTTACTTCTTGATATTTAGGACCTGGTTCA  >  NC_007793/1605975‑1606169

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