Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,056,273	G→A	P11P (CCG→CCA)	fau →	5‑formyltetrahydrofolate cyclo‑ligase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,056,273	0	G	A	100.0%	136.3 / NA	42	P11P (CCG→CCA)	fau	5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/19);  total (23/19)

GCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAA  >  NC_000913/3056141‑3056402                                                                                                                                     |                                                                                                                                  gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                             <  2:636246/139‑1 (MQ=255) gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                             <  1:1494168/139‑1 (MQ=255)                        tCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                                      <  2:1687176/139‑1 (MQ=255)                                        aCCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGatgatg                                                                                      >  1:1114949/1‑139 (MQ=255)                                           aGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACt                                                                                   >  1:1398333/1‑139 (MQ=255)                                                tGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                             >  2:672852/1‑92 (MQ=255)                                                tGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                             >  1:1388579/1‑92 (MQ=255)                                                tGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                             <  2:1388579/92‑1 (MQ=255)                                                tGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                             <  1:672852/92‑1 (MQ=255)                                                       aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaac                                                                                                                               <  1:140448/83‑1 (MQ=255)                                                       aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaac                                                                                                                               >  2:140448/1‑83 (MQ=255)                                                       aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGaa                                                                                                                    <  1:1651481/94‑1 (MQ=255)                                                       aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGaa                                                                                                                    >  2:1651481/1‑94 (MQ=255)                                                       aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtg                                                                       >  2:32130/1‑139 (MQ=255)                                                                    ccACTGACATTATCCCGACAAGAAATCCTCAAAATGATTCGGCAACGTCCTCGTGCGTTAATGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGACGACTTATCCCCCGGTGGTGATGGCACATa                                                          >  2:1093757/1‑139 (MQ=255)                                                                            aTTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                                  >  1:1234236/1‑139 (MQ=255)                                                                                tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGTGTTAACGCCAGaacaa                                                                                                                             >  1:1382669/1‑60 (MQ=255)                                                                                tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGTGTTAACGCCAGaacaa                                                                                                                             <  2:1382669/60‑1 (MQ=255)                                                                                tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                                         <  1:1611091/112‑1 (MQ=255)                                                                                tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                                         >  2:1611091/1‑112 (MQ=255)                                                                                tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAt                                              >  1:1585160/1‑139 (MQ=255)                                                                                  ccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTc                                            <  2:1114949/139‑1 (MQ=255)                                                                                    gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                    <  1:594279/113‑1 (MQ=255)                                                                                    gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                    >  2:594279/1‑113 (MQ=255)                                                                                      cAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTa                                                                                 <  1:634110/98‑1 (MQ=255)                                                                                      cAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTa                                                                                 >  2:634110/1‑98 (MQ=255)                                                                                               cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGcaca                                                            <  1:104154/110‑1 (MQ=255)                                                                                               cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGcaca                                                            >  2:104154/1‑110 (MQ=255)                                                                                               cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATg                               >  2:77959/1‑139 (MQ=255)                                                                                               cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATg                               >  2:456533/1‑139 (MQ=255)                                                                                                             gCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAcc                 <  2:682349/139‑1 (MQ=255)                                                                                                                 cgtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCg                                                                                                    <  2:508551/52‑1 (MQ=255)                                                                                                                 cgtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCg                                                                                                    >  1:508551/1‑52 (MQ=255)                                                                                                                   tcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                    >  2:283822/1‑82 (MQ=255)                                                                                                                   tcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                    <  1:283822/82‑1 (MQ=255)                                                                                                                   tcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCa           >  2:804710/1‑139 (MQ=255)                                                                                                                     gtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACt         <  1:315260/139‑1 (MQ=255)                                                                                                                           cGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTc                                            <  2:238588/98‑1 (MQ=255)                                                                                                                           cGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTc                                            >  1:238588/1‑98 (MQ=255)                                                                                                                            gTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaa  >  1:345213/1‑139 (MQ=255)                                                                                                                                aCGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGa                     >  2:27206/1‑116 (MQ=255)                                                                                                                                aCGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGa                     <  1:27206/116‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  GCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAA  >  NC_000913/3056141‑3056402

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