Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	343,442	G→A	G187R (GGA→AGA)	yahK →	aldehyde reductase, NADPH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	343,442	0	G	A	81.8%	74.8 / 0.4	33	G187R (GGA→AGA)	yahK	aldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (5/22);  total (11/22)
Fisher's exact test for biased strand distribution p-value = 4.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACTGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATG  >  NC_000913/343198‑343667                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     ctactGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTgggggatt                                                                                                                                                                                                                                >  2:63084/1‑246 (MQ=255)   actGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACata                                                                                                                                                                                                                                <  1:395031/249‑1 (MQ=255)         aTCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGACATATGGGGaat                                                                                                                                                                                                                        >  1:149437/1‑249 (MQ=255)                   ccgcgcttccgatctgCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGccc                                                                                                                                                                                                              <  2:57580/236‑1 (MQ=255)                                  gCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGAGACATATGGGGGTTAAGCTGGccccgatcgggagagcac                                                                                                                                                                                               >  1:57580/1‑236 (MQ=255)                                     gACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACAtg                                                                                                                                                                                            <  1:123153/251‑1 (MQ=255)                                        gCCGGACGAACCGGGCCCTACTCTGGGCGGCTCCTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCCTGGACATCTTGGGGTTAAACCGGCCCACCCGAgggggggacctttgg                                                                                                                                                                                         >  1:144877/1‑242 (MQ=255)                                        gCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATgtgg                                                                                                                                                                                         <  2:13522/251‑1 (MQ=255)                                                 aCCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑a                                                                                                                                                                               <  2:127884/251‑1 (MQ=255)                                                     ggCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGAGACATATTGGGGTTTAGCCGGGCCCCGCGAGGGGGGCACATTTGGTGGGAttttcccc                                                                                                                                                                            >  1:68154/1‑246 (MQ=255)                                                                gggCGGCTACTCACAACAGTTCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCGGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGGCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCa‑a                                                                                                                                                               <  2:239488/251‑1 (MQ=255)                                                                                  gATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAg                                                                                                                                             <  1:117310/251‑1 (MQ=255)                                                                                                aGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGGCCTATTGGGGTTAAACTGGCCCACGCCATGGGGGCACATGTGGTGGCATTTTACC‑CTTCTGAGGCA‑AAACGCGATGCGGCAAAAGACCCTggggccggt                                                                                                                               >  1:135738/1‑246 (MQ=255)                                                                                                    atatGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAg                                                                                                                          <  2:185919/251‑1 (MQ=255)                                                                                                             gCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTaa                                                                                                                   <  2:149415/249‑1 (MQ=255)                                                                                                                tATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGACATATGGGGGTTAAGCTGGCCCCCCCGGTGGGGGGGCATGTGGGGGCATTT‑TCCACATTTTTGGCA‑AAAAACGAAGCGGCAAAAA‑CCCCggggggccgtgaaggttgtaaatcac                                                                                                              >  2:175212/1‑230 (MQ=255)                                                                                                                 aTTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCAc                                                                                                              <  2:226671/250‑1 (MQ=255)                                                                                                                 aTTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCAc                                                                                                              <  1:67611/250‑1 (MQ=255)                                                                                                                                      cTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGGTTAAGCCGGCCCCCGCGGTGGGGGGAAATGTGGGGGGATTT‑ACCCCTTTTGAGGCA‑AAACGCGAAGCGGCAAAAA‑CCCTGGGGGGCGGTGAAGGTGTTTACTCCCGCAAAGCCGGTGAGATggcggg                                                                                        >  1:268762/1‑250 (MQ=255)                                                                                                                                          cggcggTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGGGGTGGCAGTT‑ACCACTGCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCGGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCAt                                                                                    <  1:232732/251‑1 (MQ=255)                                                                                                                                           ggcggTGGCTCCTTTGTTGTGTGCAGGGATCACCCCGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATc                                                                                   <  2:408726/251‑1 (MQ=255)                                                                                                                                               gTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGaa                                                                               <  2:22357/251‑1 (MQ=255)                                                                                                                                                 ggcgccTTTGTTGTGGGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAaga                                                                             <  2:119376/247‑1 (MQ=255)                                                                                                                                                            tgtgtgCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGAttt                                                                   <  1:59618/250‑1 (MQ=255)                                                                                                                                                                 gCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtttt                                                             <  1:343808/251‑1 (MQ=255)                                                                                                                                                                  cAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTg                                                            <  1:253193/251‑1 (MQ=255)                                                                                                                                                                      gATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATa                                                        <  1:63084/251‑1 (MQ=255)                                                                                                                                                                        tCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGAGACATATGGGGGGTAAGCTGGCCCACGCGATGGGGGCAAATGTGGGGGGATTT‑ACCAATTCTGGGGCA‑AAAAGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGGAGGTGTTAACTCACGCAAAGCCGATGGGATGGCGGCTCATCTGAAAAGTTTTGGTTTCATTTTGAATAc                                                       >  1:366600/1‑250 (MQ=255)                                                                                                                                                                            cacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAg                                                  <  1:370536/251‑1 (MQ=255)                                                                                                                                                                                tATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCCGGCCCCCCCGATGGGGGCACATGGGGTGGCATTT‑ACCACTTCTGAGGCA‑AAAAGCGAAGCGGCAAAAA‑ACCTGGGGGCCGATGAAATTGTTAACTCACACAATGCCCATGAGAAGGCGGCTCATCTGCAAAGGTTCTATTTTATTTTGAATACAGTAGCTgc                                              >  1:17965/1‑251 (MQ=255)                                                                                                                                                                                       cgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGAAATATGGGGGTTAAGCTGGGCCACGCGATGGGGGCGCATGTGGTGGCCTTT‑ACCCCTTCTGTGGCACAAAC‑CCAAGCGGCAAAAA‑CCCCGGGGGGCGGTGAAGGTGTTAAATCCCGCAATGCCGATGAGATGGGGGCTCCTCTTAAGAGTTTTTATTTTCTTTTGAATACACGAGCTGCGCCCCAt                                       >  2:186632/1‑251 (MQ=255)                                                                                                                                                                                                aCTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGGTTAAGCTGGGCCCCGCGGTGGGGGCACATGTGGGGGGATTT‑ACCACTTCTTAGGGA‑AAAAGCGAAACGGCAAAAA‑CCCCGGGGGGCGATGGAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGGTTTGCTTTTATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                              >  1:151744/1‑251 (MQ=255)                                                                                                                                                                                                         ccgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                     <  1:175212/251‑1 (MQ=255)                                                                                                                                                                                                                            gggCGTGGTCGGCATCGGCGGTCTGAGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATg  <  2:396291/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     CTACTGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTT‑ACCACTTCTGAGGCA‑AAACGCGAAGCGGCAAAAG‑CCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATG  >  NC_000913/343198‑343667

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: ‑