Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	234,575	(A)6→5	coding (881/1113 nt)	SAUSA300_RS01050 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	234,570	0	A	.	98.4%	240.9 / ‑3.1	64	coding (876/1113 nt)	SAUSA300_RS01050	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (62/1);  total (63/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

TTGCGAATTTAACTAATTTAACTAGTAATGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATTAAAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATATTTAGCTGATTATATTCGAAGCCTTA  >  NC_007793/234485‑234664                                                                                      |                                                                                               ctGCGAATTTAACTAATTTAACTAGTAATGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATaaa                                                                                 >  1:1620737/2‑100 (MQ=255)                          tAATGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAAc                                                        >  1:2442522/1‑100 (MQ=255)                           caTGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACa                                                       >  1:4537092/2‑100 (MQ=255)                           aaTGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACa                                                       >  1:1108267/1‑100 (MQ=255)                           aaTGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACa                                                       >  1:20452/1‑100 (MQ=255)                           aaTGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACa                                                       >  1:5053090/1‑100 (MQ=255)                              gCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAgg                                                    >  1:3452652/1‑100 (MQ=255)                              gCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAgg                                                    >  1:1952955/1‑100 (MQ=255)                              gCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCAATGagatcga                                                    >  1:5205205/1‑94 (MQ=255)                                  ccAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATa                                                >  1:1799310/2‑100 (MQ=255)                                   cAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAg                                               >  1:2233616/1‑100 (MQ=255)                                     aaGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGat                                             >  1:1259094/1‑100 (MQ=255)                                     aaGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGat                                             >  1:2101093/1‑100 (MQ=255)                                     aaGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGat                                             >  1:1785849/1‑100 (MQ=255)                                     aaGCTGGCGATGATGGATTTAGAAAGTATTTTATAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGat                                             >  1:4508028/1‑100 (MQ=255)                                      ttatGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatt                                            >  1:1746784/4‑100 (MQ=255)                                      tgCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatt                                            >  1:1735205/2‑100 (MQ=255)                                      tgCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑‑AAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:1981246/2‑100 (MQ=255)                                      tccTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatt                                            >  1:4539514/3‑100 (MQ=255)                                      ggCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatt                                            >  1:2156342/2‑100 (MQ=255)                                      ctttGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatt                                            >  1:634462/4‑100 (MQ=255)                                      ctcTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatt                                            >  1:1337430/3‑100 (MQ=255)                                       gCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:3823073/1‑100 (MQ=255)                                       gCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:5720019/1‑100 (MQ=255)                                       gCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:1399077/1‑100 (MQ=255)                                       gCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:2942146/1‑100 (MQ=255)                                       gCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:266248/1‑100 (MQ=255)                                       ccTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:1658481/2‑100 (MQ=255)                                       ccTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatta                                           >  1:3803179/2‑100 (MQ=255)                                        tcggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGattat                                          >  1:966978/3‑100 (MQ=255)                                        gtGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGattat                                          >  1:1783417/2‑100 (MQ=255)                                        cTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGattat                                          >  1:3366102/1‑100 (MQ=255)                                        cTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGattat                                          >  1:56119/1‑100 (MQ=255)                                         tGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGattatt                                         >  1:2115671/1‑100 (MQ=255)                                         tGGCGATGACGGATTTAGAAAGAATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACATGGAATAGatcgga                                         >  1:4995846/1‑96 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatcggaa                                        >  1:1748462/1‑95 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGatcggaa                                        >  1:3782104/1‑95 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:5263512/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:5097451/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:4984091/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:4415568/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:4074144/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:3960420/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:1921172/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:2047411/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:3598190/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:1603158/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:1487869/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:2861996/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:2685287/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:2533808/1‑100 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAACAGatcggaa                                        >  1:1897334/1‑95 (MQ=255)                                          ggCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACA‑TAAAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTAttt                                        >  1:1641369/1‑100 (MQ=255)                                           gCGATGATGGATTTAGAGAGTATTTTGTAAATTACCAAATATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTg                                       >  1:5449650/1‑100 (MQ=255)                                           gCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTg                                       >  1:878555/1‑100 (MQ=255)                                           gCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTg                                       >  1:1701792/1‑100 (MQ=255)                                            cGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAAGAAGGAATAGATTATTTGc                                      >  1:1400198/1‑100 (MQ=255)                                              atgatgGATTTAGAAAGTATTTTGTAAATTACCACACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCa                                    >  1:5237707/1‑100 (MQ=255)                                                gatgGATTTAGAGAGCATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCAtt                                  >  1:2475254/1‑100 (MQ=255)                                                         aGAAAGTATTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATAttt                         >  1:2412241/1‑100 (MQ=255)                                                                aTTTTGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATATTTAGCTGag                  <  1:898348/100‑2 (MQ=255)                                                                    tGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATATTTAGCTGATtata              >  1:3836933/1‑100 (MQ=255)                                                                    tGTAAATTACCAAACATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATATTTAGCTGATtata              >  1:3064810/1‑100 (MQ=255)                                                                                aaCATT‑AAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATATTTAGCTGATTATATTCGAAGCCTTa  >  1:2700894/1‑100 (MQ=255)                                                                                      |                                                                                               TTGCGAATTTAACTAATTTAACTAGTAATGCAATCAAAGCTGGCGATGATGGATTTAGAAAGTATTTTGTAAATTACCAAACATTAAAAAAGATTCATAAAGGAATAACAAATCCAAGCGATGAACAAGGAATAGATTATTTGCCATTTCAATATTTAGCTGATTATATTCGAAGCCTTA  >  NC_007793/234485‑234664

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