Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,692	0	C	A	44.8%	73.3 / 48.8	58	intergenic (+318/+1055)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (5/27);  new base A (26/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 8.16e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTA  >  NC_000913/2903548‑2903838                                                                                                                                                 |                                                                                                                                                      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAacaccc                                                                                                                                                  >  1:163915/1‑148 (MQ=255) aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTc                                                                                                                                                  <  1:10962/150‑1 (MQ=255)   aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCa                                                                                                                                                 >  2:745861/1‑149 (MQ=255)           ttAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTATcc                                                                                                                                           >  2:439930/1‑147 (MQ=255)                tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTAATCT‑Cacaaaa                                                                                                                                  >  2:326136/1‑148 (MQ=255)                tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAACAACACCCAATAATCT‑Ccaaa                                                                                                                                    >  2:293133/1‑148 (MQ=255)                 acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACa                                                                                                                                  <  1:592819/149‑1 (MQ=255)                      cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGAc                                                                                                                             <  1:436200/149‑1 (MQ=255)                      cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACa                                                                                                                            <  2:727079/150‑1 (MQ=255)                          gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTATCCT‑Caaaaaaaaaaaaaaa                                                                                                                         >  2:862825/1‑148 (MQ=255)                          gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCATTATCCT‑CTAaaaaaaaaaaaaa                                                                                                                         >  2:93040/1‑148 (MQ=255)                          gTCATTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAAc                                                                                                                         <  2:759089/149‑1 (MQ=255)                            cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCCATTAACCT‑Caaaaaaaaaaaaaa                                                                                                                          >  1:560057/1‑146 (MQ=255)                                 gAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACCCATTATACT‑CAAAAAAAAAAAAAACTTCCCGAtt                                                                                                                >  1:819977/1‑151 (MQ=255)                                     aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACCCCTTAGCCT‑CAAAAAAAAAACAAAATTCCAAAtttttt                                                                                                            >  1:217865/1‑151 (MQ=255)                                          aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTTACCT‑CCAAAAAAAAAAAAACATCCAGAATTATTTCtt                                                                                                        >  1:115608/1‑150 (MQ=255)                                          aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑a                                                                                                      <  2:456009/151‑1 (MQ=255)                                             aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTATGTAAAAACACACCTTATCCG‑CAAAAAAAAAAAAAAATTcccaaatttttattt                                                                                                        >  2:693775/1‑134 (MQ=255)                                               aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCaa                                                                                                 <  1:148283/151‑1 (MQ=255)                                               aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCaa                                                                                                 <  2:68639/151‑1 (MQ=255)                                                  aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGa                                                                                              <  2:11828/151‑1 (MQ=255)                                                     cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAACACTCATTAACCT‑CCAAAAAAAAACAAAACTCCAGATTTTTTACTT‑TAACAAAGaa                                                                                             >  1:66895/1‑149 (MQ=255)                                                     cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATa                                                                                           <  2:355949/151‑1 (MQ=255)                                                         cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑a                                                                                       <  1:772007/150‑1 (MQ=255)                                                          tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACACATTATCCT‑CCAAAAAAAAACAAAATTCCAAAATTTTTTTTT‑AACCAATTaaacaaaa                                                                                       >  2:319510/1‑144 (MQ=255)                                                          tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAAAACTCATTATACT‑CCAAAAAAAAAAAAAATTTCAGAATTATTTATT‑TACCAAAGAATAaa‑aa                                                                                      >  1:2642/1‑150 (MQ=255)                                                           ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑aa                                                                                      <  1:456499/149‑1 (MQ=255)                                                            tATCAATTTTCACTTAAAGTATGAACTGGATATCTTTGTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGc                                                                                    <  1:247723/150‑1 (MQ=255)                                                               cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTaa                                                                                <  1:122049/151‑1 (MQ=255)                                                                   tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACCCACATTAACCC‑CCAAAAAAAAAAAAAAATCCCGAATTTTTAATTAAACCAAT‑AATAAA‑AAAATTtaaaa                                                                             >  1:55414/1‑150 (MQ=255)                                                                   tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTATCCC‑CAAAAAAAAGAAAAAATTCCCGCATTATTACTT‑AAACAATAAAAAAC‑AAACTTaaaaa                                                                             >  1:169691/1‑150 (MQ=255)                                                                      cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACCACACATTAGCCT‑CCAAAAAAAAAAAAAATTCCAGATTTATTACTT‑AACCAAAAAAAAAA‑AAACTTaaaaaaaa                                                                          >  1:102396/1‑148 (MQ=255)                                                                      cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAatat                                                                         <  2:777529/151‑1 (MQ=255)                                                                          tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACt                                                                      <  1:156541/150‑1 (MQ=255)                                                                              gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCCATATCCT‑CCAAACAAAAACAAACTTCCCCATTTTTTTTTT‑TACCAATGAATAAA‑AAACTTAAAAAAAATACTGtt                                                                   >  1:39606/1‑149 (MQ=255)                                                                                aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCACTTACCC‑CAAAAAAAAAACAAACTTCCCGATTTTTTTCTT‑TACCAATGAATAAC‑AAAATTTAAAAAaataatttttt                                                                 >  1:717515/1‑142 (MQ=255)                                                                                  gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCATTATCCT‑CAAAAAAAAAAAAAAAATCCAGATTTATTAATT‑TACCAAAGAAAAAA‑AAAATTAAAAAAAATACTGTtttttt                                                              >  2:641547/1‑148 (MQ=255)                                                                                       ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGaa                                                        <  1:844381/151‑1 (MQ=255)                                                                                        gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCATTAGACT‑CAAAAAAAAAAAAAAATTTCAGATTTTTTACTT‑TAACAAAGAAAAAA‑AAACTTAAAAAAAATAATGTTTTTTTTAAGaa                                                        >  1:727758/1‑150 (MQ=255)                                                                                              tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTATTCA‑CAAAACAAAAACAAAATTCCAGAATTTTTATTT‑AACAAATAAATAAC‑AAACTTAAAAAAAAAAATTTTATTTTTAaaaaaaaa                                                    >  2:28661/1‑148 (MQ=255)                                                                                               ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAtt                                                  <  2:894159/149‑1 (MQ=255)                                                                                                   tCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCATTAACCTCCAAAAAAAAAACAAA‑TTCCAAAATTATTACTT‑AAACACTGAAAAAA‑AAAAATAAAAAATATAATGTTTTTTTTAaaaaaaaaaaaaaa                                              >  2:73936/1‑143 (MQ=255)                                                                                                    cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAgg                                            <  1:854135/150‑1 (MQ=255)                                                                                                          aaCTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTAGTCC‑CCAAAAAAAAAAAAAAATCCCTATTTTTTTCTT‑AAACAAACAATAACAAAACTAAAAAAAAATATTTTTTTTTTAAAAAAAAAATTaaacgga                                           >  2:883020/1‑141 (MQ=255)                                                                                                                  atcaACAGCGTATAGAGGCGGTTATGTAAAAACAAACATTAACCC‑CAAAAAAAAAACAAAACTCCCGATTTATTACTTAAACAAAAGAAAAAA‑AAATTTAAAAAAAAAAC‑GTTTTTTTTAAAAAAAAAATAAAAAGCATTAATTTAata                               >  1:415818/1‑149 (MQ=255)                                                                                                                     aaCAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTTACCT‑CAAAAAAAAGAAAAAAATTCAGATTTTTTTTTT‑AAACAAAGAAAAAA‑AAACTTAAAAAAAAAAATTTTATTTTAAAAAAAAAATTAAAAAGACTTAATTTAATAtt                             >  2:178377/1‑148 (MQ=255)                                                                                                                     aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaa                          <  2:812679/151‑1 (MQ=255)                                                                                                                        aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                       <  1:945518/151‑1 (MQ=255)                                                                                                                            tataGAGGCGGTTATGTAAAAACACTCCTTAGTCC‑CAAAACAAAAAAAAAATTCCAGATTTTTTACTT‑AACCAAATAATAAC‑AAACTTAAAAAAAATAATTTTTTTCTTAAAAAAAAAATAACAGGCATTAATTAAATATTaaaaaaaaa                    >  2:218509/1‑148 (MQ=255)                                                                                                                                gagGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa               <  1:919511/151‑1 (MQ=255)                                                                                                                                  ggCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt             <  1:749813/151‑1 (MQ=255)                                                                                                                                  ggCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt             >  2:905230/1‑151 (MQ=255)                                                                                                                                     ggTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa           >  2:704350/1‑150 (MQ=255)                                                                                                                                        tATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAtt       >  2:374747/1‑151 (MQ=255)                                                                                                                                         aTGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAtt       <  2:214232/150‑1 (MQ=255)                                                                                                                                           gTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTc     <  2:782737/150‑1 (MQ=255)                                                                                                                                           gTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt    >  2:98548/1‑151 (MQ=255)                                                                                                                                           gTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt    <  2:115594/151‑1 (MQ=255)                                                                                                                                             aaaaCCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa  <  1:36097/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      ACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTT‑AACCAATGAATAAC‑AAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTA  >  NC_000913/2903548‑2903838

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