Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,177	0	A	C	50.0%	88.1 / 10.2	58	G308G (GGT→GGG)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (27/2);  new base C (0/29);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 3.09e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTG  >  NC_000913/2921028‑2921320                                                                                                                                                      |                                                                                                                                                aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccac                                                                                                                                                >  2:717586/1‑151 (MQ=255)      cATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                                                                                           >  1:411273/1‑151 (MQ=255)      cATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                                                                                           >  1:470051/1‑151 (MQ=255)                 tgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGAt                                                                                                                                >  2:723523/1‑151 (MQ=255)                        ccACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaacc                                                                                                                          >  2:47830/1‑150 (MQ=255)                          aCTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                        >  1:787362/1‑150 (MQ=255)                             cAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaac                                                                                                                    >  1:850445/1‑151 (MQ=255)                             cAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaa                                                                                                                     >  1:713331/1‑150 (MQ=255)                                 gTTGACGTGGTTGCAGACCCCCATCCCCATCGACAGCACCATGCCCATTCCGACCGGTTTTTCACGCGCAATGTCCAGCGATCCCGGGCGGCGCACCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                                  <  2:263971/149‑1 (MQ=255)                                       gtgtttGCAGACCCCCATCCCCACCGCCAGCACCAGGCCCATTCCGATGGTTGTTTTACGGGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                          <  2:295487/147‑1 (MQ=255)                                           ttGCAGAACCCCACCCCCATCGACAGCAACATGCTCATTACGATCGGTGTTTTACGCGAAATGTCCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAATCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                       <  1:830934/150‑1 (MQ=255)                                           ttGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgga                                                                                                        >  1:1016970/1‑149 (MQ=255)                                           tggCAGAACACCATCACCATCGACAGCAACATGCCAATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaa                                                                                                            <  1:870271/143‑1 (MQ=255)                                                   caccatccccatCGACAGCAACATCCCCATTACGACCGGTGTTTTACGCGAATTGTTCAGCGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaaa                                                                                               <  1:622588/150‑1 (MQ=255)                                                     ccaccaccatcGACAGCAACATCCCCATTACGATCGGTTTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGcc                                                                                            <  1:196835/151‑1 (MQ=255)                                                      catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                           >  2:34022/1‑151 (MQ=255)                                                         ccccagcGACAGCAACATCCCCTTTCCGACCGTTTTTTTACGCAAAATGTCCAGCGACCCCGGGCGGCCCATCACCCAGTCGAAAATAACCCCCCCCAACACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCt                                                                                        <  1:1042476/149‑1 (MQ=255)                                                            catcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                     >  1:401947/1‑151 (MQ=255)                                                               cGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                     >  2:1062532/1‑148 (MQ=255)                                                                   accAACATCCCCATTCCGATCGGTGTTTAGCGCAAAATGTCCGGCGACCCGGTGCGGCACACCAGCCACGCGAAAATACCCCCCCCCGACCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGaa                                                                                 <  1:944208/146‑1 (MQ=255)                                                                   aGCAACATGCCCATTACGATCGGTGTTTTACGCGCACTGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCg                                                                              >  2:218372/1‑151 (MQ=255)                                                                      accttGCCCATTCCGTTCGGTTTTTTACGCGCAATGTTCAGCAATCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                               <  1:1014325/143‑1 (MQ=255)                                                                      aaCATGCCCATTACGATCGGTTTTTTACGCGCAATTTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACa                                                                           <  1:422315/151‑1 (MQ=255)                                                                        cATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATc                                                                         >  1:675637/1‑151 (MQ=255)                                                                            cccATTCCGACGGGTGTTTACCGGGAAATGTCCGCCGCCCCCGTGGGGCGCACCGCCCAGTCGAAAATAACCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                     <  2:769349/151‑1 (MQ=255)                                                                                ttACGTTCGGTGTTTAACGCGAAATGTTCAGCGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGt                                                                   <  1:611213/149‑1 (MQ=255)                                                                                  aCGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGcc                                                                >  2:316321/1‑150 (MQ=255)                                                                                     aTCGGGGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCAGCTTTCAGAATCGACATCCCGCGTGCCTg                                                              >  2:360887/1‑149 (MQ=255)                                                                                      tCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAc                                                           >  2:508835/1‑151 (MQ=255)                                                                                      tCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAc                                                           >  2:941967/1‑151 (MQ=255)                                                                                        ggTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                          >  1:538382/1‑150 (MQ=255)                                                                                            ttttACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGAt                                                      >  2:1029063/1‑150 (MQ=255)                                                                                                      aaTTTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                                <  1:34022/146‑1 (MQ=255)                                                                                                      aaTTTTCAGCGATCCCGGGCGGCCCACCAGCCAGTCGGAAATAACCCCCCCCAACACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                               <  2:165245/147‑1 (MQ=255)                                                                                                      aaTGTCCAGCGACCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                             <  1:529888/149‑1 (MQ=255)                                                                                                          tccAGCGACCCCGGGCGGCGCATCAGCAAGTGGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCa                                         <  1:420277/147‑1 (MQ=255)                                                                                                            cAGCGACCCGGGGCGGCCCACCCCCCAGTCGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCa                                         <  2:507204/147‑1 (MQ=255)                                                                                                               cGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                   >  1:939910/1‑150 (MQ=255)                                                                                                                gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGgaata                                  >  2:118473/1‑150 (MQ=255)                                                                                                                   ccggggcgcccaaccagccagcCGAAAAAAACCCCCCCCACACCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                                 <  2:1513/134‑1 (MQ=255)                                                                                                                   cccGTGCGGCGCATCAGCCAGTCGAAAAAAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAgg                                       <  1:674414/142‑1 (MQ=255)                                                                                                                     cGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                              >  1:401926/1‑149 (MQ=255)                                                                                                                      gggcggcaTATCAGCCAGTGGGAAAAAATCCCCCCCAGAACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                              <  1:620318/146‑1 (MQ=255)                                                                                                                      gTGCGGCGCACCAGCCAGTCGAAAATAACCCCACCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                           <  2:650363/151‑1 (MQ=255)                                                                                                                          ggcgCATCCGCCAGCCGAAAATAACCCCACCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                       <  2:953919/151‑1 (MQ=255)                                                                                                                             gcATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                      >  2:869851/1‑149 (MQ=255)                                                                                                                                 cacccggcgaaaaaTAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg                <  2:857508/140‑1 (MQ=255)                                                                                                                                 cAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg                >  1:981275/1‑151 (MQ=255)                                                                                                                                 cAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAgg                   <  1:673493/148‑1 (MQ=255)                                                                                                                                    ccggTCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg                <  1:497373/145‑1 (MQ=255)                                                                                                                                    ccAGTCGAAAAAACCCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg                  <  1:527142/146‑1 (MQ=255)                                                                                                                                    caaGTCGAAAATACCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg                  <  2:1007959/144‑1 (MQ=255)                                                                                                                                       gccgaaaaTAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt           <  2:291987/145‑1 (MQ=255)                                                                                                                                          aaaaaaaaCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc         <  2:164911/147‑1 (MQ=255)                                                                                                                                            aaaTAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa        <  2:790442/148‑1 (MQ=255)                                                                                                                                              ataacCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt   <  2:182318/151‑1 (MQ=255)                                                                                                                                               taatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg  >  1:110368/1‑151 (MQ=255)                                                                                                                                                aatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg  >  2:635167/1‑150 (MQ=255)                                                                                                                                                aatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg  >  1:589033/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTG  >  NC_000913/2921028‑2921320

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