Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	897,504	T→C	intergenic (‑63/‑265)	SAUSA300_RS04405 ← / → SAUSA300_RS04410	DUF368 domain‑containing protein/ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	897,504	0	T	C	98.5%	180.0 / ‑5.0	68	intergenic (‑63/‑265)	SAUSA300_RS04405/SAUSA300_RS04410	DUF368 domain‑containing protein/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (54/13);  total (55/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

CTTTTAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGTCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCGGATTAGAG  >  NC_007793/897408‑897600                                                                                                 |                                                                                                  cTTTTAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCt                                                                                                <  1:3321050/100‑1 (MQ=255) cTTTTAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCt                                                                                                <  1:2909741/100‑1 (MQ=255) cTTTTAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAAGGCCCt                                                                                                <  1:4454025/100‑1 (MQ=255)  ttttAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAGAATCATAAATAGCATCATACCATATTACAAATGCCCTa                                                                                               <  1:3668373/100‑1 (MQ=255)  ttttAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCtg                                                                                               <  1:5187121/100‑2 (MQ=255)  ttatAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTa                                                                                               <  1:5716173/97‑1 (MQ=255)   tttAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAg                                                                                              <  1:264334/100‑1 (MQ=255)      aGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCtcacag                                                                                           <  1:4253555/100‑6 (MQ=255)              aaTCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATaa                                                                                   >  1:2260955/1‑100 (MQ=255)                          ttgttgCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAAt                                                                       >  1:3284297/1‑100 (MQ=255)                           agttgCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAAtt                                                                      >  1:5210877/2‑100 (MQ=255)                             ttgCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCa                                                                    >  1:2635126/1‑100 (MQ=255)                              tgCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCAt                                                                   >  1:5201487/1‑100 (MQ=255)                              tgCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCAt                                                                   >  1:2562925/1‑100 (MQ=255)                               gCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATa                                                                  >  1:1240798/1‑100 (MQ=255)                               gCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATa                                                                  <  1:628777/100‑1 (MQ=255)                               gCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATa                                                                  >  1:740706/1‑100 (MQ=255)                               gCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATa                                                                  >  1:2485631/1‑100 (MQ=255)                                cATGAATAACTCCTTCCAAACGATTGGAATAAAATCATAAATGGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATag                                                                 <  1:1768753/100‑2 (MQ=255)                                   gcttAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAAt                                                              >  1:3369992/4‑100 (MQ=255)                                   gagtAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAAt                                                              >  1:1039480/4‑100 (MQ=255)                                     aTAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATcc                                                            >  1:431819/1‑100 (MQ=255)                                       accTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCAt                                                          >  1:5051046/3‑100 (MQ=255)                                       aaCTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCAt                                                          >  1:4147147/1‑100 (MQ=255)                                       aaCTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCAt                                                          >  1:4522444/1‑100 (MQ=255)                                        aCTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCAtt                                                         >  1:4762841/1‑100 (MQ=255)                                          tCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTga                                                       >  1:680656/1‑100 (MQ=255)                                           ccTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTgag                                                      >  1:895213/1‑100 (MQ=255)                                           ccTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTgag                                                      >  1:2385391/1‑100 (MQ=255)                                           acTTTCAAACGATTGGAACAAAATCATAAATAGCATCATACCATATTACAAATGCCCTGGTGAAATGATAACATATTTTAAATTCATAAAATCCATTgag                                                      >  1:1252113/2‑100 (MQ=255)                                              ttCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAGATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGaaa                                                   >  1:2232373/1‑100 (MQ=255)                                              ttCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAGATGCCCCAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGaaa                                                   >  1:4814737/1‑100 (MQ=255)                                              ttCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGaaa                                                   >  1:5530163/1‑100 (MQ=255)                                              ttCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGaaa                                                   >  1:617305/1‑100 (MQ=255)                                              ttCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGaaa                                                   >  1:3574044/1‑100 (MQ=255)                                              ttCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGaaa                                                   >  1:3151993/1‑100 (MQ=255)                                                         ggAATAAAATCATAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACt                                        >  1:3085870/1‑100 (MQ=255)                                                                 actctAAATAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCa                                >  1:692839/5‑100 (MQ=255)                                                                        aTAGCACCATACCATATTACAAATGTCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTatat                         >  1:1020162/1‑100 (MQ=255)                                                                         tAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATAtt                        >  1:476787/1‑100 (MQ=255)                                                                         tAGCATCATACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATAtt                        >  1:4436495/1‑100 (MQ=255)                                                                            cctcatACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATAttttt                     >  1:3056995/3‑100 (MQ=255)                                                                                atACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAa‑g                <  1:5192117/100‑1 (MQ=255)                                                                                atACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTCAAa‑g                <  1:779896/100‑1 (MQ=255)                                                                                 tACCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTCAAAa‑ga               >  1:1044321/1‑100 (MQ=255)                                                                                  cccATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAa‑gag              >  1:675963/2‑100 (MQ=255)                                                                                  aCCATATTACAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAa‑gag              >  1:3119157/1‑100 (MQ=255)                                                                                          aCAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAG‑TAGATCGGaag      >  1:319126/1‑98 (MQ=255)                                                                                          aCAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCGGAtt      <  1:2812515/100‑1 (MQ=255)                                                                                          aCAAATGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTAAAA‑GAGATCGGaag      >  1:1579192/1‑98 (MQ=255)                                                                                           aaaaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCAGAtcg     >  1:2839050/2‑98 (MQ=255)                                                                                            aaaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAg‑gagagatcggaag    >  1:436257/1‑92 (MQ=255)                                                                                            aaaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGATCGGaagag    >  1:1303883/1‑96 (MQ=255)                                                                                            aaaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCAGAtcgg    >  1:5220716/1‑97 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAa‑gagagatcggaaga   >  1:1218790/1‑91 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGATCGGaagagc   >  1:1995413/1‑95 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCAGAtcgga   >  1:3515959/1‑96 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTAAAA‑GAGATCGGaagagc   >  1:4541010/1‑95 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTCTAAA‑GAGATCGGaagagc   >  1:1548453/1‑95 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTCTAAA‑GAGATCGGaagagc   >  1:3244979/1‑95 (MQ=255)                                                                                             aaTGCCCTAGTGAAATGATAACATATTTTAAATGCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCGGATTaga   >  1:1189698/1‑100 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTATA‑GAGATCGGaagagca  >  1:5632423/1‑94 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTACA‑GAGATCGGaagagca  >  1:505222/1‑94 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAATGAGATCGGaagagc   >  1:5311562/1‑95 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGATCGGaagagca  >  1:4026078/1‑94 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCGGATTagag  >  1:5289386/1‑100 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCGGATTagag  >  1:2302848/1‑100 (MQ=255)                                                                                              aTGCCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTCTAAA‑GAGATCGGaagagca  >  1:2451097/1‑94 (MQ=255)                                                                                                 |                                                                                                  CTTTTAGAATATTAATCCATTTAAATTGTTGCATGAATAACTCCTTTCAAACGATTGGAATAAAATCATAAATAGCATCATACCATATTACAAATGTCCTAGTGAAATGATAACATATTTTAAATTCATAAAATCCATTGAGAAATTATGTGCACTTATTATCATTTATATTTTTAAA‑GAGAGCGGATTAGAG  >  NC_007793/897408‑897600

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