Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,144,826	(A)7→8	coding (624/627 nt)	SAUSA300_RS10920 →	nitroreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,144,819	1	.	A	96.6%	154.6 / ‑1.2	59	coding (617/627 nt)	SAUSA300_RS10920	nitroreductase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (55/2);  total (57/2)
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.04e-01

TTATTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAAT‑‑AAAAAAATAATCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGT  >  NC_007793/2144738‑2144891                                                                                   |                                                                          ttATTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCaaggc                                                         <  1:1204147/100‑4 (MQ=255)  tATTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAatata                                                        >  1:4506284/1‑100 (MQ=255)  tATTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGCAAGATCGTTTCGTAATA‑AAAAAAATAATCaacaga                                                        >  1:1392702/1‑96 (MQ=255)   aTTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACGTTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAatatat                                                       >  1:5806125/1‑100 (MQ=255)   aTTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAatatat                                                       >  1:4009822/1‑100 (MQ=255)    ttGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATAtt                                                      >  1:653518/1‑100 (MQ=255)     tGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTg                                                     >  1:1668320/1‑100 (MQ=255)     tGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTg                                                     >  1:5665932/1‑100 (MQ=255)     tGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTg                                                     >  1:4311302/1‑100 (MQ=255)         tgaaTGCCATTTGGTGATATTCGTGAAGCTGCTGTCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAg                                                 >  1:2751005/3‑100 (MQ=255)         cAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATAAAAAAAAATAATCAATATATTGTTa                                                  >  1:5962773/1‑100 (MQ=255)         cAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATAAAAAAAAATAATCAATATATTGTTa                                                  >  1:5635624/1‑100 (MQ=255)         cAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAg                                                 >  1:6217822/1‑100 (MQ=255)         cAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAg                                                 >  1:231013/1‑100 (MQ=255)         cAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAg                                                 >  1:828896/1‑100 (MQ=255)         cAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGATTCGTAATA‑AAAAAAATAATCAatataaagatcg                                                 >  1:267059/1‑93 (MQ=255)          aaaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATAAAAAAAAATAATCAATATATTGttgg                                                 >  1:5265708/1‑98 (MQ=255)          aaaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAgg                                                >  1:4658761/1‑100 (MQ=255)          aaaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAgg                                                >  1:2243728/1‑100 (MQ=255)          aaaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAgg                                                >  1:3351211/1‑100 (MQ=255)           taTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAgat                                               >  1:824260/2‑98 (MQ=255)           aaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTTGTAATA‑AAAAAAATAATCAatttagagatcgga                                               >  1:4560621/1‑91 (MQ=255)           aaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGa                                               >  1:2624664/1‑100 (MQ=255)           aaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGa                                               >  1:4359605/1‑100 (MQ=255)           aaTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGa                                               >  1:3040824/1‑100 (MQ=255)            ttGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaa                                              >  1:737933/2‑100 (MQ=255)            aTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaa                                              >  1:6446994/1‑100 (MQ=255)            aTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaa                                              >  1:3861821/1‑100 (MQ=255)            aTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaa                                              >  1:5695174/1‑100 (MQ=255)            aTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTCCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaa                                              >  1:1084935/1‑100 (MQ=255)            aTGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTAGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaa                                              >  1:4608034/1‑100 (MQ=255)             tGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaac                                             >  1:6200632/1‑99 (MQ=255)             tGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaac                                             >  1:4812191/1‑99 (MQ=255)             tGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAt                                             >  1:5021750/1‑100 (MQ=255)             tGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAt                                             >  1:3383117/1‑100 (MQ=255)             cgCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAt                                             >  1:2619658/2‑100 (MQ=255)              gCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGagat                                            >  1:137297/1‑97 (MQ=255)              gCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaagg                                            >  1:6161405/1‑98 (MQ=255)              gCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaacg                                            >  1:409816/1‑98 (MQ=255)              gCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaacg                                            >  1:6059406/1‑98 (MQ=255)              gCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGaacg                                            >  1:6602745/1‑98 (MQ=255)              gCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAAT‑‑AAAAAAATAATCAATATATTGTTAGGaagat                                           >  1:6296123/1‑97 (MQ=255)                cATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAgggg                                          >  1:698285/1‑100 (MQ=255)                 ctatGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGCTAGGaagatcg                                         >  1:3301519/4‑95 (MQ=255)                   ttGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAACAGGACa                                       >  1:3784580/1‑100 (MQ=255)                   ctgctGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGACa                                       >  1:4751419/5‑100 (MQ=255)                    tGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAACGGGACAg                                      >  1:4711906/1‑100 (MQ=255)                    tGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAAGGGACAg                                      >  1:801682/1‑100 (MQ=255)                     ggTGATATTCGTGAAGCTGCTGGCGAAAAAACATCTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGACaag                                     >  1:403442/1‑98 (MQ=255)                      gTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGaacaga                                    >  1:957242/1‑95 (MQ=255)                          tgcccGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAAT‑‑AAAAAAATAATCAATATATTGTTAGGAACGGGACAGatcggaa                               >  1:5131189/5‑94 (MQ=255)                           atTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGACaagatcgga                               >  1:5045489/1‑92 (MQ=255)                           atTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAGGGGACAGatcggaa                               >  1:3796232/1‑94 (MQ=255)                           atTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAACAGGACAGatcggaa                               >  1:565734/1‑94 (MQ=255)                             tCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGACAGatcggaaga                             >  1:175305/1‑92 (MQ=255)                              cGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGAAAGatcgaaagag                            >  1:6424267/1‑91 (MQ=255)                              cGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAAGGGAACAGatcggaagag                            >  1:4821392/1‑91 (MQ=255)                              cGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAACGGGACAGatcggaagag                            >  1:3673473/1‑91 (MQ=255)                                                        tAATCCAGTGGAAGATCGTTTCGTAATA‑AAAAAAATAATCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATgtaat  <  1:1549589/100‑3 (MQ=255)                                                                                   |                                                                          TTATTGCTCAAATGCCATTTGGTGATATTCGTGAAGCTGCTGGCGAAAAAACATTTAATCCAGTGGAAGATCGTTTCGTAAT‑‑AAAAAAATAATCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGT  >  NC_007793/2144738‑2144891

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