Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,710	C→A	S14R (AGC→AGA)	dusA →	tRNA‑dihydrouridine synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,710	0	C	A	100.0%	220.6 / NA	63	S14R (AGC→AGA)	dusA	tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/31);  total (32/31)

GTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261573‑4261844                                                                                                                                          |                                                                                                                                       gttATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAgagc                                                                                                                                      >  2:510079/1‑140 (MQ=255)            gTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaa                                                                                                                           >  1:198469/1‑140 (MQ=255)               ccccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaaaaa                                                                                                                        >  1:55970/1‑140 (MQ=255)               ccccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGaaaaaa                                                                                                                        >  1:119549/1‑140 (MQ=255)                 ccGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACt                                                                                                                      >  2:105287/1‑140 (MQ=255)                  cGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTg                                                                                                                     >  1:84380/1‑140 (MQ=255)                    cTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGAc                                                                                                                   <  1:258777/140‑1 (MQ=255)                        tGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTc                                                                                                               <  1:177050/140‑1 (MQ=255)                                   aaTAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGt                                                                                                    <  2:397210/140‑1 (MQ=255)                                   aaTAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGt                                                                                                    <  1:84566/140‑1 (MQ=255)                                        tGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGtt                                                                                               >  2:182713/1‑140 (MQ=255)                                        tGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGtt                                                                                               >  1:219391/1‑140 (MQ=255)                                         gATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGttt                                                                                              >  2:260540/1‑140 (MQ=255)                                          aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  1:189580/1‑140 (MQ=255)                                          aTGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTa                                                                                             >  2:560059/1‑140 (MQ=255)                                             cTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCg                                                                                          >  2:601852/1‑140 (MQ=255)                                                     cataacCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCa                                                                                  >  1:201450/1‑140 (MQ=255)                                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:541978/1‑140 (MQ=255)                                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  1:456658/1‑140 (MQ=255)                                                          ccccATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACAGGAGTGGTCGGTTTAGCGTTGCACCAATGCt                                                                             >  2:241739/1‑140 (MQ=255)                                                            ccATGATATATCGATAGAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGAGGTCGGTTTAGCGTTGCACCAATGCTCg                                                                           >  2:530977/1‑140 (MQ=255)                                                               tGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACt                                                                        >  1:131056/1‑140 (MQ=255)                                                                gATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTg                                                                       <  1:365620/140‑1 (MQ=255)                                                                 atatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgg                                                                      <  2:198469/140‑1 (MQ=255)                                                                 atatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgg                                                                      <  1:531630/140‑1 (MQ=255)                                                                  tatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgga                                                                     <  1:305376/140‑1 (MQ=255)                                                                  tatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgga                                                                     <  2:131056/140‑1 (MQ=255)                                                                  tatatCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTgga                                                                     <  1:232984/140‑1 (MQ=255)                                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  >  2:375665/1‑140 (MQ=255)                                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  >  2:303242/1‑140 (MQ=255)                                                                     atCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTggacgg                                                                  >  1:232119/1‑140 (MQ=255)                                                                         aTAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacag                                                              >  1:287783/1‑140 (MQ=255)                                                                          tAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacaga                                                             >  2:48751/1‑140 (MQ=255)                                                                          tAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacaga                                                             >  1:404358/1‑140 (MQ=255)                                                                            aaaTAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAGCCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGgacagaca                                                           >  2:156996/1‑140 (MQ=255)                                                                                       aCATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTAtt                                                >  1:336752/1‑140 (MQ=255)                                                                                           ttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTa                                                  >  2:377581/1‑134 (MQ=255)                                                                                           ttGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTa                                                  <  1:377581/134‑1 (MQ=255)                                                                                                  tGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTg                                     >  1:12718/1‑140 (MQ=255)                                                                                                   gCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGc                                    >  2:591411/1‑140 (MQ=255)                                                                                                    cACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCt                                   <  2:119549/140‑1 (MQ=255)                                                                                                     aCGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGTCGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCtt                                  <  2:212119/140‑1 (MQ=255)                                                                                                     aCGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCtt                                  <  2:232119/140‑1 (MQ=255)                                                                                                          aaTTCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCg                             <  2:336752/140‑1 (MQ=255)                                                                                                            ttCTGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCa                           <  2:149876/140‑1 (MQ=255)                                                                                                               tGAAATGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATa                        <  2:533031/140‑1 (MQ=255)                                                                                                                 aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      >  2:65222/1‑140 (MQ=255)                                                                                                                 aaaTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACg                      <  2:320640/140‑1 (MQ=255)                                                                                                                   aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    <  1:600498/140‑1 (MQ=255)                                                                                                                   aTGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACAGGAGGGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGtt                    >  2:114522/1‑140 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  1:182713/140‑1 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:339668/140‑1 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  1:185982/140‑1 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:591539/140‑1 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  1:241739/140‑1 (MQ=255)                                                                                                                    tGCAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTg                   <  2:561066/140‑1 (MQ=255)                                                                                                                      cAAAAAATCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCt                 >  2:382750/1‑140 (MQ=255)                                                                                                                         aaaaTCAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGta              <  2:12718/140‑1 (MQ=255)                                                                                                                              cAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCg         >  2:544340/1‑140 (MQ=255)                                                                                                                              cAACCAAACCAGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCg         <  2:583527/140‑1 (MQ=255)                                                                                                                                accaaaccaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGaa       <  1:601852/140‑1 (MQ=255)                                                                                                                                     accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  1:331262/140‑1 (MQ=255)                                                                                                                                     accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  1:419803/140‑1 (MQ=255)                                                                                                                                     accaGAGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGt  <  2:189580/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       GTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGCACGGTAATTCTGAAATGCAAAAAATCAACCAAACCAGCGCAATGCCTGAAAAAACTGACGTTCACTGGAGTGGTCGGTTTAGCGTTGCACCAATGCTCGACTGGACGGACAGACATTGCCGCTATTTCTTGCGTCTGCTTTCCCGCAATACGTTGCTGTATACCGAAATGGT  >  NC_000913/4261573‑4261844

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