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%	203.8 / NA	63	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 (30/33);  total (30/33)

TTTAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTT  >  NC_000913/1742439‑1742695                                                                                                                                         |                                                                                                                         tttAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcc                                                                                                                        >  2:316711/1‑139 (MQ=255)     aacattcgTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCct                                                                                                                    >  1:1293463/8‑139 (MQ=255)           tGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCGGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgg                                                                                                              >  1:1422989/1‑139 (MQ=255)              ttACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCggggg                                                                                                           <  2:835536/139‑1 (MQ=255)               tACCTATTATCTGAATGCCTGTGCAGGTGCTTTCGGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGc                                                                                                          <  2:1422989/139‑1 (MQ=255)                    attatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                      <  2:746421/122‑1 (MQ=255)                    attatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                      >  1:746421/1‑122 (MQ=255)                       atCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAg                                                                                                  >  2:460339/1‑139 (MQ=255)                       atCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAg                                                                                                  >  1:1466886/1‑139 (MQ=255)                                ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                         <  2:1293463/139‑1 (MQ=255)                                ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                         <  1:316711/139‑1 (MQ=255)                                ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                         <  1:1678239/139‑1 (MQ=255)                                    tgCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcc                                                                                                                        >  2:817394/1‑104 (MQ=255)                                     gCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCtgtg                                                                                    <  1:460339/139‑1 (MQ=255)                                      cAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                   >  1:1720716/1‑139 (MQ=255)                                        ggTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCt                                                                                 <  1:1291166/139‑1 (MQ=255)                                         gTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                      >  1:734631/1‑101 (MQ=255)                                         gTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                      <  2:734631/101‑1 (MQ=255)                                         gTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTg                                                                                >  2:1328418/1‑139 (MQ=255)                                         gTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTg                                                                                >  1:494852/1‑139 (MQ=255)                                            cTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGtt                                                                             >  1:1430013/1‑139 (MQ=255)                                              ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCggggg                                                                                                           <  1:514581/107‑1 (MQ=255)                                              ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCggggg                                                                                                           >  2:514581/1‑107 (MQ=255)                                              ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGc                                                                                                          <  1:224668/108‑1 (MQ=255)                                              ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGc                                                                                                          >  2:224668/1‑108 (MQ=255)                                              ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGa                                                                           >  2:1013058/1‑139 (MQ=255)                                              ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGa                                                                           >  1:1269834/1‑139 (MQ=255)                                                  gcgcCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCa                                                                       <  1:1328418/139‑1 (MQ=255)                                                       cGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCa                                                                  >  1:674416/1‑139 (MQ=255)                                                         tGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                   <  2:298076/120‑1 (MQ=255)                                                         tGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                   >  1:298076/1‑120 (MQ=255)                                                                cAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGAtt                                                                         >  1:61165/1‑123 (MQ=255)                                                                cAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGAtt                                                                         <  2:61165/123‑1 (MQ=255)                                                                   ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcc                                                                                                                        <  2:1156185/73‑1 (MQ=255)                                                                   ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgcc                                                      >  1:927467/1‑139 (MQ=255)                                                                        gCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAg                                                 >  2:1432187/1‑139 (MQ=255)                                                                        gCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAg                                                 >  1:1280122/1‑139 (MQ=255)                                                                          aGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGca                                               <  2:1720716/139‑1 (MQ=255)                                                                           ggTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcac                                              <  2:1280122/139‑1 (MQ=255)                                                                           ggTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcac                                              <  1:1013058/139‑1 (MQ=255)                                                                                tgtTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                                         <  2:233429/139‑1 (MQ=255)                                                                                  ttCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                      >  2:533767/1‑60 (MQ=255)                                                                                  ttCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                      <  1:533767/60‑1 (MQ=255)                                                                                     tcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc                                    <  2:494852/139‑1 (MQ=255)                                                                                         gcgGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTAt                                <  2:927467/139‑1 (MQ=255)                                                                                         gcgGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTAt                                <  2:1466886/139‑1 (MQ=255)                                                                                            gtgtTGAAAACGGCCTTCGAGTGACTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                                         <  1:541990/127‑1 (MQ=255)                                                                                            gtgtTGAAAACGGCCTTCGAGTGACTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                                         >  2:541990/1‑127 (MQ=255)                                                                                               ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCa                          <  1:1190905/139‑1 (MQ=255)                                                                                               ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCa                          <  2:1269834/139‑1 (MQ=255)                                                                                                    aaCGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCTGCCAGCACACGTTCTACCGTATCTACCACTGcc                     <  2:674416/139‑1 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGTCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCAt                                                                      <  1:1158912/87‑1 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCAt                                                                      >  2:1158912/1‑87 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGt                                          <  2:147838/115‑1 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGt                                          >  1:147838/1‑115 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                                      <  1:848951/119‑1 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                                      >  2:848951/1‑119 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGa                  >  1:311059/1‑139 (MQ=255)                                                                                                       ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGa                  >  2:1638716/1‑139 (MQ=255)                                                                                                          cTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGtt               <  2:1430013/139‑1 (MQ=255)                                                                                                                       gCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGAttt  <  2:311059/139‑1 (MQ=255)                                                                                                                        cTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgc                                                       >  2:1123683/1‑85 (MQ=255)                                                                                                                        cTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgc                                                       <  1:1123683/85‑1 (MQ=255)                                                                                                                                         |                                                                                                                         TTTAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTT  >  NC_000913/1742439‑1742695

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