Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,019,080	C→A	intergenic (‑3/+128)	SAUSA300_RS10140 ← / ← SAUSA300_RS10145	glutamine amidotransferase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,019,080	0	C	A	100.0%	156.0 / NA	54	intergenic (‑3/+128)	SAUSA300_RS10140/SAUSA300_RS10145	glutamine amidotransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (2/52);  total (2/52)

CAGTATTAAAGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGCAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACAAGATATAAAAAATATTA  >  NC_007793/2018981‑2019163                                                                                                    |                                                                                    cAGTATTAAAGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGc                                                                                     <  1:973097/100‑1 (MQ=255)         aaGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCAcg                                                                             <  1:4641919/100‑2 (MQ=255)          aGCCTGCATTCTTTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCt                                                                            <  1:616067/100‑1 (MQ=255)          aGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCt                                                                            <  1:1501010/100‑1 (MQ=255)          aGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCt                                                                            <  1:5260431/100‑1 (MQ=255)          aGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCt                                                                            <  1:5133425/100‑1 (MQ=255)          aGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCt                                                                            <  1:4665313/100‑1 (MQ=255)          aGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCt                                                                            <  1:3490215/100‑1 (MQ=255)           gCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCtt                                                                           <  1:1885422/100‑1 (MQ=255)            ccTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCttg                                                                          <  1:6608203/100‑2 (MQ=255)            ccTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCtgt                                                                          <  1:5365825/100‑3 (MQ=255)             cTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCtttg                                                                         <  1:6200615/100‑2 (MQ=255)                cATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCtat                                                                      <  1:6083572/100‑1 (MQ=255)                cATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCtat                                                                      <  1:3947463/100‑1 (MQ=255)                cATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCtat                                                                      <  1:2970328/100‑1 (MQ=255)                  ttCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCtatat                                                                    <  1:3041853/100‑1 (MQ=255)                     ctcttccgatctTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATc                                                                 <  1:1984142/91‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:3497107/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:5983172/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:573114/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:2809091/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:3355592/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:4731761/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:3971298/100‑1 (MQ=255)                       tAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAt                                                               <  1:3775965/100‑1 (MQ=255)                        aaTGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCtagatcggaag                                                              >  1:654477/1‑94 (MQ=255)                        aaTGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAtt                                                              <  1:4872850/100‑1 (MQ=255)                        aaTGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCAtt                                                              <  1:595795/100‑1 (MQ=255)                          tGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTgg                                                            <  1:1875070/100‑2 (MQ=255)                           gCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGat                                                           <  1:6458864/100‑1 (MQ=255)                           gCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGat                                                           <  1:6523962/100‑1 (MQ=255)                            ccTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGatt                                                          <  1:5572327/100‑2 (MQ=255)                            ccTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGatg                                                          <  1:1032523/100‑2 (MQ=255)                            ccTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGatg                                                          <  1:6503074/100‑2 (MQ=255)                            ccTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGatc                                                          <  1:5142373/100‑2 (MQ=255)                            ccTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGatc                                                          <  1:583030/100‑2 (MQ=255)                             ctctTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGattt                                                         <  1:5233494/100‑3 (MQ=255)                             ctctTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGatga                                                         <  1:6355052/100‑3 (MQ=255)                                    gggCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACt                                                  <  1:1636544/100‑1 (MQ=255)                                      gCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCTCCTTTCTATATATCATTGATATAATTACtat                                                <  1:2194/100‑1 (MQ=255)                                            gctcttccgatCTTTAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACag                                          <  1:4898172/91‑2 (MQ=255)                                              atTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGt                                        <  1:635585/100‑1 (MQ=255)                                                   atatCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATc                                   <  1:4957027/100‑1 (MQ=255)                                                   atatCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATc                                   <  1:5359212/100‑1 (MQ=255)                                                      tCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGt                                <  1:4519831/100‑1 (MQ=255)                                                                  tttGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATAc                    <  1:6299645/100‑1 (MQ=255)                                                                   ttGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACa                   <  1:5146102/100‑1 (MQ=255)                                                                   ttGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACa                   <  1:509779/100‑1 (MQ=255)                                                                   ttGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACa                   <  1:2980706/100‑1 (MQ=255)                                                                   ttGCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACa                   <  1:2313471/100‑1 (MQ=255)                                                                   ttGCTAGAATAATTGCTACTTTTTAAGTCATTGAAAATCACCTTTCTATAGATCATGGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACa                   <  1:1053017/100‑1 (MQ=255)                                                                   cttccgatctaatTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGTCAAGTAAATCAGTGATTAAACATACa                   <  1:6075211/91‑1 (MQ=255)                                                                     gCTAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACAAg                 <  1:4130099/100‑1 (MQ=255)                                                                       tAGAATAATTGCTACTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATGCAAGat               <  1:6066994/100‑1 (MQ=255)                                                                                    aCTTTTTTAGTCATTGAAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACAAGATATAAAAAATATTa  >  1:3715764/1‑100 (MQ=255)                                                                                                    |                                                                                    CAGTATTAAAGCCTGCATTCTCTAATGCCTCTTTAGGGCTTGAATATTCTATATCTTCAAATTCGTTTGCTAGAATAATTGCTACTTTTTTAGTCATTGCAAATCACCTTTCTATATATCATTGATATAATTACTATAGACAAGTAAATCAGTGATTAAACATACAAGATATAAAAAATATTA  >  NC_007793/2018981‑2019163

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