Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	97.1%	123.7 / NA	35	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (13/20);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATT  >  NC_000913/3485962‑3486440                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              gTAATACATTGATGTACTGCATGTAGGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAGAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAtt                                                                                                                                                                                                                                        <  2:218920/251‑1 (MQ=255)  tAATACATTGATGTACTGCATGTATGCAAAGGACGTAACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTc                                                                                                                                                                                                                                       >  2:537004/1‑251 (MQ=255)                 cgtgaTGTATTCAAAGGACGTCACATTCCCGTGCAGTAAAGTTGATAGCCCCTGCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGGTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAGCCGCGCATGTTGCTTGGCAAACCGCAAACAGGCCCGACTCTCGAATGGTTCGTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAg                                                                                                                                                                                                                       <  1:530024/247‑1 (MQ=255)                       gTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGaa                                                                                                                                                                                                                  >  1:451623/1‑250 (MQ=255)                           gCAAAGGCCGTCACATTCCCGTGCAGTACAGTTGATACCCCCTTCCCAGGTAGCGGGAATCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGCCTCGGGAGAAAGCTTATAACATAGGATAACCGCTCATGGTGCTTGTCAAACCGCAACCAGCCCCGACTCTCGTATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGATCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCt                                                                                                                                                                                                             <  1:1425/251‑1 (MQ=255)                                      cacaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATc                                                                                                                                                                                                   <  1:367935/250‑1 (MQ=255)                                             ccGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAgg                                                                                                                                                                                           <  1:73115/251‑1 (MQ=255)                                                      aCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCGGGAGAAAGTTTATAACAGAGGATAACCGCGCGTGGGGCTTGGCAAACCTCAAACAGACCCGACGCTCGAATGGTTCTTGTCTCATTGCCACAGTCATAAGTACCCATCCAGGAGCAAGCTTAGTCACCAGGT‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTgg                                                                                                                                                                                   <  1:318078/250‑1 (MQ=255)                                                            gATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGtactac                                                                                                                                                                                                      >  1:25740/1‑225 (MQ=255)                                                            gATAGCCCCTCCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACGCTCGAATGGTTCTTGTCTCATTGCCACATGCATAATTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCt                                                                                                                                                                            <  1:45292/251‑1 (MQ=255)                                                             aTAGCCCCTCCCCAGGTGGCGGGAAGCATATTCCGGCAATCCGGAGACAGCGGCTTTATCGGGCTCTGGAGAAAGCTTAGAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGTCTCTGTGGCAGTGCt                                                                                                                                                                            <  2:70364/250‑1 (MQ=255)                                                                   ccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaa                                                                                                                                                                      >  2:210599/1‑250 (MQ=255)                                                                        ccAGGTAACGGGCAGCGTATTTCGGCACTCCAGAGACACCGGCGTTATCTGGATCTGGAGAAAGCTTATACCAGCGGATAACCGCGCATGGCGCTTGGCAAACAGCCAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTACAGGCTCTGTGGCAGTGCTGATCACAGa                                                                                                                                                                   >  2:63110/1‑248 (MQ=255)                                                                                  ggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaa                                                                                                                                                       >  2:10047/1‑250 (MQ=255)                                                                                               ggCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCAGGGTGCTTGGCAAACCGCAAACGGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGATCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGtacta                                                                                                                                                                                                       <  1:223575/189‑1 (MQ=255)                                                                                               ggCAATCCAGAGACAGCGGAGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGtacta                                                                                                                                                                                                       >  2:223575/1‑189 (MQ=255)                                                                                                gCAACCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAACCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCtctc                                                                                                                                        <  2:211574/251‑1 (MQ=255)                                                                                                     ccAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATc                                                                                                                                   >  1:56045/1‑251 (MQ=255)                                                                                                           aCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                                                                                             >  2:529566/1‑251 (MQ=255)                                                                                                             aGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAg                                                                                                                           <  1:10047/251‑1 (MQ=255)                                                                                                                  cgTTATCTGGCCTTGGAGAAAGTTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGCCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGa                                                                                                                      <  2:196723/251‑1 (MQ=255)                                                                                                                         gggccctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTAtt                                                                                                               <  2:451623/246‑1 (MQ=255)                                                                                                                           gCTCTGGAGAAAGCTTATAACAGAGGAAAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGAGCAAAGACGAAGATGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTgg                                                                                                             <  1:43090/251‑1 (MQ=255)                                                                                                                                ggAGAAAGCTTATAACAGAGGATAACCCCGCATGGTGCTGGGCAAACCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGTG‑TGAAAAAACGGAAACGCTGGACTGCATCGTTAAAGGCTCGGTTGCAGTGCTGATCAAAGACGAAGATGGTAAAGAAATTATCCTCTCCAATCTGAATCAGGGTGATTTTATTGGCGAAc                                                                                                        <  1:63110/251‑1 (MQ=255)                                                                                                                                 gagaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACt                                                                                                       <  1:316443/251‑1 (MQ=255)                                                                                                                                    aaaGCTTATAACAGAGGATAAACGCGCATGGTGCTTGGCAAACCGCAAACCGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTggg                                                                                                    >  1:171992/1‑251 (MQ=255)                                                                                                                                      aGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGcc                                                                                                  >  1:183295/1‑251 (MQ=255)                                                                                                                                             aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑a                                                                                          <  2:183295/251‑1 (MQ=255)                                                                                                                                                                           aaaCCGCAAACAAACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACg                                                            <  2:171992/251‑1 (MQ=255)                                                                                                                                                                                   aaCAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACTTAGCGCATGGGTACGTGCGaaaa                                                    <  2:182444/251‑1 (MQ=255)                                                                                                                                                                                            gACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCATCCAAGAGCAAGCTTATTCCCCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACGGGCGAAAACCGCCtgtg                                           <  1:106164/251‑1 (MQ=255)                                                                                                                                                                                             aCTCTCGAATGGTTTTTGTCTCCTGGCCCGATTCATAAGTCCCCACCCAAGAGAATGCTGATTC‑CCAGGGGTGAAAACGCGGAACCGCGGTTCTAAATCGTAAAAGGCTCTGGGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAGGATCCTCTCCTATCTGAATCGGGTTGATTTTATGGGAGAAC‑GGGCCTTTTTGGAGGAGGGCAAGGAGCGTAGCGCATGGGTACGTGGAAAAACCGCCtgtg                                           <  2:167847/250‑1 (MQ=255)                                                                                                                                                                                                    aaTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTATAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCt                                   >  2:461018/1‑251 (MQ=255)                                                                                                                                                                                                          tCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAAtt                             >  1:12681/1‑251 (MQ=255)                                                                                                                                                                                                                                     aCCCATCCAAGAGCAAGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt  >  1:159542/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              GTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGG‑TGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTG‑AAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATT  >  NC_000913/3485962‑3486440

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