Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,056,402	A→G	E54E (GAA→GAG)	fau →	putative 5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,056,402	0	A	G	100.0%	145.2 / NA	41	E54E (GAA→GAG)	fau	putative 5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/23);  total (18/23)

CGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGG  >  NC_000913/3056272‑3056538                                                                                                                                   |                                                                                                                                         cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTgg                                                                                                                                 >  1:777102/1‑140 (MQ=255)           gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTa                                                                                                                       <  2:146641/140‑1 (MQ=255)           gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTa                                                                                                                       <  1:309905/140‑1 (MQ=255)            cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTTTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTaa                                                                                                                      <  2:721940/140‑1 (MQ=255)             agGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAg                                                                                                                     <  2:9635/140‑1 (MQ=255)               gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAgcg                                                                                                                   <  1:221171/140‑1 (MQ=255)               gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAgcg                                                                                                                   <  2:777102/140‑1 (MQ=255)                          acaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                                        >  1:814009/1‑140 (MQ=255)                             aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAg                                                                                                     <  1:330688/140‑1 (MQ=255)                                  cTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                                                <  2:598299/140‑1 (MQ=255)                                    aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCa                                                                                              <  2:619686/140‑1 (MQ=255)                                    aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCa                                                                                              <  2:860832/140‑1 (MQ=255)                                               gaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTg                                                                                   >  1:336146/1‑140 (MQ=255)                                                    aTCCCCCGGTGGTGTTGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCa                                                                                              <  2:807610/124‑1 (MQ=255)                                                    aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCa                                                                                              >  1:807610/1‑124 (MQ=255)                                                    aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGc                                                                                  >  1:180658/1‑136 (MQ=255)                                                    aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGc                                                                                  <  2:180658/136‑1 (MQ=255)                                                    aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGt                                                                              >  2:172502/1‑140 (MQ=255)                                                         ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAttt                                                                         >  2:116996/1‑140 (MQ=255)                                                              ggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGt                                                                    <  2:814009/140‑1 (MQ=255)                                                                tgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTc                                                                  >  1:940673/1‑140 (MQ=255)                                                                         aTACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTAc                                                         >  2:1026007/1‑140 (MQ=255)                                                                               tCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCg                                                   >  2:62622/1‑140 (MQ=255)                                                                                  cTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCaa                                                >  2:472913/1‑140 (MQ=255)                                                                                     tATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGc                                             >  1:920245/1‑140 (MQ=255)                                                                                         ccTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAAc                                         <  1:116996/140‑1 (MQ=255)                                                                                         ccTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAAc                                         >  2:807892/1‑140 (MQ=255)                                                                                         ccTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAAc                                         <  1:472913/140‑1 (MQ=255)                                                                                             tctTTTGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                                     <  2:148073/140‑1 (MQ=255)                                                                                               ttttGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGaa                                          <  2:358788/133‑1 (MQ=255)                                                                                               ttttGATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGaa                                          >  1:358788/1‑133 (MQ=255)                                                                                                   gATGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGaa                               <  1:1026007/140‑1 (MQ=255)                                                                                                    aTGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                              >  2:363184/1‑140 (MQ=255)                                                                                                    aTGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                              >  2:4292/1‑140 (MQ=255)                                                                                                    aTGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                              <  2:336146/140‑1 (MQ=255)                                                                                                    aTGGCGAACTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                              >  1:950579/1‑140 (MQ=255)                                                                                                           aCTCGACACCCAGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGa                       <  2:367438/140‑1 (MQ=255)                                                                                                                      aGCCACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAg            >  2:363663/1‑140 (MQ=255)                                                                                                                        ccACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGcc          <  2:310526/140‑1 (MQ=255)                                                                                                                        ccACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGcc          <  1:4292/140‑1 (MQ=255)                                                                                                                        ccACTCATAGAGCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGcc          <  1:170838/140‑1 (MQ=255)                                                                                                                                agagCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTgg  <  1:888717/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                         CGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCCAAAATTGG  >  NC_000913/3056272‑3056538

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