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%	121.8 / NA	39	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 (21/18);  total (21/18)

TAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTCA  >  NC_000913/1742441‑1742697                                                                                                                                       |                                                                                                                           tAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGACAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                        <  2:525728/139‑1 (MQ=255) tAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                        <  1:123894/139‑1 (MQ=255)         tGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgg                                                                                                                >  2:794760/1‑139 (MQ=255)         tGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgg                                                                                                                >  2:77794/1‑139 (MQ=255)         tGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgg                                                                                                                >  1:101796/1‑139 (MQ=255)         tGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgg                                                                                                                >  1:791821/1‑139 (MQ=255)                    tatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCa                                                                                                     <  1:794760/139‑1 (MQ=255)                    tatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCa                                                                                                     <  1:722817/139‑1 (MQ=255)                    tatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCa                                                                                                     >  1:256378/1‑139 (MQ=255)                     atCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgggg                                                                                                              <  1:403762/129‑1 (MQ=255)                     atCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgggg                                                                                                              >  2:403762/1‑129 (MQ=255)                     atCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCATTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAg                                                                                                    >  1:585681/1‑139 (MQ=255)                      tCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGa                                                                                                   <  1:652558/139‑1 (MQ=255)                           aTGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAgg                                                                                              >  1:414511/1‑139 (MQ=255)                            tGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGc                                                                                             >  1:380932/1‑139 (MQ=255)                              ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                           <  1:40029/139‑1 (MQ=255)                                            ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGa                                                                             >  2:626246/1‑139 (MQ=255)                                            ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGa                                                                             >  2:410488/1‑139 (MQ=255)                                               cgcgCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCGCGCGGTGTTGAAAACGGCCTTCGAGTGGGTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTc                                                                          >  1:172885/1‑139 (MQ=255)                                                 cgcCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCNGNGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCAt                                                                        >  2:777210/1‑139 (MQ=255)                                                       tGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCAtt                                                                  >  2:400298/1‑139 (MQ=255)                                                         atatTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCAtttt                                                                <  2:565318/139‑1 (MQ=255)                                                              cAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGt                                                           <  1:84301/139‑1 (MQ=255)                                                                      gCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAg                                                   >  2:108159/1‑139 (MQ=255)                                                                              tgtTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                                           <  1:108159/139‑1 (MQ=255)                                                                                 tctcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                                        <  2:414511/139‑1 (MQ=255)                                                                                 tctcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                                        <  1:77794/139‑1 (MQ=255)                                                                                   tcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc                                      <  2:172885/139‑1 (MQ=255)                                                                                               gAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcc                                                                                                                          <  2:315459/43‑1 (MQ=255)                                                                                               gAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACt                          >  2:760291/1‑139 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgggg                                                                                                              <  2:159724/49‑1 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgggg                                                                                                              >  1:159724/1‑49 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGa                    >  2:30667/1‑139 (MQ=255)                                                                                                       ccTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGt                  >  1:706682/1‑139 (MQ=255)                                                                                                                     gCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGAttt    <  1:410488/139‑1 (MQ=255)                                                                                                                      cTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTc   >  1:707615/1‑139 (MQ=255)                                                                                                                       tAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTCa  <  2:714905/139‑1 (MQ=255)                                                                                                                                ttCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGAt          <  2:388670/122‑1 (MQ=255)                                                                                                                                ttCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGAt          >  1:388670/1‑122 (MQ=255)                                                                                                                                       |                                                                                                                           TAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTCA  >  NC_000913/1742441‑1742697

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