Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	961,344	A→G	L48L (TTA→TTG)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,344	0	A	G	100.0%	227.1 / NA	63	L48L (TTA→TTG)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/31);  total (32/31)

TTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACT  >  NC_000913/961211‑961479                                                                                                                                      |                                                                                                                                        ttGCCCCGGTTATTACCATTGATGGCCCAAGGGGTGCAGGGAAAGGCGCCTTGGGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGATGCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATcgc                                                                                                                                   <  2:384220/140‑3 (MQ=255)   gCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCa                                                                                                                                 <  2:337081/140‑1 (MQ=255)         gTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCatgttga                                                                                                                           >  2:1027979/1‑140 (MQ=255)           tattaCCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCatgttgatg                                                                                                                         <  2:716411/140‑1 (MQ=255)             ttaCCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGTATTGCATCACCatgttgatgtt                                                                                                                       <  1:362291/140‑1 (MQ=255)             ttaCCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCatgttgatgtt                                                                                                                       <  1:367753/140‑1 (MQ=255)              taCCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCatgttgatgttg                                                                                                                      >  1:722130/1‑140 (MQ=255)                   ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCAc                                                                                                                                   >  2:564495/1‑122 (MQ=255)                   ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCAc                                                                                                                                   <  1:564495/122‑1 (MQ=255)                   ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCg                                                                                                                 >  2:419378/1‑140 (MQ=255)                   ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCg                                                                                                                 >  2:1006288/1‑140 (MQ=255)                        ggCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGa                                                                                                            >  1:216160/1‑140 (MQ=255)                        ggCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGa                                                                                                            >  1:306725/1‑140 (MQ=255)                          cccAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATg                                                                                                          <  1:153971/140‑1 (MQ=255)                             aaGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATgcgc                                                                                                       <  1:1006288/140‑1 (MQ=255)                                   tGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTAc                                                                                                 <  1:419378/140‑1 (MQ=255)                                       gggAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCt                                                                                             <  1:625464/140‑1 (MQ=255)                                        ggAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTg                                                                                            <  2:216160/140‑1 (MQ=255)                                           aaGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                         <  2:488011/140‑1 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                         >  1:1004908/1‑138 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                         <  2:95258/138‑1 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                         <  2:1004908/138‑1 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                         >  1:95258/1‑138 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                                                       >  1:896375/1‑140 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                                                       >  2:234377/1‑140 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                                                       >  1:994011/1‑140 (MQ=255)                                             ggCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGATGATGTTGCGTCGGAACATGCGCTGGTACCGCTGGCATc                                                                                       >  1:753676/1‑140 (MQ=255)                                                 ccTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCAt                                                                                   <  2:10306/140‑1 (MQ=255)                                                  cTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATc                                                                                  >  2:381609/1‑140 (MQ=255)                                                  cTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATc                                                                                  >  1:774169/1‑140 (MQ=255)                                                    tgtgTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTg                                                                                >  1:35535/1‑140 (MQ=255)                                                                tGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGtttt                                                                    >  1:308029/1‑140 (MQ=255)                                                                tGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGtttt                                                                    >  1:567344/1‑140 (MQ=255)                                                                cggCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGtt                                                                      >  1:747313/2‑138 (MQ=255)                                                                cggCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGtt                                                                      <  2:747313/137‑1 (MQ=255)                                                                    ggAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTgtgt                                                                <  1:510124/140‑1 (MQ=255)                                                                      aaGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTgtgt                                                                >  2:77891/1‑138 (MQ=255)                                                                      aaGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTgtgt                                                                <  1:77891/138‑1 (MQ=255)                                                                       aGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGCACGTTTTGTGTCGa                                                             >  2:585195/1‑140 (MQ=255)                                                                           ttGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCaa                                                         <  2:812417/140‑1 (MQ=255)                                                                            tGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAAt                                                        >  1:853885/1‑140 (MQ=255)                                                                            tGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGCACGTTTTGTGTCGACCAAt                                                        >  1:917305/1‑140 (MQ=255)                                                                                aTGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCa                                                    <  1:740316/140‑1 (MQ=255)                                                                                   gCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATc                                                 <  2:994011/140‑1 (MQ=255)                                                                                     aTCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTg                                               >  1:949803/1‑140 (MQ=255)                                                                                      tCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCTCCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTgg                                              <  2:910218/140‑1 (MQ=255)                                                                                        tgctgGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGaa                                            >  2:1041647/1‑140 (MQ=255)                                                                                         gctgGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAg                                           <  2:308029/140‑1 (MQ=255)                                                                                         gctgGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAg                                           <  2:35535/140‑1 (MQ=255)                                                                                         gctgGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAg                                           >  2:701653/1‑140 (MQ=255)                                                                                                tCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCt                                    <  1:69241/140‑1 (MQ=255)                                                                                                tCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCACTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCt                                    <  2:722130/140‑1 (MQ=255)                                                                                                             aTCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCt                                    >  1:829656/1‑127 (MQ=255)                                                                                                             aTCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCt                                    <  2:829656/127‑1 (MQ=255)                                                                                                             aTCGCGTACTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGAt                       <  2:774169/140‑1 (MQ=255)                                                                                                                    aCTGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAgcg                >  2:400318/1‑140 (MQ=255)                                                                                                                      tGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAgcggc              >  2:1026597/1‑140 (MQ=255)                                                                                                                      tGGCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAgcggc              <  1:1027979/140‑1 (MQ=255)                                                                                                                        gCATTGGCGGCATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGa            >  1:230530/1‑140 (MQ=255)                                                                                                                             ggcggcATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTc       <  1:701653/140‑1 (MQ=255)                                                                                                                             ggcggcATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTc       >  1:221714/1‑140 (MQ=255)                                                                                                                               cggcATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGt     <  1:381609/140‑1 (MQ=255)                                                                                                                                  cATTGCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACt  >  1:229497/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                        TTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACT  >  NC_000913/961211‑961479

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