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%	159.7 / NA	50	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 (25/25);  total (25/25)

TAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTCAATGTTGACG  >  NC_000913/1742441‑1742706                                                                                                                                       |                                                                                                                                    tAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                                                 <  1:755836/139‑1 (MQ=255)             tACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGc                                                                                                                     >  2:1405403/1‑139 (MQ=255)                  attatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAttt                                                                                                                <  2:84167/139‑1 (MQ=255)                         gAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAttt                                                                                                                >  1:888400/1‑132 (MQ=255)                         gAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAttt                                                                                                                <  2:888400/132‑1 (MQ=255)                         gAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCa                                                                                                         <  2:1684996/139‑1 (MQ=255)                          aaTGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAg                                                                                                        >  2:402882/1‑139 (MQ=255)                              ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                                    <  2:1713102/139‑1 (MQ=255)                                 gtgCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCtg                                                                                                 >  1:110873/1‑139 (MQ=255)                                  tgCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCtgt                                                                                                <  2:110873/139‑1 (MQ=255)                                       gTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTg                                                                                           >  1:955258/1‑139 (MQ=255)                                       gTGCTTTCCGCACCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAt                                                                                                                  <  2:362481/116‑1 (MQ=255)                                       gTGCTTTCCGCACCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAt                                                                                                                  >  1:362481/1‑116 (MQ=255)                                         gCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGt                                                                                         <  1:1405403/139‑1 (MQ=255)                                          cTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcc                                                                                                                                   >  2:42945/1‑96 (MQ=255)                                            ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGa                                                                                      >  1:1755942/1‑139 (MQ=255)                                             tCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                                    <  1:1606726/124‑1 (MQ=255)                                             tCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                                                    >  2:1606726/1‑124 (MQ=255)                                                gcgcCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCa                                                                                  <  2:955258/139‑1 (MQ=255)                                                       tGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCAtt                                                                           >  1:758250/1‑139 (MQ=255)                                                                gCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                              <  1:1260495/111‑1 (MQ=255)                                                                gCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                              >  2:1260495/1‑111 (MQ=255)                                                                 cttttgCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAtt                                                                                                                 >  1:1194124/6‑91 (MQ=255)                                                                 ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAtt                                                                                                                 <  2:1194124/91‑1 (MQ=255)                                                                 ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                              <  2:1353935/110‑1 (MQ=255)                                                                 ctctGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                              >  1:1353935/1‑110 (MQ=255)                                                                           tCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCacac                                                       >  2:29991/1‑139 (MQ=255)                                                                                    cACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc                                               >  1:1938819/1‑138 (MQ=255)                                                                                    cACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTAcc                                               <  2:1938819/138‑1 (MQ=255)                                                                                     aCGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGACAGCACACGTTCTACCGt                                             >  2:1383972/1‑139 (MQ=255)                                                                                       gcgGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTAt                                           <  2:1755942/139‑1 (MQ=255)                                                                                           tgtTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTAc                                       <  1:402882/139‑1 (MQ=255)                                                                                             ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAttt                                                                                                                <  2:1129897/64‑1 (MQ=255)                                                                                             ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGAttt                                                                                                                >  1:1129897/1‑64 (MQ=255)                                                                                             ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCa                                     >  2:666740/1‑139 (MQ=255)                                                                                                     ggtcTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGa                             >  2:1126758/4‑139 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                              >  1:801130/1‑74 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                                              <  2:801130/74‑1 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcaca                                                        >  1:1889529/1‑112 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcaca                                                        <  2:1889529/112‑1 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTCCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                                                 >  2:1973296/1‑119 (MQ=255)                                                                                                     ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTCCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                                                 <  1:1973296/119‑1 (MQ=255)                                                                                                       ccTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCg                                                                                                                          >  1:1097221/1‑44 (MQ=255)                                                                                                       ccTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCg                                                                                                                          <  2:1097221/44‑1 (MQ=255)                                                                                                               tGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGAt                   <  1:1126758/139‑1 (MQ=255)                                                                                                                     gCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGAttt             <  1:29991/139‑1 (MQ=255)                                                                                                                     gCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGCTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGAttt             <  2:70475/139‑1 (MQ=255)                                                                                                                                ttCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTCAATGTTGAcg  >  2:573897/1‑139 (MQ=255)                                                                                                                                 tCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcac                                                         >  1:1525538/1‑83 (MQ=255)                                                                                                                                 tCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcac                                                         <  2:1525538/83‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    TAAACGAATGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGAGTTTGTGGATCGATTTCAATGTTGACG  >  NC_000913/1742441‑1742706

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