Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,742,575	G→A	intergenic (+14/+26)	cfa → / ← ribC	cyclopropane fatty acyl phospholipid synthase, SAM‑dependent/riboflavin synthase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,742,575	0	G	A	100.0%	108.0 / NA	35	intergenic (+14/+26)	cfa/ribC	cyclopropane fatty acyl phospholipid synthase, SAM‑dependent/riboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/21);  total (14/21)

TTTAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCA  >  NC_000913/1742439‑1742671                                                                                                                                         |                                                                                                 tttAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcc                                                                                                >  2:227112/1‑139 (MQ=255)        gAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCtctc                                                                                         <  2:1025706/139‑1 (MQ=255)               tACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGc                                                                                  <  2:121089/139‑1 (MQ=255)                           gAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCa                                                                      >  1:123868/1‑139 (MQ=255)                            aaTGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAg                                                                     >  1:31894/1‑139 (MQ=255)                            aaTGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAg                                                                     >  1:342193/1‑139 (MQ=255)                            aaTGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAg                                                                     >  2:112365/1‑139 (MQ=255)                               gCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCt                                                                  >  2:573302/1‑139 (MQ=255)                                ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                 <  2:31894/139‑1 (MQ=255)                                  tgtgCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCt                                                               >  1:729345/1‑139 (MQ=255)                                      cAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                           <  2:123868/139‑1 (MQ=255)                                       aGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGcc                                                          <  2:1200059/139‑1 (MQ=255)                                            cTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                           >  1:1118917/1‑133 (MQ=255)                                            cTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                           <  2:1118917/133‑1 (MQ=255)                                                ccGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGAtt                                                 <  1:313416/139‑1 (MQ=255)                                                   cgcCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCa                                               <  2:556111/138‑1 (MQ=255)                                                   cgcCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCa                                               >  1:556111/1‑138 (MQ=255)                                                                cAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGt                                 <  1:419768/139‑1 (MQ=255)                                                                  gCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgc                               >  2:730673/1‑139 (MQ=255)                                                                   ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCtgtg                                                            >  1:892851/1‑109 (MQ=255)                                                                   ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCtgtg                                                            <  2:892851/109‑1 (MQ=255)                                                                   ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgcc                              >  2:1051437/1‑139 (MQ=255)                                                                   ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgcc                              <  2:769063/139‑1 (MQ=255)                                                                    tctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                 >  1:298769/1‑103 (MQ=255)                                                                    tctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                 <  2:298769/103‑1 (MQ=255)                                                                           ggTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcac                      <  1:245274/139‑1 (MQ=255)                                                                              cGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACg                   <  1:112365/139‑1 (MQ=255)                                                                                   tctcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa              <  1:573302/139‑1 (MQ=255)                                                                                    ctcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAc             <  2:342193/139‑1 (MQ=255)                                                                                     tcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc            <  1:1051437/139‑1 (MQ=255)                                                                                     tcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc            <  2:729345/139‑1 (MQ=255)                                                                                     tcACGCGGTGTTGAAAACGGCCTTCGAGTGACTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc            <  1:92371/139‑1 (MQ=255)                                                                                             tgtTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                              >  2:1365327/1‑49 (MQ=255)                                                                                             tgtTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                              <  1:1365327/49‑1 (MQ=255)                                                                                               ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCa  <  1:759267/139‑1 (MQ=255)                                                                                                                                         |                                                                                                 TTTAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCA  >  NC_000913/1742439‑1742671

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