Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,628,332	T→C	intergenic (+69/‑20)	ydfG → / → rspR	3‑hydroxy acid dehydrogenase YdfG/DNA‑binding transcriptional repressor RspR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,628,332	0	T	C	100.0%	137.4 / NA	38	intergenic (+69/‑20)	ydfG/rspR	3‑hydroxy acid dehydrogenase YdfG/DNA‑binding transcriptional repressor RspR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

AAGGCGATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACTCACTCACAAGCAAGAACGAATGACCGTCGAAACGCAACTTAATCCCACACAGCCTGTCAATCAGCAGATTTATCGTATTCTTCGTCGCGACATTGTCCATTGCCTGATTGCTCCAGGCACACCGTTGTCGGAAAAAGAAGTTTCTGTTCGTTTCAATGTGTCACGCCAGCCGGTTCGTGAAGCCTTTATTAAACTGGCGGAAAACGGCCTGATTCAAATTCGTCCGCAACGTGG  >  NC_000913/1628088‑1628566                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           aaGGCGATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCAATGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAAcccactca                                                                                                                                                                                                                                      <  2:103353/251‑1 (MQ=255)      gATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtac                                                                                                                                                                                                                                  >  2:18182‑M1/1‑247 (MQ=255)              tAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagc                                                                                                                                                                                                                         >  2:22543‑M1/1‑239 (MQ=255)                   cagcaaaaacctaTCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggc                                                                                                                                                                                                                      >  2:337690‑M1/5‑234 (MQ=255)                                 cAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtcc                                                                                                                                                                                                       <  2:110812‑M1/250‑31 (MQ=255)                                  aaaaTTCCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgt                                                                                                                                                                                                     <  1:392169‑M1/251‑33 (MQ=255)                                    aaTACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttc                                                                                                                                                                                                   >  2:162470‑M1/1‑217 (MQ=255)                                         cGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgagcagcggcaggacctttaggtccgttctcaat                                                                                                                                                                                              <  1:153769‑M1/251‑40 (MQ=255)                                                ttGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgtt                                                                                                                                                                                                    <  2:47142‑M1/238‑34 (MQ=255)                                                ttGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgtt                                                                                                                                                                                                    >  1:47142‑M1/1‑205 (MQ=255)                                                      ccAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaact                                                                                                                                                                                 <  1:22543‑M1/251‑53 (MQ=255)                                                               gtctgcgatgccggccggtgggtgtcaccgctgcctgctcCCGTCAATAACAATACCCTGGAAATGATGCCGTTTGCCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTACTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacaggaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgatt                                                                                                                                                                         <  2:244015‑M1/234‑61 (MQ=255)                                                                      aaGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATCCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttagggccgttctcaataccaaattcaacttcctgattctcagtta                                                                                                                                                                 <  2:388493‑M1/251‑69 (MQ=255)                                                                        gCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATTCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaatcccaaattcaacttcctgattctcagttaat                                                                                                                                                               <  2:131609‑M1/251‑71 (MQ=255)                                                                              gtgggggtgtcaaCGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttg                                                                                                                                                         <  2:350861‑M1/247‑77 (MQ=255)                                                                               ggtgggtgTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttg                                                                                                                                                         <  2:112089‑M1/250‑77 (MQ=255)                                                                                           cGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctct                                                                                                                                            <  1:441908‑M1/251‑90 (MQ=255)                                                                                           cGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaacgttttgaaatcgttgctct                                                                                                                                            <  2:20060‑M1/251‑90 (MQ=255)                                                                                               gcctgcTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaa                                                                                                                                         <  1:119686‑M1/250‑93 (MQ=255)                                                                                                                 caataCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttcacaataccaaattcaacttcctgattctcagttaatgatttgaaatcgttgctctgaattgctgagaaatagacaa                                                                                                                       >  2:109855‑M1/1‑140 (MQ=255)                                                                                                                          ggAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatcttt                                                                                                             <  1:337690‑M1/251‑121 (MQ=255)                                                                                                                          ggAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatcttt                                                                                                             >  1:222635‑M1/1‑131 (MQ=255)                                                                                                                          ggAAATGATGCCGGTGACCTAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAAATTTTATACCCTGCGTAACTGCCGTGTTATTGCTTGTCACAAAAAAGTTGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcatgacctttaggtccgttctcaataccaaattcaacttccagattctcagttaatgtttttaaatcgttgctctgaattgctgagaaatagacaaacgcatcttt                                                                                                             <  2:105229‑M1/251‑121 (MQ=255)                                                                                                                          ggAAAGGATGCCGGTTACCCATAGCTATGCCTGACTTAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCTGGTTATTTCTTGTCACAAAAAATTGGTAGACTCATGCATTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaatgcaacttcttgattctgagttaattttttgaaatcgttgctctgaattgctgagaaattgacaaactcatctt                                                                                                              <  1:394302‑M1/250‑120 (MQ=255)                                                                                                                                      ggTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatcttt                                                                                                 <  1:418261‑M1/251‑133 (MQ=255)                                                                                                                                                    aTGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaa                                                                                    <  1:441222‑M1/250‑146 (MQ=255)                                                                                                                                                                atttccaCCTTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcgtcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaaccttttt                                                                       <  2:256662‑M1/248‑159 (MQ=255)                                                                                                                                                                  gTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttc                                                                      >  1:11015‑M1/1‑91 (MQ=255)                                                                                                                                                                   tCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttc                                                                      >  2:101845‑M1/1‑90 (MQ=255)                                                                                                                                                                         gTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaa                                                              <  2:222635‑M1/251‑168 (MQ=255)                                                                                                                                                                                     tataCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttccct                                                   >  2:405247‑M1/1‑72 (MQ=255)                                                                                                                                                                                      ataCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatcttgcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaa                                                 <  1:149919‑M1/251‑181 (MQ=255)                                                                                                                                                                                             gCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtc                                          <  1:18182‑M1/251‑188 (MQ=255)                                                                                                                                                                                               gTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgcggagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcat                                        <  1:164195‑M1/251‑190 (MQ=255)                                                                                                                                                                                                          gTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagaca                             <  1:162470‑M1/251‑201 (MQ=255)                                                                                                                                                                                                           ttATTGCTTGTCACAAAAAAGAGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacacacacatctttgctgccatctttcggcgtgatgaaaccaacacctttttcagggttaaaccatttcactaaaccagtcattttgttagaca                             >  2:226365‑M1/1‑50 (MQ=255)                                                                                                                                                                                                                          aaaaaGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttg             >  1:125827‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                                    aGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaagaaatt     >  1:431097‑M1/1‑25 (MQ=255)                                                                                                                                                                                                                                     gACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggaactttaggtccgatctcaataccaaattcaacttcctgattctcagataatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaagaaattagc  >  2:70818‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           AAGGCGATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACTCACTCACAAGCAAGAACGAATGACCGTCGAAACGCAACTTAATCCCACACAGCCTGTCAATCAGCAGATTTATCGTATTCTTCGTCGCGACATTGTCCATTGCCTGATTGCTCCAGGCACACCGTTGTCGGAAAAAGAAGTTTCTGTTCGTTTCAATGTGTCACGCCAGCCGGTTCGTGAAGCCTTTATTAAACTGGCGGAAAACGGCCTGATTCAAATTCGTCCGCAACGTGG  >  NC_000913/1628088‑1628566

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