Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	41.8%	86.2 / 13.3	55	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (9/23);  new base G (23/0);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.98e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCG  >  NC_000913/1157125‑1157420                                                                                                                                                   |                                                                                                                                                          atCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTaaa                                                                                                                                                        >  2:346861/1‑150 (MQ=255)   ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAtt                                                                                                                                                      <  2:436974/150‑1 (MQ=255)   ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTa                                                                                                                                                     <  2:517459/151‑1 (MQ=255)    cAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTaaaat                                                                                                                                                      >  2:140058/1‑147 (MQ=255)    cAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAtt                                                                                                                                                      >  1:467298/1‑149 (MQ=255)       aTTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTgg                                                                                                                                                 <  2:32087/151‑1 (MQ=255)         tGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGAt                                                                                                                                               >  2:396080/1‑151 (MQ=255)          ggTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTatt                                                                                                                                                   >  2:146178/1‑144 (MQ=255)             cGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTTGGtc‑t                                                                                                                                            >  1:843127/1‑149 (MQ=255)             cGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATaattggtatcg                                                                                                                                           >  1:245687/1‑141 (MQ=255)                     gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTa                                                                                                                                   <  1:132380/150‑1 (MQ=255)                        cAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTATTgggttttggggttt                                                                                                                                     >  1:450387/1‑134 (MQ=255)                            ccGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑t                                                                                                                         <  2:401905/151‑1 (MQ=255)                                  aTCGTGCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGAGTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTcc                                                                                                                    <  1:583726/150‑1 (MQ=255)                                  aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTcc                                                                                                                    <  1:815330/150‑1 (MQ=255)                                     gTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTCTTGATC‑TCGGGTTTTAAAT‑ATCC‑tttt                                                                                                                      >  1:781326/1‑145 (MQ=255)                                     gTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGATC‑TCGGGTTTTTCAT‑CTTC‑tttt                                                                                                                      >  1:937818/1‑145 (MQ=255)                                       ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCAt                                                                                                               <  1:187030/150‑1 (MQ=255)                                        cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATtg                                                                                                             <  2:686840/151‑1 (MQ=255)                                        cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATC‑TTGGATTTTATAT‑CTTC‑TTTTTCGGCAtt                                                                                                              >  1:402813/1‑150 (MQ=255)                                         aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATACTGGATC‑TCGGGATTTTCAT‑CTTC‑tttttctggctt                                                                                                              >  1:567693/1‑141 (MQ=255)                                         aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATTCTTGGTT‑TCGGGTTTTAAAT‑TTCC‑TTTTTCGGCAtt                                                                                                              >  1:1022864/1‑149 (MQ=255)                                             ccGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTC‑TCGGGTTTTAAAA‑CCCC‑CTTTTCGGCATTTTGaac                                                                                                        >  2:535481/1‑149 (MQ=255)                                                 gACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGGTTTTACAT‑CTCC‑TTTTTCCGCATTGTGACCTTcc                                                                                                    >  2:373624/1‑151 (MQ=255)                                                  aCGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑CCGGATTTTACCT‑CTCC‑TTTTCCGGCATTGTGCCCTAcc                                                                                                    <  1:522726/150‑1 (MQ=255)                                                    gCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGt                                                                                                  <  2:162590/150‑1 (MQ=255)                                                          gCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACATAGGACAGAGAAACGGCGGGGAAAATTATGGGTG‑TTGGGTTTTAAATCCCCC‑TTTTCGGGCATT‑TGGCCCTCCGCAATGCgg                                                                                           >  2:118080/1‑151 (MQ=255)                                                           ccGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATT‑TCGGATTTTAAAT‑CTCC‑TTTTTCGGCATTGTGACCTTCCGTAATGCgg                                                                                           >  1:389485/1‑150 (MQ=255)                                                              aTACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTC‑TCGGGTTTTACAT‑CTCC‑TTTTTCCGCCTTTTGACCTTTCGTAATGCgggg                                                                                         >  2:564719/1‑147 (MQ=255)                                                               tACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATT‑TCGGGTTTTAAAT‑CTCC‑TTTTTCCGCCTTTTGACCTTCCGTAAAGCgggg                                                                                         >  2:605517/1‑146 (MQ=255)                                                                    tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTC‑TCGGGTTTTAAAT‑TTCC‑TTTTTCGGCCTTTTGACCCTCCGTAAAGCgggg                                                                                         >  2:893690/1‑141 (MQ=255)                                                                     ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgc                                                                                  <  1:727588/149‑1 (MQ=255)                                                                                gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATAATGGATC‑TCGGATTTTAAAT‑CTTC‑TTTTTCGGCATTTTGAGCTTCCGTAAATCGGGGGAAATGCGCGATGCTGcgcg                                                                     >  2:512858/1‑151 (MQ=255)                                                                                      aaaaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTAATGGATC‑TCGGGTTTTAAAT‑CTCC‑TTTTCCGGCCCTGTTACCTTTCCTAATGCGGAGGAAATGCGCCGTGCTGcgcg                                                                     >  1:904913/1‑145 (MQ=255)                                                                                             gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCt                                                         <  2:1023716/150‑1 (MQ=255)                                                                                             gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑g                                                       <  2:593513/151‑1 (MQ=255)                                                                                                ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGAt                                                    <  1:89096/151‑1 (MQ=255)                                                                                                  aTTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCg                                                  <  2:345468/151‑1 (MQ=255)                                                                                                  aTTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATAATGGGTC‑TCGGAGTTTACAT‑CTCC‑TTTTTCGGGATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCg                                                  >  2:31768/1‑151 (MQ=255)                                                                                                    tGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATC‑TCGGGTTTTAAAT‑CTCCTTTTTCCGGC‑TTGTGGCCTTCCGTAATGCGGAGCAAATGCGCGCTGCTGCGCGTTTTTTTCCCCT‑GGAGCgg                                                 >  1:75199/1‑150 (MQ=255)                                                                                                     gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATC‑TCCGCTTTTAAAT‑CTCC‑CTTTTCGGCCATGTGTCCTTCCCTAATGCGGGGCAACTGCGCGGTGCTTCGCGTTTTGTTCCCCC‑GGGTCgg                                                 >  1:611528/1‑149 (MQ=255)                                                                                                            cGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACtggtg                                        <  2:47920/151‑1 (MQ=255)                                                                                                                 tacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTATGGGTCTTGGGATTTTA‑AT‑CTCC‑TTTTTCGGCGTTTTGATCTTCCCTAAATCGGTGGAAACGCGCGATGCTTTGCGTTATGTTCCCCC‑CGATCGGGTTTtggtggg                                      >  1:137704/1‑147 (MQ=255)                                                                                                                  acacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACt                                  <  1:970913/151‑1 (MQ=255)                                                                                                                   cacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTg                                 <  2:798548/151‑1 (MQ=255)                                                                                                                         tttACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATC‑TCGGATTTTAAAT‑TTCC‑TTTTTCGGCCTTGTGACCTTCCGTAAAGCGGAGCAACTGCGCGGTGCTGCGCGGTATTTCCCCCC‑GGGTCGGTTACTGGGGGGAACTGACTCac                           >  2:71741/1‑151 (MQ=255)                                                                                                                           tACAGAGGACAGAGAAACGGCGGGTAAATTACTGGGTC‑TCGGGTTTTTCAT‑CTCC‑TTTTTCGGCATTGTGACCTTTCCTAATTCGGAGCAACTGCGCGATGCCGCCCGTTTTTTTCCCCC‑GGGTTGGTTAATGGGGGGAAATGACTCaac                          >  2:397068/1‑148 (MQ=255)                                                                                                                           tACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCacc                          <  1:192115/150‑1 (MQ=255)                                                                                                                           tACAGAGGACAGAGAAACGGCGGGTAAATTAATGGGTT‑TCGGGTTTTATAT‑CTCC‑TTTTTCGGCATTGTGACCTTCCGTTATTCGGAGGAACTGCGCGATGCTGCGCGTTTTTTTCCCCC‑GGATCGGGTACTGGGGGAAACTGACCCacc                          >  1:118886/1‑150 (MQ=255)                                                                                                                             cAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCacctta                       <  2:129715/151‑1 (MQ=255)                                                                                                                              agagGACAGAGAAACGGCGGGGAAAATACTTGGTC‑TCGGGTTTTTAAT‑ATCC‑CTTTTCCGCACTGTGAACTTTCCTAATGCGGGGCAACAGCGCGGTGCTGCGCGTTTTGTTCCCCTCGGTTGGGTTCTTGGGGGAAAC‑GAATCacc                          >  1:791370/1‑147 (MQ=255)                                                                                                                                    cAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTgcg                <  1:595920/151‑1 (MQ=255)                                                                                                                                              gcggGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGc      >  2:547164/1‑151 (MQ=255)                                                                                                                                              gcggGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGc      >  2:759446/1‑151 (MQ=255)                                                                                                                                               cggGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCa     <  2:245690/151‑1 (MQ=255)                                                                                                                                                   tAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCg  >  1:414871/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                          ATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATC‑TCGGATTTTACAT‑CTCC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCT‑GGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCG  >  NC_000913/1157125‑1157420

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