Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,125,592	A→T	I49L (ATA→TTA)	SAUSA300_RS10825 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,125,592	0	A	T	100.0%	137.4 / NA	47	I49L (ATA→TTA)	SAUSA300_RS10825	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (44/3);  total (44/3)

TTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATC‑AAAAAACGTCTACAAGTGTAGACGTTGA  >  NC_007793/2125499‑2125669                                                                                              |                                                                               ttAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATg                                                                          >  1:3946617/1‑100 (MQ=255)  tAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATgg                                                                         >  1:1545788/1‑100 (MQ=255)  tAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATgg                                                                         >  1:91423/1‑100 (MQ=255)  tAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATgg                                                                         >  1:611212/1‑100 (MQ=255)  tAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATgg                                                                         >  1:2724872/1‑100 (MQ=255)   aaCTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGc                                                                        >  1:4008131/1‑100 (MQ=255)   aaCTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGc                                                                        >  1:5096061/1‑100 (MQ=255)   aaCTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGc                                                                        <  1:4942135/100‑1 (MQ=255)   aaCTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGc                                                                        >  1:1642673/1‑100 (MQ=255)       ctggccaaTCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATa                                                                    >  1:3924970/7‑100 (MQ=255)       aGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATa                                                                    >  1:259872/1‑100 (MQ=255)       aGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATa                                                                    >  1:2870462/1‑100 (MQ=255)       aGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATa                                                                    >  1:4545307/1‑100 (MQ=255)       aGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATa                                                                    >  1:3398410/1‑100 (MQ=255)             aaTCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCGCATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTaa                                                              >  1:5110772/1‑100 (MQ=255)              tgtgcccATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAAc                                                             >  1:4092325/6‑100 (MQ=255)                cAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAGTACTTGGAATTTATTTCACAAAAATTAGATGGCGCTAAGTTAACCa                                                           >  1:825004/1‑100 (MQ=255)                cAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAACTTAGATGGCGATAAGTTAACCa                                                           >  1:3178761/1‑100 (MQ=255)                cAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCa                                                           >  1:3344977/1‑100 (MQ=255)                 gaCCATTCACTTTTGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCaa                                                          >  1:5624240/2‑100 (MQ=255)                    ccttCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAaaga                                                       >  1:2029173/3‑100 (MQ=255)                    cATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAaaga                                                       <  1:639987/100‑1 (MQ=255)                         aCTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGAtt                                                  >  1:2921946/1‑100 (MQ=255)                             aTGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCAtt                                              >  1:4714742/1‑100 (MQ=255)                                gcttACCTAAAGATTTATCGTCACATAATGCGCGTAAGTACTTGGAATTTATTTCACCGAAATTAGAGGGGGATAAGTTAACCAAAGAAGATTCATTATg                                           >  1:3729702/3‑100 (MQ=255)                                 gTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATCTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGa                                          >  1:4758260/1‑100 (MQ=255)                                   tACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGAtt                                        >  1:3992069/1‑100 (MQ=255)                                     ccTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGAtttt                                      >  1:4324813/1‑100 (MQ=255)                                              ttATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aa                            >  1:4003007/1‑100 (MQ=255)                                              ttATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aa                            >  1:5279555/1‑100 (MQ=255)                                               gcggggcacaTAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGCTTTTACTAAAC‑aaa                           >  1:2717669/7‑100 (MQ=255)                                                 acGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaa                         >  1:2536742/2‑100 (MQ=255)                                                  cGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaa                        >  1:897557/1‑100 (MQ=255)                                                     cacaTAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacg                     >  1:4551210/1‑97 (MQ=255)                                                     cacaTAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacg                     >  1:2280214/1‑97 (MQ=255)                                                        gtAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcct                  >  1:3240433/2‑94 (MQ=255)                                                          aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctac                >  1:412968/1‑92 (MQ=255)                                                          aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctac                >  1:3889151/1‑92 (MQ=255)                                                          aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcccag                >  1:4519552/1‑92 (MQ=255)                                                          aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTAct                >  1:3470874/1‑99 (MQ=255)                                                           ctGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctagg               <  1:4191690/99‑10 (MQ=255)                                                            tGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctacta              >  1:780784/1‑90 (MQ=255)                                                             gcgcGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTACTAgt             >  1:1873586/1‑100 (MQ=255)                                                                 gTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAACAAAAAAACGCCTACTAgtgtg          >  1:1192966/1‑99 (MQ=255)                                                                   caaTACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAACAAAAAAACGCCTACTAgtgtgaa        >  1:5339638/2‑97 (MQ=255)                                                                      caCTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTACTAgtgtgaaaacg    >  1:5324164/2‑93 (MQ=255)                                                                        cTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTACTAgtgtgaaaacgta  >  1:225185/1‑91 (MQ=255)                                                                                              |                                                                               TTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATC‑AAAAAACGTCTACAAGTGTAGACGTTGA  >  NC_007793/2125499‑2125669

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