Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	69,513	G→C	Q103E (CAA→GAA)	SAUSA300_RS00310 ←	carbamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	69,513	0	G	C	100.0%	147.6 / NA	53	Q103E (CAA→GAA)	SAUSA300_RS00310	carbamate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (6/47);  total (6/47)

TTTTGTATAAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTGTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAGGGAAAGG  >  NC_007793/69421‑69609                                                                                             |                                                                                                 ttttGTATAAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCt                                                                                           <  1:5025876/100‑1 (MQ=255)      tataAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGcc                                                                                      <  1:1292181/100‑1 (MQ=255)       ataAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGccc                                                                                     <  1:1764368/100‑1 (MQ=255)       ataAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCAGCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGccc                                                                                     <  1:5251267/100‑1 (MQ=255)        taAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCa                                                                                    <  1:4098527/100‑1 (MQ=255)         aaaaTAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCag                                                                                   <  1:784127/100‑2 (MQ=255)          gatctgtCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCAtt                                                                                  <  1:4827464/94‑1 (MQ=255)          aaaTAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGATCTTGTGTTACGAGCGTAACAACTTCTTTATCTATGCCCaca                                                                                  <  1:4932331/100‑3 (MQ=255)          aaaTAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCAtt                                                                                  <  1:985555/100‑1 (MQ=255)          aaaTAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCAtt                                                                                  <  1:1467382/100‑1 (MQ=255)          aaaTAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCAtg                                                                                  <  1:3997160/100‑2 (MQ=255)          aaaTAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCAtc                                                                                  <  1:2331468/100‑2 (MQ=255)               atCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATtgaatg                                                                             <  1:1204202/100‑1 (MQ=255)               atCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATAGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATtgaagt                                                                             <  1:5587791/100‑3 (MQ=255)                tCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATtgaatga                                                                            <  1:210925/100‑1 (MQ=255)                tCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATtgaatga                                                                            <  1:4838002/100‑1 (MQ=255)                   aaTTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGt                                                                         >  1:5676456/1‑100 (MQ=255)                      tGGTTTAGTAGGGTTACTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCg                                                                      >  1:1637653/1‑100 (MQ=255)                           tagtagGGTTATTGAAAGCAGAATCATCACTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTg                                                                 >  1:4705901/1‑100 (MQ=255)                              cagGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTaa                                                              >  1:847930/2‑100 (MQ=255)                                      ttGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCtt                                                      <  1:4025390/100‑1 (MQ=255)                                        gAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGc                                                    <  1:614808/100‑1 (MQ=255)                                            gCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCt                                                <  1:1684195/100‑1 (MQ=255)                                             cAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTg                                               <  1:47385/100‑1 (MQ=255)                                              aGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTgg                                              <  1:4068903/100‑2 (MQ=255)                                                 atcatcGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAg                                           <  1:5224083/100‑1 (MQ=255)                                                   catcGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGcc                                         <  1:2275264/100‑1 (MQ=255)                                                    atcGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCa                                        <  1:5018965/100‑1 (MQ=255)                                                     tcGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCaa                                       <  1:3350773/100‑1 (MQ=255)                                                     gatcTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCaa                                       <  1:1730372/97‑1 (MQ=255)                                                        cTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGcggggac                                    <  1:4040523/100‑7 (MQ=255)                                                          tGCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATaa                                  >  1:2917067/1‑100 (MQ=255)                                                           gCAACTTGAACTTGTGTTACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCaagggaag                                 <  1:2315855/100‑7 (MQ=255)                                                                           ttACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGcacc                 <  1:3879288/100‑1 (MQ=255)                                                                            tACTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGcacca                <  1:4061694/100‑1 (MQ=255)                                                                             aCTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGctccaa               <  1:3157016/100‑2 (MQ=255)                                                                             aCTAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGcaccag               <  1:971643/100‑2 (MQ=255)                                                                               tAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACcacag             <  1:4116992/100‑2 (MQ=255)                                                                               tAGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACAt             <  1:582202/100‑1 (MQ=255)                                                                                aGCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACAtt            <  1:2478533/100‑1 (MQ=255)                                                                                 gCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTc           <  1:5214521/100‑1 (MQ=255)                                                                                 gCGTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTc           <  1:2454959/100‑1 (MQ=255)                                                                                   gTAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAg         >  1:2662661/1‑100 (MQ=255)                                                                                    tAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAgg        <  1:5483627/100‑1 (MQ=255)                                                                                    tAACAACTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACATTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAgg        <  1:1404708/100‑1 (MQ=255)                                                                                     aacaacTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAggg       <  1:1941498/100‑1 (MQ=255)                                                                                     aacaacTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAggg       <  1:4899744/100‑1 (MQ=255)                                                                                     aacaacTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAggg       <  1:1745522/100‑1 (MQ=255)                                                                                     aacaacTTCTTTATCTATGACCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGATTGACTCATAGCACCACATTCAggg       <  1:4293122/100‑1 (MQ=255)                                                                                      acatcTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAggta      <  1:3400308/100‑3 (MQ=255)                                                                                      acaacTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAGgga      <  1:5301158/100‑1 (MQ=255)                                                                                      acaacTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAGgga      <  1:854174/100‑1 (MQ=255)                                                                                       caacTTCTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAGggaa     <  1:4398818/100‑1 (MQ=255)                                                                                          cttcttTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAGggaaagg  <  1:5261447/100‑1 (MQ=255)                                                                                             |                                                                                                 TTTTGTATAAAATAATCCAATTGGTTTAGTAGGGTTATTGAAAGCAGAATCATCGCTTGCAACTTGAACTTGTGTTACTAGCGTAACAACTTGTTTATCTATGCCCATTGAATGAAGTTCGTTTTGTAAACTTTCTTGCATCTGATAGCCAATATAAGCTTGACTCATAGCACCACATTCAGGGAAAGG  >  NC_007793/69421‑69609

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