Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	912,256	C→T	P66S (CCA→TCA)	SAUSA300_RS04495 →	protein NagD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	912,256	0	C	T	100.0%	168.6 / NA	58	P66S (CCA→TCA)	SAUSA300_RS04495	protein NagD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (54/4);  total (54/4)

TCCTCATTTATACGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAACCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAGTGGTTTAA  >  NC_007793/912171‑912352                                                                                      |                                                                                                 tCCTCATTTATACGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTa                                                                                    <  1:3634388/100‑1 (MQ=255) accTCATTTATACGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTa                                                                                    >  1:4443890/2‑100 (MQ=255)     cATTTATACGTAACTAATAATTCAACAAATACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGt                                                                                >  1:3099914/1‑100 (MQ=255)        ttATACGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAg                                                                             >  1:2679655/1‑100 (MQ=255)         tataCGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGc                                                                            >  1:4625637/1‑100 (MQ=255)          ataCGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCg                                                                           >  1:3756805/1‑100 (MQ=255)            aCGTAACTAATAATTCAACAAAAACACCTGAGCAAGCAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGCTGTAACGTCAGCGtt                                                                         >  1:97459/1‑100 (MQ=255)                 gcTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCa                                                                    >  1:977482/2‑100 (MQ=255)                          tCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGAtt                                                           >  1:3575608/1‑100 (MQ=255)                          ccaacaaAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGAtt                                                           >  1:3982888/2‑100 (MQ=255)                           caacaaAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATta                                                          >  1:148880/1‑100 (MQ=255)                           caacaaAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATta                                                          >  1:848491/1‑100 (MQ=255)                           caacaaAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATta                                                          >  1:480013/1‑100 (MQ=255)                           caacaaAAACACCTGAGCAAGTAACTAAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATta                                                          >  1:2190408/1‑100 (MQ=255)                           caacaaAAACACCTGAGCAAGAAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATta                                                          >  1:5457302/1‑100 (MQ=255)                               aaaaaCACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATAtt                                                      >  1:5579464/1‑100 (MQ=255)                                taaaCACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATAttt                                                     >  1:63497/2‑100 (MQ=255)                                aaaaCACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATAttt                                                     >  1:2388231/1‑100 (MQ=255)                                 caaCACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTc                                                    >  1:2073867/2‑100 (MQ=255)                                 aaaCACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTc                                                    >  1:947776/1‑100 (MQ=255)                                 aaaCACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTc                                                    >  1:1102275/1‑100 (MQ=255)                                  cccacCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCa                                                   >  1:4340449/3‑100 (MQ=255)                                   acacCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAg                                                  >  1:4854873/1‑100 (MQ=255)                                     acCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGaa                                                >  1:2396159/1‑100 (MQ=255)                                      gcacaccAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCGGAAc                                               >  1:4197984/7‑100 (MQ=255)                                         gAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAAt                                            >  1:5212696/1‑100 (MQ=255)                                             aaGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCAcc                                        >  1:3459360/1‑100 (MQ=255)                                                 aaCTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGa                                    >  1:3996886/1‑100 (MQ=255)                                                  aCTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAg                                   >  1:1885538/1‑100 (MQ=255)                                                   cTGAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGc                                  <  1:741937/100‑1 (MQ=255)                                                     taaaaaTTACGTGAAATGCACCTTGATGCTCAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCAt                                >  1:5513666/2‑100 (MQ=255)                                                     gAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAAGatcga                                >  1:3663260/1‑96 (MQ=255)                                                     gAAAAATTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACCGCTGATTATATTCCAGAACAATCACCAGGAGCAt                                >  1:5299437/1‑100 (MQ=255)                                                     gAAAAATTACGTGAAATGCACATTGACGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAgatc                                >  1:3731248/1‑97 (MQ=255)                                                        tcaTTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAg                             >  1:4601444/3‑100 (MQ=255)                                                         aaTTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAgg                            <  1:2012713/100‑2 (MQ=255)                                                         aaTTACGTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGCCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGt                            >  1:946519/1‑100 (MQ=255)                                                              agTGAAATGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGtatata                       >  1:2624090/2‑100 (MQ=255)                                                                  aaaTGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCGCCAGGAGCATCAGTATATATGtt                   <  1:509348/100‑1 (MQ=255)                                                                   caTGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTa                  >  1:4804833/2‑100 (MQ=255)                                                                   aaTGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTa                  >  1:4435856/1‑100 (MQ=255)                                                                    aTGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAg                 >  1:1770870/1‑100 (MQ=255)                                                                    aTGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAg                 >  1:1893693/1‑100 (MQ=255)                                                                     tGCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAga                >  1:2709621/1‑99 (MQ=255)                                                                      gCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGt               >  1:1367946/1‑100 (MQ=255)                                                                      gCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGt               >  1:3414731/1‑100 (MQ=255)                                                                      gCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGt               >  1:605250/1‑100 (MQ=255)                                                                      gCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGt               >  1:2286617/1‑100 (MQ=255)                                                                      gCACATTGATGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGt               >  1:1237909/1‑100 (MQ=255)                                                                            gcaTGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGCGGGAga         >  1:3076200/3‑99 (MQ=255)                                                                              aTGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAGTgg       >  1:65695/1‑100 (MQ=255)                                                                               tGCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAgatcg      >  1:56581/1‑96 (MQ=255)                                                                               ggCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAGTGGt      >  1:4139665/2‑100 (MQ=255)                                                                                gCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGagagatcgg     >  1:2172427/1‑95 (MQ=255)                                                                                gCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAGTGGtt     >  1:1541969/1‑100 (MQ=255)                                                                                gCTAAATCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAGTGGtt     >  1:1770571/1‑100 (MQ=255)                                                                                   aaaTCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGATCGGaagagcacac  >  1:3583664/1‑91 (MQ=255)                                                                                   aaaTCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAAGTGGGAgatcggaag  >  1:2577463/1‑92 (MQ=255)                                                                                      |                                                                                                 TCCTCATTTATACGTAACTAATAATTCAACAAAAACACCTGAGCAAGTAACTGAAAAATTACGTGAAATGCACATTGATGCTAAACCAGAAGAGGTTGTAACGTCAGCGTTAGCCACTGCTGATTATATTTCAGAACAATCACCAGGAGCATCAGTATATATGTTAGGTGGGAGTGGTTTAA  >  NC_007793/912171‑912352

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: ‑