Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,960,127	G→C	5.9%	A22A (GCG→GCC)	argS →	arginyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,960,127	0	G	C	5.9%	244.2 / 5.4	68	A22A (GCG→GCC)	argS	arginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base C (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TTAACTGGTGCGGGCAATTTTTGCTCGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACAATTAGCAGAGCAGGTGCTGACTCATCTGG  >  NC_000913/1960012‑1960260                                                                                                                    |                                                                                                                                      ttAACTGGTGCGGGCAATTTTTGCTCGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCg                                                                                                                                <  1:186475/123‑1 (MQ=255) ttAACTGGTGCGGGCAATTTTTGCTCGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCg                                                                                                                                >  2:186475/1‑123 (MQ=255)                     ttGCTCGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCg                                                                                                                   >  4:1088009/1‑116 (MQ=255)                     ttGCTCGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCg                                                                                                                   <  3:1088009/116‑1 (MQ=255)                          cGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGc                                                                               <  4:431612/147‑1 (MQ=255)                          cGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGc                                                                               >  3:431612/1‑147 (MQ=255)                            cTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCgt                                                                                                                  >  2:77421/1‑110 (MQ=255)                            cTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCgt                                                                                                                  <  1:77421/110‑1 (MQ=255)                                  caATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCat                                                                      >  2:891506/1‑148 (MQ=255)                                     tGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCaaa                                                                                                       >  2:584235/1‑112 (MQ=255)                                     tGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCaaa                                                                                                       <  1:584235/112‑1 (MQ=255)                                      gTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCatgatg                                                                  >  4:291346/1‑148 (MQ=255)                                          ggTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGtt                                                                                                                     >  2:447215/1‑93 (MQ=255)                                          ggTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGtt                                                                                                                     <  1:447215/93‑1 (MQ=255)                                          ggTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAg                                                              >  1:934579/1‑148 (MQ=255)                                           gTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGtt                                                                                                                     >  1:666759/1‑92 (MQ=255)                                           gTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGtt                                                                                                                     <  2:666759/92‑1 (MQ=255)                                           gTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGt                                                             <  3:291346/148‑1 (MQ=255)                                            tATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCa                                                               <  3:1105681/145‑1 (MQ=255)                                            tATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCa                                                               >  4:1105681/1‑145 (MQ=255)                                              ttCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGc                                                          <  2:934579/148‑1 (MQ=255)                                                  ggTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCatgatg                                                                  >  1:1095021/1‑136 (MQ=255)                                                  ggTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCatgatg                                                                  <  2:1095021/136‑1 (MQ=255)                                                     gAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCgg                                                                                           <  4:840177/108‑1 (MQ=255)                                                     gAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCgg                                                                                           >  3:840177/1‑108 (MQ=255)                                                     gAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATgg                                                                 >  4:1072135/1‑134 (MQ=255)                                                     gAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATgg                                                                 <  3:1072135/134‑1 (MQ=255)                                                       atatTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACt                                                 <  1:891506/148‑1 (MQ=255)                                                           tCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAagttca                                                                                                 <  3:542836/96‑1 (MQ=255)                                                           tCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAagttca                                                                                                 >  4:542836/1‑96 (MQ=255)                                                           tCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGt                                             >  3:516411/1‑148 (MQ=255)                                                            cAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTaaaaaa                                                   >  4:148273/1‑141 (MQ=255)                                                            cAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTaaaaaa                                                   <  3:148273/141‑1 (MQ=255)                                                             aGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGa                                                                                        >  1:374978/1‑103 (MQ=255)                                                             aGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGa                                                                                        <  2:374978/103‑1 (MQ=255)                                                             aGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTAt                                           <  4:516411/148‑1 (MQ=255)                                                                 tcttctCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAagttc                                                                                                  <  4:571477/89‑1 (MQ=255)                                                                 tcttctCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAagttc                                                                                                  >  3:571477/1‑89 (MQ=255)                                                                   ttctCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTaaaaa                                                    >  2:1072069/1‑133 (MQ=255)                                                                   ttctCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTaaaaa                                                    <  1:1072069/133‑1 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGAtt                                                                                                                                   <  4:671727/48‑1 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGAtt                                                                                                                                   >  3:671727/1‑48 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCatgatg                                                                  <  3:814991/113‑1 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCatgatg                                                                  >  4:814991/1‑113 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGc                                                                >  3:226279/1‑115 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGc                                                                <  4:226279/115‑1 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACa                               >  3:759103/1‑148 (MQ=255)                                                                         cAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACa                               >  1:734422/1‑148 (MQ=255)                                                                           gAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACAAt                             <  4:759103/148‑1 (MQ=255)                                                                             aaaaTTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGGTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACAATTa                           >  2:609593/1‑148 (MQ=255)                                                                                  tCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTg                                                >  2:559175/1‑122 (MQ=255)                                                                                  tCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTg                                                <  1:559175/122‑1 (MQ=255)                                                                                  tCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAAGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGt                                                             >  2:3711/1‑109 (MQ=255)                                                                                  tCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAAGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGt                                                             <  1:3711/109‑1 (MQ=255)                                                                                       cAGGCCATGATTGCTGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACAATTAGCAGAGCAgg                 <  1:609593/148‑1 (MQ=255)                                                                                           ccATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTgg                                               >  4:1110650/1‑114 (MQ=255)                                                                                           ccATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTgg                                               <  3:1110650/114‑1 (MQ=255)                                                                                              tGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACa                               <  1:838081/127‑1 (MQ=255)                                                                                              tGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACa                               >  2:838081/1‑127 (MQ=255)                                                                                                    ccgCAGGCGCGCCTGCCGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTa                                                                             >  4:994097/3‑75 (MQ=255)                                                                                                    ccgCAGGCGCGCCTGCCGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTa                                                                             >  4:210959/3‑75 (MQ=255)                                                                                                    ccgCAGGCGCGCCTGCCGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTa                                                                             <  3:210959/73‑1 (MQ=255)                                                                                                    ccgCAGGCGCGCCTGCCGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTa                                                                             <  3:994097/73‑1 (MQ=255)                                                                                                    cGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTg                                                <  2:668935/104‑1 (MQ=255)                                                                                                    cGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTg                                                >  1:668935/1‑104 (MQ=255)                                                                                                      gCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGa                                 <  4:308821/117‑1 (MQ=255)                                                                                                      gCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGa                                 >  3:308821/1‑117 (MQ=255)                                                                                                      gCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACAATTAGCAGAGCAGGTGCTGACTCATCTgg  <  2:734422/148‑1 (MQ=255)                                                                                                                    |                                                                                                                                      TTAACTGGTGCGGGCAATTTTTGCTCGCTTCATCAATGTAAGGTATTCCGGTGAATATTCAGGCTCTTCTCTCAGAAAAAGTCCGTCAGGCCATGATTGCGGCAGGCGCGCCTGCGGATTGCGAACCGCAGGTTCGTCAGTCAGCAAAAGTTCAGTTCGGCGACTATCAGGCTAACGGCATGATGGCAGTTGCTAAAAAACTGGGTATGGCACCGCGACAATTAGCAGAGCAGGTGCTGACTCATCTGG  >  NC_000913/1960012‑1960260

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