Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,624,971	(T)7→8	coding (1280/1335 nt)	SAUSA300_RS13505 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,624,964	1	.	T	92.5%	104.6 / 2.0	40	coding (1287/1335 nt)	SAUSA300_RS13505	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base T (0/37);  total (0/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

TTTTTGTCCTTCTCTTGTTACTTTTTCTGTCCCTGGTGCTAAATCCGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAAC‑‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCATTTGCATTT  >  NC_007793/2624878‑2625057                                                                                        |                                                                                                 tttttGTCCTTCTCTTGTTACTTTTTCTGTCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAAC‑‑‑ttttttttctaag                                                                                   <  1:2418159/100‑7 (MQ=255)      tccgatctcTTGCTACTTTTTCTGTCCCTGGTACTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGt                                                                              <  1:4267691/95‑1 (MQ=17)               ttGTTACTTTTTCTGTCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGtttt                                                                     <  1:605307/100‑1 (MQ=255)                   tattttttCTGTCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGgc                                                                 <  1:5574732/97‑1 (MQ=18)                   tACTTTTTCTGTCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAAC‑‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGgcg                                                                <  1:1865089/100‑1 (MQ=255)                    tcTTTTTCTGTCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGgcg                                                                <  1:2795328/99‑1 (MQ=255)                             gTCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAAC‑‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTTTGgtgt                                                      <  1:4719307/100‑1 (MQ=255)                              tCCCTGGTGCTAGATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGgtgt                                                      <  1:4604127/100‑1 (MQ=25)                              tCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGGATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGgtgt                                                      <  1:2593309/100‑1 (MQ=25)                              tCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGgtgg                                                      <  1:2030566/100‑2 (MQ=255)                              tCCCTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCACTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGgtgt                                                      <  1:894736/100‑1 (MQ=25)                              tCCCTGGTGCTAAATCAGGATTAAATTTACGTTCCTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGgtgt                                                      <  1:90405/100‑1 (MQ=25)                               cccTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTa                                                     <  1:3606107/100‑1 (MQ=255)                                ccTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAGCT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTaa                                                    <  1:681245/100‑1 (MQ=25)                                 cTGGTGCTAAATCAGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTaaa                                                   <  1:5295579/100‑1 (MQ=255)                                            tcaggATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTg                                        <  1:2085782/97‑1 (MQ=255)                                                gATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTgcata                                    <  1:3799128/100‑5 (MQ=255)                                                  ttAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCtt                                  <  1:4001722/100‑1 (MQ=255)                                                   tAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCttt                                 <  1:4302392/100‑1 (MQ=255)                                                     ctcttccgatctttcctAAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCa                               <  1:1508320/91‑1 (MQ=255)                                                     aaTTTACGTTCTTTTTTGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCa                               <  1:4936969/100‑1 (MQ=255)                                                         tACGTTCTTTCTTGAACGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTc                           <  1:3963040/100‑1 (MQ=255)                                                         tACGTTCTTTCTTAAATGAATTCTCTTCAGCT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTc                           <  1:1549410/100‑1 (MQ=255)                                                           cGTTCTTTCTTGAATGGAATCTCTTCAACTTTTTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTc                           <  1:6250271/100‑1 (MQ=255)                                                                tttcttGAATGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTtctc                    <  1:4450506/100‑1 (MQ=255)                                                                       cttccgaTCTCTTCATCT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAgt             <  1:3974080/94‑2 (MQ=255)                                                                        aTGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGcc            <  1:6536416/100‑1 (MQ=255)                                                                        aTGGAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGcc            <  1:158258/100‑1 (MQ=255)                                                                          ggAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAt          <  1:514316/100‑1 (MQ=255)                                                                           gAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCaga         <  1:2943106/100‑3 (MQ=255)                                                                           gAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAtt         <  1:2949800/100‑1 (MQ=255)                                                                           gAATCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAta         <  1:3501484/100‑2 (MQ=255)                                                                            aaTCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAttt        <  1:994582/100‑1 (MQ=255)                                                                            aaTCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAttg        <  1:958248/100‑2 (MQ=255)                                                                            aaTCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAtta        <  1:1500413/100‑2 (MQ=255)                                                                            aaTCTCTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCAtag        <  1:6346334/100‑3 (MQ=255)                                                                              tctcTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGACCCTTTCAAGTCAGTTCTCATCCAGCCATTTGc      <  1:2036444/100‑1 (MQ=255)                                                                              tctcTTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCATTTGc      <  1:1438181/100‑1 (MQ=255)                                                                                 ctTCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGATCCTTTCAAGTCAGTTCTCATCCAGCCATTTGCAtt   <  1:2676383/100‑1 (MQ=255)                                                                                  ttCAACT‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCATTTGCAttt  <  1:4141602/100‑1 (MQ=255)                                                                                  ttCAAC‑‑‑TTTTTTTCCTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCATTTGCAtta  <  1:2476925/100‑2 (MQ=255)                                                                                        |                                                                                                 TTTTTGTCCTTCTCTTGTTACTTTTTCTGTCCCTGGTGCTAAATCCGGATTAAATTTACGTTCTTTCTTGAATGGAATCTCTTCAAC‑‑‑TTTTTTT‑CTAATTCTGTTATTGTTTTTGGCGCTTCTGGTGTAAAAGTAAACTCTGAACCTTTCAAGTCAGTTCTCATCCAGCCATTTGCATTT  >  NC_007793/2624878‑2625057

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