Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,056,324	T→C	Y28Y (TAT→TAC)	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,324	0	T	C	100.0%	114.0 / NA	34	Y28Y (TAT→TAC)	fau	putative 5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/24);  total (10/24)

CGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAG  >  NC_000913/3056191‑3056455                                                                                                                                      |                                                                                                                                    cGCAACTACCAGAACCCCCACTCACCTTATCCCGACATGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCg                                                                                                                               <  1:193869/140‑1 (MQ=255)    aaCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCCACAAGCCACTACCCTGATGATGACTTACCCCCCggtg                                                                                                                            >  1:237219/1‑140 (MQ=255)     aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCggtgg                                                                                                                           <  2:138745/140‑1 (MQ=255)       tACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGgtggtg                                                                                                                         <  1:2337/140‑1 (MQ=255)       tACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGgtggtg                                                                                                                         <  1:79709/140‑1 (MQ=255)                        aCATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGc                                                                                                        >  2:246010/1‑140 (MQ=255)                              tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTAtt                                                                                                  >  1:110345/1‑140 (MQ=255)                              tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTAtt                                                                                                  <  1:30168/140‑1 (MQ=255)                                 cGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCct                                                                                               <  1:385645/140‑1 (MQ=255)                                  gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctc                                                                                              <  1:366974/140‑1 (MQ=255)                                   aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCGCCGGTGGGGATGGCGCATACGATCGCTGTATTCctct                                                                                             >  1:183123/1‑140 (MQ=255)                                    cAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctctc                                                                                            >  2:479332/1‑140 (MQ=255)                                               cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCg                                                                                 <  1:511307/140‑1 (MQ=255)                                                      aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGa                                                                          <  1:639128/140‑1 (MQ=255)                                                      aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGa                                                                          <  1:373246/140‑1 (MQ=255)                                                        tCGGCAACGTCGTCGTGTGTTACCGCCGGAACAACGGCAGGAAATGGGGCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGaca                                                                        <  2:283827/140‑1 (MQ=255)                                                             aaCGTCGTCGTGCGTTACCGCCAGAACAACAACAGCAAACGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAg                                                                   <  2:189901/140‑1 (MQ=255)                                                              aCGTCGTCGTGCGTTAAGGCCGAAACAACAGGAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc                                                                  <  2:270387/140‑1 (MQ=255)                                                               cgtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGcc                                                                 <  1:334171/140‑1 (MQ=255)                                                                    tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCa                                                            <  2:1121/140‑1 (MQ=255)                                                                    ccgtGCGTTAACGCCGGAACAACAGAAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCa                                                            <  2:452724/139‑1 (MQ=255)                                                                          gTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACt                                                              <  2:415460/132‑1 (MQ=255)                                                                          gTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACt                                                              >  1:415460/1‑132 (MQ=255)                                                                                            gcagGAAATGGGTCCACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                    <  1:624261/140‑1 (MQ=255)                                                                                            gcagGAAATGGGTCAACGAGCCTCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                    <  2:286055/140‑1 (MQ=255)                                                                                            gcagGAAATGGGTCAACAAGTCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                    <  1:377649/140‑1 (MQ=255)                                                                                            gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                    <  2:110345/140‑1 (MQ=255)                                                                                            gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                    <  2:392406/140‑1 (MQ=255)                                                                                                 aaaTGGGTCAACAAGCCGCTACCCGCATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcgc                               <  1:246010/140‑1 (MQ=255)                                                                                                       gTCAACAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTatat                         >  1:598589/1‑140 (MQ=255)                                                                                                                ccGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGtt                >  2:439053/1‑140 (MQ=255)                                                                                                                 cGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGttt               >  2:549861/1‑140 (MQ=255)                                                                                                                   cTACCCGGATGATGACTCACCCCCCGGTGGTGATGGCACCTACGGTCGTTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg             <  2:279677/140‑1 (MQ=255)                                                                                                                              atgaCTTACCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAg  >  2:405455/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                    CGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAG  >  NC_000913/3056191‑3056455

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