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

ACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCC  >  NC_000913/3056138‑3056392                                                                                                                                        |                                                                                                                        aCAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaa                                                                                                                      >  1:168321/1‑139 (MQ=255) aCAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaa                                                                                                                      <  2:1390245/139‑1 (MQ=255)    gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                   <  2:9095/139‑1 (MQ=255)    gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                   <  1:1513636/139‑1 (MQ=255)       tGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAAcag                                                                                                                >  1:1433077/1‑139 (MQ=255)         cggcggCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAAcagca                                                                                                              >  1:970449/1‑139 (MQ=255)            ctgcgtctcGGAGATTCCCGTCGTAGCGGGCATCCGCCATGACGCAACTACCAGTACGCCCACTGACGTTATCCCGACAAGAAATCCGCCAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGa                                                                                                           <  2:693593/134‑1 (MQ=255)                           tCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                            <  2:506303/139‑1 (MQ=255)                           tCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                            <  2:168321/139‑1 (MQ=255)                           tCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                            >  1:499316/1‑139 (MQ=255)                                               gCCATGACGCAACTACCAGAACTCTCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACtt                                                                        >  2:1039522/1‑139 (MQ=255)                                               gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACtt                                                                        >  2:875124/1‑139 (MQ=255)                                               gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACtt                                                                        >  1:895009/1‑139 (MQ=255)                                                   tGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATcc                                                                    >  1:803790/1‑139 (MQ=255)                                                    gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaac                                                                                                                     <  2:1231224/89‑1 (MQ=255)                                                    gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaac                                                                                                                     >  1:1231224/1‑89 (MQ=255)                                                    gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                   >  1:214840/1‑139 (MQ=255)                                                    gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                   >  1:605295/1‑139 (MQ=255)                                                        cAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                               <  2:895009/139‑1 (MQ=255)                                                         aaCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                            >  2:1484331/1‑109 (MQ=255)                                                         aaCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                            <  1:1484331/109‑1 (MQ=255)                                                          aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtg                                                             >  2:309456/1‑139 (MQ=255)                                                            tACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtggt                                                           <  2:605295/139‑1 (MQ=255)                                                             aCCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                          <  2:1433077/139‑1 (MQ=255)                                                                 gAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATgg                                                      >  1:1351261/1‑139 (MQ=255)                                                                   aCTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAccc                                                                                    >  2:216097/1‑107 (MQ=255)                                                                   aCTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAccc                                                                                    <  1:216097/107‑1 (MQ=255)                                                                       ccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGAc                                                                          >  1:469737/1‑113 (MQ=255)                                                                       ccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGAc                                                                          <  2:469737/113‑1 (MQ=255)                                                                       ccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATa                                                >  1:1128849/1‑139 (MQ=255)                                                                        cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                     <  2:968453/117‑1 (MQ=255)                                                                        cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                     >  1:968453/1‑117 (MQ=255)                                                                        cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATAc                                               <  2:1128849/139‑1 (MQ=255)                                                                          cTGACTTTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACgg                                             <  2:499316/139‑1 (MQ=255)                                                                            gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATg                                                       <  1:664540/127‑1 (MQ=255)                                                                            gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATg                                                       >  2:664540/1‑127 (MQ=255)                                                                            gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTc                                           >  2:829944/1‑139 (MQ=255)                                                                            gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTc                                           >  1:1381394/1‑139 (MQ=255)                                                                               aTTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                        <  1:829944/139‑1 (MQ=255)                                                                               aTTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                        <  2:214840/139‑1 (MQ=255)                                                                                  aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTa                                     >  2:1338127/1‑139 (MQ=255)                                                                                  aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTa                                     <  2:1381394/139‑1 (MQ=255)                                                                                   tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtgg                                                            >  2:1193516/1‑115 (MQ=255)                                                                                   tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtgg                                                            <  1:1193516/115‑1 (MQ=255)                                                                                   tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAt                                    >  1:820730/1‑139 (MQ=255)                                                                                     ccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggt                                                              <  1:1418198/111‑1 (MQ=255)                                                                                     ccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggt                                                              >  2:1418198/1‑111 (MQ=255)                                                                                       gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                     >  2:1334859/1‑102 (MQ=255)                                                                                       gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                     <  1:1334859/102‑1 (MQ=255)                                                                                       gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCct                                <  2:1428264/139‑1 (MQ=255)                                                                                        aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGcaca                                                  >  1:612274/1‑120 (MQ=255)                                                                                        aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGcaca                                                  <  2:612274/120‑1 (MQ=255)                                                                                        aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACgg                                             >  2:876264/1‑125 (MQ=255)                                                                                        aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACgg                                             <  1:876264/125‑1 (MQ=255)                                                                                         cAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                   <  1:1183962/102‑1 (MQ=255)                                                                                         cAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                   >  2:1183962/1‑102 (MQ=255)                                                                                             aaaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCttt                          >  2:181864/1‑139 (MQ=255)                                                                                                     aaaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCg                  <  1:181864/139‑1 (MQ=255)                                                                                                           aTTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg            <  2:740497/139‑1 (MQ=255)                                                                                                                 cAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAccc      <  1:106181/139‑1 (MQ=255)                                                                                                                     gtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGcc  <  1:309456/139‑1 (MQ=255)                                                                                                                     gtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGcc  <  1:875124/139‑1 (MQ=255)                                                                                                                         tcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTcaaca                                                                                               >  1:1475683/1‑42 (MQ=255)                                                                                                                         tcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTcaaca                                                                                               <  2:1475683/42‑1 (MQ=255)                                                                                                                         tcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACgg                                             >  2:643843/1‑92 (MQ=255)                                                                                                                         tcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACgg                                             <  1:643843/92‑1 (MQ=255)                                                                                                                            tGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatg                                                                               <  1:1411368/55‑1 (MQ=255)                                                                                                                            tGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatg                                                                               >  2:1411368/1‑55 (MQ=255)                                                                                                                               gTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGt                                            >  2:296482/1‑87 (MQ=255)                                                                                                                               gTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGt                                            <  1:296482/87‑1 (MQ=255)                                                                                                                                 tAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                        >  2:253964/1‑89 (MQ=255)                                                                                                                                 tAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                        <  1:253964/89‑1 (MQ=255)                                                                                                                                        |                                                                                                                        ACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCC  >  NC_000913/3056138‑3056392

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