Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,519,175	Δ1 bp	coding (386/588 nt)	SAUSA300_RS07370 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,519,174	0	C	.	97.6%	155.7 / ‑2.5	42	coding (387/588 nt)	SAUSA300_RS07370	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (7/34);  total (7/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

AATCATAATATGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAACCTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGTTCAAT  >  NC_007793/1519077‑1519263                                                                                                  |                                                                                          aaTCATAATATGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAActa                                                                                         <  1:5208114/100‑3 (MQ=255)   tCATAATATGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAActagg                                                                                       <  1:5619078/100‑5 (MQ=255)     ataataTGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCTAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGAt                                                                                    <  1:1439734/100‑1 (MQ=255)     ataataTGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGAt                                                                                    <  1:5744985/100‑1 (MQ=255)        atatGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATaac                                                                                 <  1:2380271/100‑3 (MQ=255)                    aaGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCt                                                                     <  1:1097837/100‑1 (MQ=255)                     aGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAACCTAGGATAGAAGTAAATTGCCt                                                                     <  1:5588301/100‑1 (MQ=255)                     aGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTg                                                                    <  1:5143244/100‑1 (MQ=255)                     aGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTg                                                                    <  1:2333807/100‑1 (MQ=255)                      ggTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGa                                                                   <  1:2122181/100‑1 (MQ=255)                            aaTAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACGGcc                                                             <  1:3049887/100‑1 (MQ=255)                            aaTAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGcc                                                             <  1:4222955/100‑1 (MQ=255)                             aTAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCa                                                            <  1:4486839/100‑1 (MQ=255)                             aTAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCa                                                            <  1:2838559/100‑1 (MQ=255)                              tAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCaa                                                           <  1:5731994/100‑2 (MQ=255)                               aGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTCCCTGAACAGCCaga                                                          <  1:3895245/100‑3 (MQ=255)                                gTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCAttt                                                         <  1:5181156/100‑2 (MQ=255)                                 tCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATtat                                                        <  1:1701809/100‑1 (MQ=255)                                      ataATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACtgt                                                   <  1:1802950/100‑1 (MQ=255)                                        aaTCATTCAAAAACGCCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTgg                                                 <  1:3074904/100‑1 (MQ=255)                                        aaTCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTgg                                                 <  1:2251924/100‑1 (MQ=255)                                         aTCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTggg                                                <  1:3270140/100‑1 (MQ=255)                                          tCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTgggc                                               <  1:5117876/100‑2 (MQ=255)                                                 gaaaCACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACAtt                                        >  1:4434272/2‑100 (MQ=255)                                                   tgcacCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTaa                                      >  1:1401129/3‑100 (MQ=255)                                                    acacCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAAt                                     >  1:5181504/1‑100 (MQ=255)                                                         ctcttccgatcTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCaa                                <  1:3416488/90‑1 (MQ=255)                                                         cAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCaa                                <  1:4656486/100‑1 (MQ=255)                                                              catcGCTACACTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCAt                           >  1:4149538/1‑100 (MQ=255)                                                               ctcggttcacTGATTCCAATCATAGAACGTTTAAA‑CTAGGGTAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCAtt                          <  1:4126247/93‑1 (MQ=255)                                                                     acacTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGt                    <  1:2099019/100‑1 (MQ=255)                                                                     acacTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGt                    >  1:4725990/1‑100 (MQ=255)                                                                     acacTGATTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGt                    <  1:2280752/100‑1 (MQ=255)                                                                           aTTCCAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATc              <  1:2989414/100‑1 (MQ=255)                                                                            cttccgatctAGAACTTTTAAA‑CTAGGATAGAATTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCa             <  1:1849424/91‑1 (MQ=255)                                                                             tCCATTCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCAc            <  1:4744781/100‑1 (MQ=255)                                                                             gccAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCAc            >  1:5666612/2‑100 (MQ=255)                                                                               cAATCATAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCt          <  1:1918665/100‑1 (MQ=255)                                                                                   cttAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGtt      <  1:5069028/98‑1 (MQ=255)                                                                                    aTAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTgctt     <  1:5225061/100‑4 (MQ=255)                                                                                    aTAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGTTc     <  1:2269880/100‑1 (MQ=255)                                                                                    aTAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGTTc     <  1:2898599/100‑1 (MQ=255)                                                                                     tAGAACGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGTTCa    <  1:5087094/100‑1 (MQ=255)                                                                                       acaCGTTTAAA‑CTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGTTCaat  >  1:1581047/3‑100 (MQ=255)                                                                                                  |                                                                                          AATCATAATATGGAAATTGAAGGTGAAAATAGTCAATATAATCATTCAAAAACACCCAAATCATCGCTACACTGATTCCAATCATAGAACGTTTAAACCTAGGATAGAAGTAAATTGCCTGAACAGCCATTATACTGTGGGAAAACATTAATACCAAACCATTTACTGTAATATCACCTTGTTCAAT  >  NC_007793/1519077‑1519263

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