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%	160.0 / NA	51	P11P (CCG→CCA)	fau	5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/29);  total (22/29)

AGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGA  >  NC_000913/3056140‑3056401                                                                                                                                      |                                                                                                                                 aGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaaca                                                                                                                             <  2:470073/139‑1 (MQ=255)          cggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGa                                                                                                                    >  1:1100547/1‑139 (MQ=255)               tctcGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATgg                                                                                                               <  1:1297706/139‑1 (MQ=255)                          cccTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGcc                                                                                                    >  2:343046/1‑139 (MQ=255)                                ttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                              <  2:209835/139‑1 (MQ=255)                                       gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaaca                                                                                                                             <  1:233589/101‑1 (MQ=255)                                       gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaaca                                                                                                                             >  2:233589/1‑101 (MQ=255)                                       gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGCAATGGGTCAACAAGCCGCTACCCGgatga                                                                                       >  2:316004/1‑139 (MQ=255)                                       gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatga                                                                                       >  2:700914/1‑139 (MQ=255)                                       gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatga                                                                                       <  1:737521/139‑1 (MQ=255)                                       gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatga                                                                                       >  2:293487/1‑139 (MQ=255)                                        cACCAGCCATGACGCAACTACCAGAACTCCCACTGATATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatgat                                                                                      <  2:401238/139‑1 (MQ=255)                                        cACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatgat                                                                                      <  1:700914/139‑1 (MQ=255)                                             gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACtt                                                                                 >  1:288975/1‑139 (MQ=255)                                             gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACtt                                                                                 >  1:1223033/1‑139 (MQ=255)                                                  gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                            <  2:1100547/139‑1 (MQ=255)                                                     gCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacga                                                                                                                            <  2:777758/88‑3 (MQ=255)                                                     gCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                            >  1:777758/1‑88 (MQ=255)                                                      cAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                                        <  2:1223033/139‑1 (MQ=255)                                                      cAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                                        <  1:343046/139‑1 (MQ=255)                                                       aaCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggt                                                                       >  1:30960/1‑139 (MQ=255)                                                        aCTACCAGAACTCCCACTGACATTATCCCTACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACTGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtg                                                                      <  1:875727/139‑1 (MQ=255)                                                        aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGTATGACGACTTATCCTCCggtg                                                                      >  2:623732/1‑139 (MQ=255)                                                           aCCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                   <  1:293487/139‑1 (MQ=255)                                                           aCCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                   <  1:316004/139‑1 (MQ=255)                                                           aCCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                   <  2:1217229/139‑1 (MQ=255)                                                            ccAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGa                                                                  <  2:288975/139‑1 (MQ=255)                                                                   tCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGcaca                                                           >  2:743342/1‑139 (MQ=255)                                                                     ccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGga                                                                                          <  1:746372/106‑1 (MQ=255)                                                                     ccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGga                                                                                          >  2:746372/1‑106 (MQ=255)                                                                      cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATgg                                                                                                               <  1:942912/84‑1 (MQ=255)                                                                      cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATgg                                                                                                               >  2:942912/1‑84 (MQ=255)                                                                      cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATAc                                                        <  1:489490/139‑1 (MQ=255)                                                                             aTTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                                 <  1:702248/139‑1 (MQ=255)                                                                                aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTa                                              <  2:81395/139‑1 (MQ=255)                                                                                 tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAt                                             >  1:175202/1‑139 (MQ=255)                                                                                   ccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGCACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGTTGGCACATACGGTCGCTGTATTc                                           <  1:623732/139‑1 (MQ=255)                                                                                    cGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTcc                                          <  1:208250/139‑1 (MQ=255)                                                                                     gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCct                                         >  2:1247039/1‑139 (MQ=255)                                                                                      aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatg                                                                                        >  2:213943/1‑91 (MQ=255)                                                                                      aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatg                                                                                        <  1:213943/91‑1 (MQ=255)                                                                                      aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctc                                        >  2:957391/1‑139 (MQ=255)                                                                                                cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATg                              <  1:1020227/139‑1 (MQ=255)                                                                                                               cAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAccc               >  2:1023003/1‑139 (MQ=255)                                                                                                                  cgtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc            <  1:957391/139‑1 (MQ=255)                                                                                                                    tcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCa          <  2:30960/139‑1 (MQ=255)                                                                                                                       tcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc            >  2:317148/1‑134 (MQ=255)                                                                                                                       tcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc            <  1:317148/134‑1 (MQ=255)                                                                                                                            cGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                              >  2:1051808/1‑47 (MQ=255)                                                                                                                            cGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                              <  1:1051808/47‑1 (MQ=255)                                                                                                                            cGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGa  >  1:609357/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                 AGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGA  >  NC_000913/3056140‑3056401

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