Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,968,683	G→A	M256I (ATG→ATA)	rfe →	UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,968,683	0	G	A	97.1%	104.9 / ‑6.4	35	M256I (ATG→ATA)	rfe	UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/11);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATT  >  NC_000913/3968438‑3968929                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       aTGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGAGCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGatat                                                                                                                                                                                                                                                   <  1:196244/251‑1 (MQ=255)    aTTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATAtggg                                                                                                                                                                                                                                                >  1:223333/1‑250 (MQ=255)      tttGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATAtggtgg                                                                                                                                                                                                                                              >  1:163831/1‑251 (MQ=255)                        aaCCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACAGTCGcc                                                                                                                                                                                                                            >  1:99501/1‑251 (MQ=255)                               gagctaaTCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTaa                                                                                                                                                                                                                     <  2:194855/246‑1 (MQ=255)                                    aaTCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTaa                                                                                                                                                                                                                     >  1:194855/1‑246 (MQ=255)                                              tttGCGCTGATCGCCGCCATCCTGCCATACATCATGCTCAACCGTGGTATCCTGGGTCGCCGCTAAAACGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACAGCTTTGTGGATAATAGCCATTCCGCTAATAGATCTGGTGGCGATTATGAACCGTCGCCCGCGTCAAGGCATGAGACCAtt                                                                                                                                                                                                      >  2:41963/1‑251 (MQ=255)                                                 gCGATGATCGCCGCCATCCGGCCATACATCATGCTTAACCTTGGTATCCTGGGGCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTTATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCTTtctc                                                                                                                                                                                                   <  1:185876/251‑1 (MQ=255)                                                      gatCGCCGCCATCCTGCCATACATCATGCTTCACCTTGGTATCCTGGGTCGCCGCTACAAAGACTTTATAGGTGATGCGGGCAGTACGCTGATTGGATTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGAATACTCGCCATTCCTCTAATAGATCTGGTGGAGCGTATGTACCGTCGCCTGCATAAAGGCATGAGCCCATTCACTCCTg                                                                                                                                                                                              >  2:23229/1‑251 (MQ=255)                                                            cgccATCCTGACATACATCATGCTTAAACTTGGTATCCTGGGACGCCGCTACAAAATCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGCTCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCCTCAGCCCGGTTACCGCTTTGTGGCTAATCGCCATTCCGCTAATAGATCTGGTGGCGTTTATGTACCGTCGCCTGCGTAAAGGCCTAACCCCATTCTCTCCTGACCGTc                                                                                                                                                                                        >  2:227919/1‑251 (MQ=255)                                                                              catGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTgt                                                                                                                                                                      >  1:36136/1‑250 (MQ=255)                                                                                tGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATc                                                                                                                                                                    >  1:240456/1‑251 (MQ=255)                                                                                    tAACCTGGGTATCCTGGGTCGCCGCTTCAAAATCTTTATGGGTGATGCGGGCAGTACGCTGATGGGTTTTACCGTTATCTGGATCCTGCTCGAAACTACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTTGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCGCTCCTGACCGTCAGCATATTCACCATTTGATCAt                                                                                                                                                                  <  2:223333/249‑1 (MQ=255)                                                                                            gTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCAATTCTCTCCTGACCGTCAGCCAATTCACCATTTGCTCATGCGTGCCggg                                                                                                                                                        >  2:290876/1‑251 (MQ=255)                                                                                              aTCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGtt                                                                                                                                                      <  1:195957/251‑1 (MQ=255)                                                                                                 cTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGCTAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGtttcc                                                                                                                                                   >  2:1490/1‑249 (MQ=255)                                                                                                                aaaGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTg                                                                                                                                                                                              <  1:150794/193‑1 (MQ=255)                                                                                                                aaaGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGCTATCTGGATCCTGCTCGAAACGACCCAGGGCACAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAATGGCATGAGCCCATTCTCTCCTg                                                                                                                                                                                              >  2:150794/1‑193 (MQ=255)                                                                                                                               gATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGACACGACCCAGGGCAAAACCCATCCAATCAGCCCGGTTACCGCTTTGTGGATAATCGCCCTTCCGCTAATAGATATGGTGGCGATTATGTAACGTCGCCTGCGTAAAGGCATGCGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTACCGTCAGGCATTTGTGcgtggctcgg                                                                                                                      >  2:233277/1‑241 (MQ=255)                                                                                                                               gATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCGCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTg                                                                                                                             <  1:233277/243‑1 (MQ=255)                                                                                                                                aTGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCtt                                                                                                                    <  1:250710/251‑1 (MQ=255)                                                                                                                                                             gTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGt                                                                                       >  1:226270/1‑251 (MQ=255)                                                                                                                                                                     gATCCTGCTCGAAACGACCCAGGGCAAATCCCAGCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGAGTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAg                                                                               <  2:103209/251‑1 (MQ=255)                                                                                                                                                                      aTCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGa                                                                              <  1:213358/251‑1 (MQ=255)                                                                                                                                                                                                aaCCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTggg                                                    >  1:325834/1‑251 (MQ=255)                                                                                                                                                                                                  cccatcccatCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTTGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTc                                                  >  1:84424/1‑251 (MQ=255)                                                                                                                                                                                                    catcccatcAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCAt                                                >  1:187385/1‑251 (MQ=255)                                                                                                                                                                                                     atcccatcAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTAAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGTCGTGCTGGCAGAATATTCTCTTTTTGTCCCGGGGTGGGTAATg                                               >  2:45770/1‑251 (MQ=255)                                                                                                                                                                                                       cccatcAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGATTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGcg                                             >  2:22904/1‑250 (MQ=255)                                                                                                                                                                                                         catcAGACCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCCGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGGGTGGGTCATGCTgg                                           >  2:126820/1‑251 (MQ=255)                                                                                                                                                                                                                            ttGTGGATAATCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGCGTGGGTCATGCTGGTGCTCTTTTTGCTAGCAtt                        >  1:278688/1‑251 (MQ=255)                                                                                                                                                                                                                                 gATAATCGCCATGCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATGCTTcc                   <  1:134308/251‑1 (MQ=255)                                                                                                                                                                                                                                     aTCGCCATTCCGCTAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCctct                >  1:162184/1‑250 (MQ=255)                                                                                                                                                                                                                                       cGCCATTCCGCTAATAGATATGGTGGCGATTATTTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCAGTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTAt              <  2:163831/250‑1 (MQ=255)                                                                                                                                                                                                                                                  tAATAGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCAtt  >  1:299247/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ATGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCCATTCCGCTAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATT  >  NC_000913/3968438‑3968929

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