Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,413	0	A	C	38.9%	82.3 / 16.0	55	V370G (GTG→GGG)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (25/8);  new base C (0/21);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.12e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCAC  >  NC_000913/73271‑73561                                                                                                                                                  |                                                                                                                                                     gATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGAccc                                                                                                                                              >  1:152714/1‑151 (MQ=255) gATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGAccc                                                                                                                                              <  2:1094402/151‑1 (MQ=255)     gCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGCGTTGTTGAGCAGGAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGcc                                                                                                                                           >  1:272454/1‑150 (MQ=255)      ccATTAACGCATTGGTGAAAATCACAATCCCGTCGGCAGATGGTGGCAGCCCGATAGTGTTTTGGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCCCCAGCACCTGACCCGCCa                                                                                                                                          <  2:843154/150‑1 (MQ=255)              gCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCt                                                                                                                                 >  1:1137063/1‑151 (MQ=255)                   ggTGAAAATCAAAATGCCGTCAGCAGTTTGGGGCAGGCCGATAGTTTTTTTGAGAAGAAAAAAGAAGCCGTTAGCCAGCACAACCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCcg                                                                                                                             <  2:48332/150‑1 (MQ=255)                   ggTAAAAAACAAAGTGCCGCAAGCAGTTTGGGCAAGGCCGATCGGGTAGTTGAGAAGAAAAAAGAAGCGGGTAGCCACCACAACCCCGGGC‑ATGGCGAGG‑ATCAACAT‑GCCCCTCATCCCCAGCACCTGACCCGCCAGCATTTTCTg                                                                                                                                <  2:67649/147‑1 (MQ=255)                     ataaaaACACAGTGCCGAAAGGGGAGTGTGGCAGGCCAATAGTGTTTTGGAGAAGTAACAAGACGCGGGTACCCTTCCCATCCCCCGGC‑TGGTGGAGG‑ACCAAACG‑GCCGCCTCTCCCCACCACCGCAGCCGCCAGCATTTTCTGCcgcg                                                                                                                           <  2:1110955/148‑1 (MQ=255)                        aaaTCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGccc                                                                                                                        >  2:6362/1‑150 (MQ=255)                          aTCACAATGCCGTCAGCAGATTGTGGCAGGCGGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCa                                                                                                                     <  1:168428/151‑1 (MQ=255)                                 tGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgc                                                                                                               >  2:1003589/1‑150 (MQ=255)                                 ggctgTCACCATATTGTGGCGGGCCAATTGTTTTTTTGCGCAGTAAAAAAGGACCGGTGCCCCTCACCCCCCCGGTC‑CTGGCGGGT‑AACACTGT‑GCCACCCATCCCCACCATCCTTCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcg                                                                                                                <  2:1144546/145‑1 (MQ=255)                                  gCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgcg                                                                                                              >  2:28764/1‑150 (MQ=255)                                   ccGTCAGCAGATTGGGGCAGGCCGATAGTGTTGTTGAGCAGAAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc                                                                                                            <  1:896626/151‑1 (MQ=255)                                    cgtgggaagATTGTGGCGGGCGGTAAGTGTTTTTGACCGTAAAAAAGAAGCCGGTAGCCACCCCAATCCCCGGC‑TGGGCGAGG‑TTAAACTT‑GCCGCTCTTCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCa                                                                                                                     <  1:693579/134‑1 (MQ=255)                                     gTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTg                                                                                                          >  1:510821/1‑151 (MQ=255)                                      tCAGCAGATTGTGGCAGGCCGATAGTGCTGTTGAGCAGTAAAAAGAAGCCGGGAGCCAGCACAATCCCCGGC‑GTGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGc                                                                                                         >  2:235639/1‑151 (MQ=255)                                            gATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCAGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt                                                                                                    >  2:494429/1‑150 (MQ=255)                                              ttGTGGCGGGCCGATAGTTTTGTTGAGCATAAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt                                                                                                    <  1:432071/148‑1 (MQ=255)                                              ttGGGGCAGGCGGTTAGTTTTGTGGAGCATAAAAAAGAGCCCGGTACCCACCACACCCCCCGGC‑ATGGGGAGG‑ATCACCAT‑GCCGCTCATCTCCAGCCCCTAACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt                                                                                                    <  2:1016551/148‑1 (MQ=255)                                                  ggcgggcgGATAGTGTTGTGGAGCAGAAAAAAGAGGCCGGTAGCCACCACAACCCCGGGC‑AGGGCGAGG‑ATCAACAT‑GCCGCTCATCCCCAGCCCCGGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTg                                                                                                          <  2:628521/138‑1 (MQ=255)                                                  ggCAGGCCGATAGTGTTTTTGAGCAGAAAAAAGAAGCCGGTACCCAGCACAATCCCCGGC‑TTGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagca                                                                                                <  1:410629/148‑1 (MQ=255)                                                  ggCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑AGGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagcagc                                                                                              <  1:581268/150‑1 (MQ=255)                                                          gATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc                                                                                     >  2:1097362/1‑151 (MQ=255)                                                              gtgtTGTTGAGCAAAAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCa                                                                                    <  1:1083606/148‑1 (MQ=255)                                                                gttgttGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACt                                                                               >  2:819128/1‑151 (MQ=255)                                                                gttgtgGAGCAGAAAAAAGAAGCCGGTGGCCAGCACAACCCCCGGC‑TTGGCGAGG‑ATAAACAG‑GCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg                                                                                   <  2:261554/147‑1 (MQ=255)                                                                    ttgGGCAGAAAAAAGAAGCCGGTACCCACCACAACCCCCGGC‑ATGGGGAGG‑ATCAACAG‑GCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTAc                                                                             <  1:850557/149‑1 (MQ=255)                                                                      gAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcac                                                                         >  2:796957/1‑151 (MQ=255)                                                                      gAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑AGGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCCCCGGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc                                                                           <  1:220380/149‑1 (MQ=255)                                                                      aaaaataaaaaaaaaGCCGTCCGCCAGCAC‑ACCCCCGGA‑CTTGCCAGGAATAAACTG‑CCCGCGCATTTCCAGCTCCCGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGca                                                                          <  2:176375/143‑1 (MQ=255)                                                                        gCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcaca                                                                        >  2:16431/1‑150 (MQ=255)                                                                        cgggaaaaaaGAAGCGGCTGGCAAGCAC‑ATCCCCGGCAAGGGCGAGG‑A‑CAAAATCGCCATAATTCCACAGCCCCAGGCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACt                                                                               <  2:594251/138‑1 (MQ=255)                                                                             aaaaaGAACCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATacac                                                                  <  1:901197/151‑1 (MQ=255)                                                                              aaaaGAACCCGGTAGCCAGCACAACCCCCGGC‑ATGGCGAGG‑ATCAACAG‑GCCGCTCACCCCCAGCACCGACCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACAcc                                                                 <  1:364546/151‑1 (MQ=255)                                                                                    agCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTgccgc                                                           >  2:1051114/1‑151 (MQ=255)                                                                                     gCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTgccgcc                                                          >  1:1098422/1‑151 (MQ=255)                                                                                      ccGGTAGCCGCCACAATCCCCGTG‑ATGGCGGGG‑ACAACCAT‑GCCGACCATCACCAGCCCCGGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCa                                                         <  2:441474/151‑1 (MQ=255)                                                                                      ccGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCa                                                         <  1:1075900/151‑1 (MQ=255)                                                                                          tAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgcgc                                                     >  1:468534/1‑151 (MQ=255)                                                                                           aGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgcgc                                                     >  2:1096079/1‑150 (MQ=255)                                                                                             ccaccaCAACCCCGGGC‑AGGGCGAGC‑ACCAACAT‑CCCCCCCTCCCCCGCCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgcgc                                                     <  2:416524/148‑1 (MQ=255)                                                                                               aGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTCCGCACAATACACCTGCCGCCAGCGCAATAc                                                >  1:690113/1‑151 (MQ=255)                                                                                                    accccccGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGc                                              <  1:1152297/145‑1 (MQ=255)                                                                                                    aaTCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAAc                                           >  2:1062964/1‑151 (MQ=255)                                                                                                       cccgGGC‑AGGGCGAGA‑ATCAACAT‑GCCGCTCACCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAg                                        <  1:743023/151‑1 (MQ=255)                                                                                                                tGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcac                                  >  2:615945/1‑149 (MQ=255)                                                                                                                 ggggAGG‑ATCAACAG‑GCCGCTCATCCCCAGCCCCGGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcaca                                 <  2:875541/146‑1 (MQ=255)                                                                                                                 ggCGAGG‑ACCAACAG‑GCCGCTCACCCCCAGCCCCGGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcac                                  <  2:98919/148‑1 (MQ=255)                                                                                                                   cGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgcgc                             >  2:224410/1‑151 (MQ=255)                                                                                                                      gg‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTg                          >  2:210607/1‑151 (MQ=255)                                                                                                                      gg‑ATCAACAT‑GCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCt                           <  1:471050/150‑1 (MQ=255)                                                                                                                                    cGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCgg              >  1:161948/1‑151 (MQ=255)                                                                                                                                     gCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGt             >  1:1046588/1‑151 (MQ=255)                                                                                                                                                gCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCAc  >  2:370543/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGC‑ATGGCGAGG‑ATCAACAT‑GCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCAC  >  NC_000913/73271‑73561

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