Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	44.6%	95.1 / 16.3	57	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (28/3);  new base C (0/25);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.18e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAT  >  NC_000913/4537960‑4538255                                                                                                                                                        |                                                                                                                                                   aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCAc                                                                                                                                                   >  2:317616/1‑151 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCAcc                                                                                                                                                  <  1:923108/150‑1 (MQ=255)     ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                               >  1:868568/1‑151 (MQ=255)              cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                                      <  2:324893/151‑1 (MQ=255)                  aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgcc                                                                                                                                  >  1:469937/1‑151 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                                >  1:436219/1‑150 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                             >  2:850403/1‑150 (MQ=255)                                       gaCGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTAAAATCAATTTCAAATACGGCACCCGTCCGCAATCCTGTTTT‑GGCTTCG‑CCATTAATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCa                                                                                                              <  2:403161/150‑1 (MQ=255)                                       gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                               >  1:277904/1‑149 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                    >  2:962197/1‑151 (MQ=255)                                                    ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                  >  2:210278/1‑149 (MQ=255)                                                         ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                            >  2:380409/1‑150 (MQ=255)                                                                aatttattcCCGGTAAAATCAATTTAAAATCCGGCACCCGTTCGAATCCCTGTTTT‑GGTTTCC‑CCATTAATAAGCCAGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                                      <  2:830050/146‑1 (MQ=255)                                                                aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                    >  1:233832/1‑151 (MQ=255)                                                                 aaTTATTCCCGGTAAAATAAATTTAAAATACGCCATCCGTTC‑CAAACCCGGTTTGGGCTCCG‑CAATTAAAAAGCCAGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                   <  2:634387/151‑1 (MQ=255)                                                                     attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                               >  2:901099/1‑151 (MQ=255)                                                                      ttcccGGTAAAATCAAGTTCAAATCCGGCATCCGTTCGCAACCCTGTTTG‑GGCTTCC‑CCATTAAAAAGCCAGGTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                               <  2:441449/147‑1 (MQ=255)                                                                      ttaCGGGTGAAATCAAGTACAAATACGGCATCCGTCCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                               <  2:41672/150‑1 (MQ=255)                                                                       taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                             >  2:682280/1‑151 (MQ=255)                                                                       taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                             >  2:826676/1‑151 (MQ=255)                                                                         cgggTAAAATCAATTTAAAATACGGCATCCGTTCGCAATCCTGTTTG‑GGCTTCG‑CCATAAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                           <  2:485237/149‑1 (MQ=255)                                                                         ccGGTAAAATAAATTTAAAATCCGGCATCCGTTCGAAATCCTGGTTT‑GGTTTCC‑CAATTAATAAGCCAGTTTTTATCCCCTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                            <  2:432405/150‑1 (MQ=255)                                                                               aaaTCAAGTAAAAATCCGTCATCCTTTCGAAATCCTGGTTT‑GGCTCCG‑CCATTAAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                     <  2:771258/151‑1 (MQ=255)                                                                                 aTCAAGTTAAAATACGGCATCCTTTCGCAATCCTGGTTT‑GGCTTCA‑CAATTAATAACCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                   <  2:24989/151‑1 (MQ=255)                                                                                    aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                  >  2:546083/1‑149 (MQ=255)                                                                                    aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                  >  2:509669/1‑149 (MQ=255)                                                                                          aaaTCCGGCTTCCGTTCCCAACCTTGTTTT‑GGCTTCG‑CAATAAATAAGCAAGTTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                           <  2:868584/150‑1 (MQ=255)                                                                                          aaaTACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAAAAAGCCAGGTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                             <  2:348734/148‑1 (MQ=255)                                                                                                 gCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                   >  2:819656/1‑151 (MQ=255)                                                                                                     cctttCGCAACCTTGGTTT‑GGTTTCG‑CTTTAAATAACCCAGGTTTTACCCCTTTAATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                               <  1:853044/148‑1 (MQ=255)                                                                                                     ccGTTCGCAATCCTGTTTT‑GGCTTCG‑CCATTAATAACCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                               <  2:98113/151‑1 (MQ=255)                                                                                                     ccGTTCGCAATCCTGGTTT‑GGTTTCG‑CCATTAAAAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                <  2:254680/150‑1 (MQ=255)                                                                                                     ccGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                               >  1:299650/1‑151 (MQ=255)                                                                                                      cGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                              >  1:894031/1‑151 (MQ=255)                                                                                                       gTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                              >  2:931306/1‑150 (MQ=255)                                                                                                        ttCGCAATCCTGTTTT‑GGCTTCG‑CCATTAATAAGCCGGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                              <  2:437013/149‑1 (MQ=255)                                                                                                        ttCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                             >  2:791136/1‑150 (MQ=255)                                                                                                                ccTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                    >  2:89523/1‑151 (MQ=255)                                                                                                                 cgtgTTT‑GGTTTCCCCCATTAATAGCCCGGTTTTTACCGCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                              <  1:206908/138‑1 (MQ=255)                                                                                                                    gTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                >  2:390832/1‑151 (MQ=255)                                                                                                                     ttt‑GGCTTCG‑CAATTAAGAAGCCGGTTTTTATCCCCTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                    <  1:112826/146‑1 (MQ=255)                                                                                                                      tt‑GGCTCCC‑CCTTAAAAAAGCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                               <  2:889957/150‑1 (MQ=255)                                                                                                                          gCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                           >  1:596500/1‑151 (MQ=255)                                                                                                                             tCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                        >  2:94419/1‑151 (MQ=255)                                                                                                                              cG‑CAATTAAAAACCCAGTTTTAACCACTTTAATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                          <  1:967315/148‑1 (MQ=255)                                                                                                                                 ccATTAATAAGCCAGTTTTTATCCCCTTTACTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                          <  2:467114/146‑1 (MQ=255)                                                                                                                                 caaTTAATAAGCCAGGTTTTACCCCCTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                         <  2:113195/145‑1 (MQ=255)                                                                                                                                   aTTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                    >  1:716159/1‑150 (MQ=255)                                                                                                                                      aataaGCCAGTTTTTATCCCCTTTTTTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                <  2:382129/151‑1 (MQ=255)                                                                                                                                      aaaaacccAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                  <  1:572956/143‑1 (MQ=255)                                                                                                                                      aaaaaccagctttttATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCTAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                   <  1:913910/137‑1 (MQ=255)                                                                                                                                      aaaaaGCCAGGTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                 <  2:617463/147‑1 (MQ=255)                                                                                                                                              aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc        >  1:360025/1‑151 (MQ=255)                                                                                                                                               gtttttATCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc        <  1:74754/148‑1 (MQ=255)                                                                                                                                                 ttttACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa       <  2:510030/149‑1 (MQ=255)                                                                                                                                                   ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg    >  2:866767/1‑150 (MQ=255)                                                                                                                                                    tATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt  >  2:743523/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                   ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCG‑CCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAT  >  NC_000913/4537960‑4538255

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