Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,483,254	0	G	A	47.3%	46.6 / 35.1	58	S29L (TCA→TTA)	emrK	multidrug resistance efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (0/26);  total (18/40)
Fisher's exact test for biased strand distribution p-value = 8.06e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTT  >  NC_000913/2483105‑2483400                                                                                                                                                      |                                                                                                                                                    cTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTg                                                                                                                                                     <  2:474851/150‑1 (MQ=255)       gAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGc                                                                                                                                              <  1:112322/151‑1 (MQ=255)       gAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGc                                                                                                                                              >  1:996903/1‑151 (MQ=255)        aaCGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCa                                                                                                                                             <  1:850730/151‑1 (MQ=255)         aCGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaa                                                                                                                                            <  2:946953/151‑1 (MQ=255)           gTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaata                                                                                                                                          <  1:368270/151‑1 (MQ=255)              gTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaat                                                                                                                                           >  2:250897/1‑147 (MQ=255)                        ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTa                                                                                                                               <  1:757717/148‑1 (MQ=255)                         tATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCg                                                                                                                            <  1:999430/150‑1 (MQ=255)                           tGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGcc                                                                                                                          >  1:1193091/1‑150 (MQ=255)                            gATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCa                                                                                                                         <  1:996648/150‑1 (MQ=255)                              ttAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGGCCAATAGGCATAGGCACCTCAAAACGCAATAAATAAAA‑CTACCGCCAAt                                                                                                                       >  2:677525/1‑150 (MQ=255)                                       gTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACAAATAGGCATAGGCACCTGAAAACGCAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGt                                                                                                             <  2:925093/151‑1 (MQ=255)                                                ccTGAGTTTTGTGCAGAATTTGACTCTGCTTTCCCGTTAAAATAGTCGTCATCTGTACGAATCTTTTCTTCAAATTCCATGGCCAAATAGGCATAGCCCCCTAAAAACGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTttcttc                                                                                                    <  2:526867/151‑1 (MQ=255)                                                ccTGAGATTTGTGCAAAAATTGGATCTACATCCCCCGAGAAATAGGCGTCATCTGTACTAATCATCTTTTCAAATCCCATTCACCAATAGACATAGCCACCTAAAAACGCAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTttctt                                                                                                     <  1:545578/150‑1 (MQ=255)                                                   gagaCTTGTGAAAAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCTTGTCTTCTAATTCCATTGACCAATAGGAAAAGGCACCTAAAAACGAAAAAAATAAAC‑CTACCGCCAATAAAGAAAAGTATTTtcttct                                                                                                   <  2:51894/149‑1 (MQ=255)                                                        ttgtaaaaaaaTTGATTCTGCTTTCCCCGTAACATAGGCTTCTTCTTTTAAAAACATGTTTTATAATTCCATGACCCAAAAGGAAAAGACCCCT‑AAAACCCAAAAAAAAAAATCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                                                                              <  1:661957/141‑1 (MQ=255)                                                        ttGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                             >  2:713632/1‑150 (MQ=255)                                                        ttGTGCAAAAATTGGTTCTCCTTTCCCCGTGAAATGGGCGTCATCTGTACAAATCATGTCTTCTAATCCCATTGCCAAAGAGGCATAGGCCCCTAAAACCGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                                                                              <  2:114593/149‑1 (MQ=255)                                                        ttGTGCAAAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATTTGTACTAATCATTTCTTATAATCCCATTCACCAATAGACATAGGCCCCTAAAACCGCAAAAAAAAACA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                             <  1:953523/150‑1 (MQ=255)                                                         tgtgCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAAt                                                                                           >  2:133799/1‑151 (MQ=255)                                                           tgCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGt                                                                                         >  2:1109616/1‑151 (MQ=255)                                                              aaaatttGTATCTGTATTCCCCGTAACATAGGCGTCATTTGTACAAATAATGTTTTATAATTCCATTGACCAATAGGAAAAGGCACTTAAAAACGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                                                                            <  1:77652/140‑1 (MQ=255)                                                              aaaatttGGATCTCCATTCCCCGTGACATAGGCTTCATTTGTATTAATCTTTTTTTCTAATCCAATTAACCAAAAGGCATGGCCCCCTAAAACCGAAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                      <  1:906928/146‑1 (MQ=255)                                                              aaaaatgctaTCCGCATCCCACGGGACGTGGGCTTCATCTGTTCTAACCATTTTTTCTAATCCCTTTCCACAATAGAAAAAGCACCCTAAAACCGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGttt                                                                                       <  1:758214/141‑1 (MQ=255)                                                              aGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATCCCATTGACAAATAGGCATAGGCACCTAAAAACGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGt                                                                                         <  2:435303/148‑1 (MQ=255)                                                                aaaTTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttttt                                                                                    >  1:1008314/1‑151 (MQ=255)                                                                   ttGGATCTGCTTCCCCCGTGACATAGGCGTCATCTGTACTAATCATTTTTTCTAATTCCATTGCCAAATAGGCATAGCCACCTAAAACCGAAAAAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                                                                                  <  2:4016/150‑1 (MQ=255)                                                                     ggATCAGAATCCCCCGTACAAAAGGCGTCTTCTGTACAAATCATTTCTTATAATCCCATGGCCAAATAGAAAAAGCCACCAAAAAGCGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTa                                                                                   <  1:370457/147‑1 (MQ=255)                                                                       aTCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGa                                                                              >  1:115025/1‑150 (MQ=255)                                                                              tcccgcGTGACTTGGGCTTCATCTTTACTAATCTTTTCTTCTAATTCTATACCCAAATGGAGATGGGCACTTGAAAACGCAAAAAAGAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTa                                                                          <  2:1160886/142‑1 (MQ=255)                                                                              tcccccGTGACATAGGCGTCATCTGTACTAATAATTTCTTAAAATTCCATTGACCAATAGGCATAGCCACCTAAAACCGAAATAAAAAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATc                                                                       <  1:82657/148‑1 (MQ=255)                                                                              tcccccGTGACATAGCCGTCATCTGTACAAAACATGTTTTCTAATCCGATTGACAAAGAGGCAGTGGCGCCAGAAAGCGCAATAAATAAAC‑CTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAt                                                                            <  2:404694/143‑1 (MQ=255)                                                                              tcccccGAGAAATAGGCGCCATTTTTAATATTCATGTTTTTTAATTCAATTGACCAATAGTAAAAGCCACATAAAATCGAAATAAAAAAAA‑TTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                           <  1:617456/144‑1 (MQ=255)                                                                                ccccGTGACATAGGCGTCATCTTTACTAATCATTTCTTCTAATTCCATTGACCAATAGGCATAGGCCCCGAAAAACGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCt                                                                      <  2:918814/149‑1 (MQ=255)                                                                                ccaggggACATAGGCTCCATCTGCACAAATTATGTTTATTAATCCCATTGTCCAAGAGGAATAGGCACGTCAAACCGCAATAAATAAAC‑CTGCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCt                                                                      <  2:983085/143‑1 (MQ=255)                                                                                    gTGACAAAGGCGTCATCTGTACTAATCATGTCTTCTAATTCAATAGACAAAAAGGAATAGGCACCTAAAAACGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc                                                                 <  1:686059/150‑1 (MQ=255)                                                                                       aCATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt                                                              >  2:944349/1‑150 (MQ=255)                                                                                         agaGGCGTTATCTATTCAAATCATGTCTTCTAATCCCATGGACAAATGGACAGAGGCCCAAAAAACCGAAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtat                                                            <  2:913900/148‑1 (MQ=255)                                                                                         aTAGGCTTCATTTTTACAAACCATTTCTTCAAATTCCATTGACCAATAGGCATAGGCCCCTGAAACCGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                           <  1:86735/151‑1 (MQ=255)                                                                                            ggCGTCTTCTGTACTAATAATTTCTTTAAATTCCATTGACCAATAGACAAAGGCACCTAAAAACGCAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt                                                              <  1:233688/145‑1 (MQ=255)                                                                                                tattgtatACAAATAATGTCTTATATTACCATCAACAAAAAGAAAAAGGCACATAAAAACTAAATAAAAAAAA‑TTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtattat                                                         <  1:152914/139‑1 (MQ=255)                                                                                                     tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                               >  1:681291/1‑151 (MQ=255)                                                                                                     tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                               >  1:1115849/1‑151 (MQ=255)                                                                                                     tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                               >  1:749300/1‑151 (MQ=255)                                                                                                      taccTAAACATGTCTTCTATTTCCTTTAACAAATAGGCATAAACACCGAAAATCGAAAAAAATAATA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAt                                                  <  1:596366/144‑1 (MQ=255)                                                                                                      gTACTAATAATTTTTTCAAATTCCATTGACAAATAGGAATAGGCACCTAAAAACCAAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAttt                                                <  2:558887/149‑1 (MQ=255)                                                                                                        aCTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtctc                                             >  1:1121328/1‑150 (MQ=255)                                                                                                              cATGTTTTCTAATCCCATTAACCAATAGGCATAGGCCCCTAAAAACAAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGa                                        <  2:892192/149‑1 (MQ=255)                                                                                                              cATGTCTACTAATCCCATTGACCAATAGGCATAGGCACCAGAAAACTAAATAAATAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGa                                        <  1:1123966/149‑1 (MQ=255)                                                                                                                tGTCTTATAATCCCATTCACCAATAGACAAAGGCACTTAAAACCGAAAAAAAAAAAA‑TTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATg                                      <  1:172546/149‑1 (MQ=255)                                                                                                                    ttctAATTCCATGGCCAAATAGTAAAAGGCACCTAAAACCACAAAAAATAAAA‑CTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGa                                     <  2:235007/146‑1 (MQ=255)                                                                                                                                tGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattgaat                     >  1:168361/1‑150 (MQ=255)                                                                                                                                 gACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatg                    >  2:231167/1‑150 (MQ=255)                                                                                                                                        aggcataggcaCCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGc            >  1:508062/1‑151 (MQ=255)                                                                                                                                         ggcataggcaCCAAAAAACGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGca           <  1:946394/151‑1 (MQ=255)                                                                                                                                              agccaCCTAAAACCGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGca           <  1:274400/143‑1 (MQ=255)                                                                                                                                                 caCCTAAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttct   <  1:517139/151‑1 (MQ=255)                                                                                                                                                   ccTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttctt  >  2:609116/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTT  >  NC_000913/2483105‑2483400

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