Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	646,500	G→A	P27P (CCC→CCT)	citT ←	citrate/succinate antiporter, citrate carrier

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,500	0	G	A	100.0%	236.7 / NA	73	P27P (CCC→CCT)	citT	citrate/succinate antiporter, citrate carrier
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/37);  total (36/37)

GTAATTACTGCCAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCGGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAAAT  >  NC_000913/646365‑646628                                                                                                                                        |                                                                                                                                 gTAATTACTGCCAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCaggg                                                                                                                               >  2:1324330/1‑139 (MQ=255)  tAATTACTGCCAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGa                                                                                                                              <  1:348944/139‑1 (MQ=255)   aaTTACTGCCAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTACGGCGGCATACCGTCAGGGAc                                                                                                                             >  1:455969/1‑139 (MQ=255)     ttACTGCCAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAg                                                                                                                           <  1:394805/139‑1 (MQ=255)           ccAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGAc                                                                                                                             <  2:506420/131‑1 (MQ=255)           ccAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGAc                                                                                                                             >  1:506420/1‑131 (MQ=255)           ccAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGAt                                                                                                                       >  2:1508425/1‑137 (MQ=255)           ccAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGAt                                                                                                                       <  1:1508425/137‑1 (MQ=255)           ccAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATaa                                                                                                                     >  2:1469615/1‑139 (MQ=255)           ccAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATaa                                                                                                                     >  1:49222/1‑139 (MQ=255)             aaTAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGa                                                                                                                   <  2:455969/139‑1 (MQ=255)               tAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAg                                                                                                                           >  1:716794/1‑129 (MQ=255)               tAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAg                                                                                                                           <  2:716794/129‑1 (MQ=255)                    cAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGa                                                                                                            <  2:49222/139‑1 (MQ=255)                       aTAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGacac                                                                                                         <  2:302/139‑1 (MQ=255)                       aTAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGacac                                                                                                         >  1:302/1‑139 (MQ=255)                             aCCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGAc                                                                                                                             >  1:880360/1‑113 (MQ=255)                             aCCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGAc                                                                                                                             <  2:880360/113‑1 (MQ=255)                             aCCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCG‑AGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATcac                                                                                                  >  1:499991/1‑139 (MQ=255)                              ccGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATcac                                                                                                  <  1:1788056/139‑1 (MQ=255)                                    tAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAg                                                                                            >  1:1192995/1‑139 (MQ=255)                                         cTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTgggg                                                                                       >  1:615358/1‑139 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACAccc                                                                                                       <  2:1038780/120‑1 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACAccc                                                                                                       >  1:1038780/1‑120 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATcac                                                                                                  <  1:1591198/125‑1 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATcac                                                                                                  >  2:1591198/1‑125 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCa                                                                                    >  2:755246/1‑139 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGAAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATg                                                                                                             <  2:1438571/114‑1 (MQ=255)                                            aTCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGAAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATg                                                                                                             >  1:1438571/1‑114 (MQ=255)                                              cGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGa                                                                                                                              <  2:1481997/95‑1 (MQ=255)                                              cGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGa                                                                                                                              >  1:1481997/1‑95 (MQ=255)                                               gCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCG‑AGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAg                                                                                <  2:499991/139‑1 (MQ=255)                                                 tGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGTCATGCCTGCGGCGGCATACCGTCAGGGa                                                                                                                              >  2:1106053/1‑92 (MQ=255)                                                 tGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTCCCACAAACAAAGCGAAGTAATGTCATGCCTGCGGCGGCATACCGTCAGGGa                                                                                                                              <  1:1106053/92‑1 (MQ=255)                                                           aTTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCatatat                                                                     >  2:32236/1‑139 (MQ=255)                                                                cTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCt                                                                >  2:647539/1‑139 (MQ=255)                                                                 tCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCtt                                                               <  2:1192995/139‑1 (MQ=255)                                                                      gATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCt                                                          <  1:647539/139‑1 (MQ=255)                                                                           tGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCaa                                                                                   >  2:262681/1‑109 (MQ=255)                                                                           tGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCaa                                                                                   <  1:262681/109‑1 (MQ=255)                                                                           tGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGa                                                     <  2:615358/139‑1 (MQ=255)                                                                           tGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGa                                                     <  1:750068/139‑1 (MQ=255)                                                                               gACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCAt                                                                                                     >  2:1518302/1‑87 (MQ=255)                                                                               gACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCAt                                                                                                     <  1:1518302/87‑1 (MQ=255)                                                                               gACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATa                                                 <  1:1324330/139‑1 (MQ=255)                                                                                       ttGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGaca                                                                                                          <  2:1014830/74‑1 (MQ=255)                                                                                       ttGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGaca                                                                                                          >  1:1014830/1‑74 (MQ=255)                                                                                          ccACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCg                                      <  1:1469615/139‑1 (MQ=255)                                                                                            acaAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGaa                                    >  2:584805/1‑139 (MQ=255)                                                                                                    aGCGAAGTAGTGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAata                            >  2:1076242/1‑139 (MQ=255)                                                                                                    aGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAata                            >  2:1283788/1‑139 (MQ=255)                                                                                                          gTAATGCCATGCCTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCATCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCa                      <  1:1283788/139‑1 (MQ=255)                                                                                                              tgccatgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGacac                                                                                                         <  1:574304/52‑1 (MQ=255)                                                                                                              tgccatgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGacac                                                                                                         >  2:574304/1‑52 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACAccc                                                                                                       >  1:206035/1‑49 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACAccc                                                                                                       <  2:206035/49‑1 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAg                                                                                >  2:143879/1‑72 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAg                                                                                <  1:143879/72‑1 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTaaa                                                       <  2:1292887/97‑1 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTaaa                                                       >  1:1292887/1‑97 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCa             >  1:1220069/1‑139 (MQ=255)                                                                                                                   tgccTGCGGCGGCATACCGTCAGGGACA‑GGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCaa            >  2:1496883/1‑139 (MQ=255)                                                                                                                     ccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTcc                                                                           >  1:1023755/1‑75 (MQ=255)                                                                                                                     ccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTcc                                                                           <  2:1240781/75‑1 (MQ=255)                                                                                                                     ccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTcc                                                                           <  2:1023755/75‑1 (MQ=255)                                                                                                                     ccTGCGGCGGCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTcc                                                                           >  1:1240781/1‑75 (MQ=255)                                                                                                                         cggcggCATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAATTTAATATTTCCAAATTTATCAAGTGCt       >  1:501468/1‑139 (MQ=255)                                                                                                                            cggcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTaa    <  2:501468/139‑1 (MQ=255)                                                                                                                            cggcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTaa    <  1:32236/139‑1 (MQ=255)                                                                                                                            cggcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTaa    <  1:1121095/139‑1 (MQ=255)                                                                                                                            cggcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTaa    <  1:1076242/139‑1 (MQ=255)                                                                                                                              gcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATc              <  1:794041/127‑1 (MQ=255)                                                                                                                              gcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATc              >  2:794041/1‑127 (MQ=255)                                                                                                                              gcATACCGTCAGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAaat  >  1:179008/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                 GTAATTACTGCCAATAACGCAAATAGTAACCGCAATAAAACTGATCGCTGTTGCCGGAATTGGCTCGAGGATCATGCCGACAATCATTGCCACAAACACAGCGAAGTAATGCCATGCCTGCGGCGGCATACCGTCGGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAAAT  >  NC_000913/646365‑646628

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