Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,690,352	A→G	L713L (CTT→CTC)	SAUSA300_RS13825 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,690,352	0	A	G	98.0%	123.7 / ‑4.8	49	L713L (CTT→CTC)	SAUSA300_RS13825	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/31);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

TTGCACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCAAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGTTAATGCAGG  >  NC_007793/2690258‑2690447                                                                                               |                                                                                                ttgCACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGc                                                                                            <  1:4317061/100‑1 (MQ=255) ttgCACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGc                                                                                            <  1:4855078/100‑1 (MQ=255)   gCACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCAt                                                                                          <  1:183577/100‑1 (MQ=255)    cACCTTGTGGATATTCTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCAtt                                                                                         <  1:4970789/100‑2 (MQ=255)    cACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCAta                                                                                         <  1:4749282/100‑2 (MQ=255)    cACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATg                                                                                         <  1:146340/100‑1 (MQ=255)       cTTGTGTGTATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTcc                                                                                      <  1:3992877/100‑1 (MQ=255)         tgtgGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCtt                                                                                    >  1:2027897/1‑100 (MQ=255)          gtgGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTc                                                                                   >  1:1833875/1‑100 (MQ=255)           tgGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCa                                                                                  <  1:3228271/100‑1 (MQ=255)                        aaGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTg                                                                     <  1:5276409/100‑1 (MQ=255)                         aGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCTTTCAATATCAACTCGTgg                                                                    <  1:788070/100‑1 (MQ=255)                             tgtgTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGaata                                                                <  1:4427681/100‑1 (MQ=255)                             tgtgTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGaata                                                                >  1:3395859/1‑100 (MQ=255)                                gtCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGaatacgg                                                             <  1:3899113/100‑4 (MQ=255)                                 tCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGaatagggc                                                            <  1:5243531/100‑5 (MQ=255)                                     aTGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGAtt                                                        >  1:199596/1‑100 (MQ=255)                                          ctcttccgaTCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAAcc                                                   <  1:4839952/92‑1 (MQ=255)                                             ttccgaTCTCTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATg                                                <  1:59153/95‑1 (MQ=255)                                              tACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGc                                               >  1:1586557/1‑100 (MQ=255)                                               ccgaTCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCt                                              <  1:1013949/97‑1 (MQ=255)                                               aCTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCt                                              >  1:4646710/1‑100 (MQ=255)                                                tactCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTg                                             <  1:2183243/97‑1 (MQ=255)                                                ctttCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTg                                             <  1:5648417/97‑1 (MQ=255)                                                 tATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTgg                                            <  1:2432643/100‑1 (MQ=255)                                                 tATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTgg                                            <  1:217506/100‑1 (MQ=255)                                                 tATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTgg                                            <  1:3735832/100‑1 (MQ=255)                                                  aTCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTggc                                           <  1:4104062/100‑2 (MQ=255)                                                  aTCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGt                                           <  1:38861/100‑1 (MQ=255)                                                   tCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTggag                                          <  1:3540550/100‑3 (MQ=255)                                                      ttAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAAt                                       >  1:3700786/1‑100 (MQ=255)                                                      ttAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAAt                                       >  1:503653/1‑100 (MQ=255)                                                        aGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATAc                                     >  1:2377100/1‑100 (MQ=255)                                                                    cTGTGATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCtttt                         >  1:4338599/1‑100 (MQ=255)                                                                     tgtgATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTc                        >  1:3090455/1‑100 (MQ=255)                                                                        gATTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCaa                     <  1:3157662/100‑1 (MQ=255)                                                                         aTTCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCcaa                    <  1:5138118/100‑1 (MQ=255)                                                                          ttCAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCacac                   <  1:180633/100‑1 (MQ=255)                                                                           tCAACCGTTTTATACTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACa                  >  1:2876714/1‑100 (MQ=255)                                                                            taaCCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACaa                 >  1:1540624/2‑100 (MQ=255)                                                                            cggccGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACaa                 >  1:546057/4‑100 (MQ=255)                                                                            cAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACaa                 >  1:3057817/1‑100 (MQ=255)                                                                            cAACCGTTTTATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACaa                 >  1:2259164/1‑100 (MQ=255)                                                                             aaCCGTTTTATACTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAAt                >  1:829965/1‑100 (MQ=255)                                                                                    ttATATTTTTCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGttt         <  1:2701890/100‑2 (MQ=255)                                                                                        atTTTTCAAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGTTAATGc     <  1:4066294/100‑1 (MQ=255)                                                                                        aagtttCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGTTAATGc     <  1:4598989/97‑1 (MQ=255)                                                                                           tttCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGTTAATGCAgg  <  1:4888262/100‑1 (MQ=255)                                                                                           tttCGAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGTTAATGCAgg  <  1:5046911/100‑1 (MQ=255)                                                                                               |                                                                                                TTGCACCTTGTGGATATACTTTAAAGGTTGTGTCGTATGTTTCCTTACTATCTTTAGCTTCAGATTCCTGTGATTCAACCGTTTTATATTTTTCAAGTGCATGTCCTTCAATATCAACTCGTGGAATAATGCGATTCAACCATGCTGGTAAATACCACGAACCTTTTCCAAACAATTTCGTTAATGCAGG  >  NC_007793/2690258‑2690447

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