Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,774,947	G→A	G544S (GGT→AGT)	yfjW →	uncharacterized protein YfjW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,774,947	0	G	A	91.9%	104.2 / 0.3	37	G544S (GGT→AGT)	yfjW	uncharacterized protein YfjW
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (20/14);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

CAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCC  >  NC_000913/2774710‑2775183                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             caggttAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTGTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAAt                                                                                                                                                                                                                                 <  2:491604/247‑1 (MQ=255)        gTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGCAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCtat                                                                                                                                                                                                                           >  2:379271/1‑250 (MQ=255)             aaCGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAatt                                                                                                                                                                                                                     >  1:11781/1‑251 (MQ=255)                                                             gaatggAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTg                                                                                                                                                                       <  1:400510/245‑1 (MQ=255)                                                                          aTGGACTCGTTAATAAAAGCAAAAAAATCAAAATATATGGCACCGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTc                                                                                                                                                        >  2:814830/1‑251 (MQ=255)                                                                                    taataaAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGc                                                                                                                                              >  1:192919/1‑251 (MQ=255)                                                                                             aaaaaaaaTCAAAATATATGGCACAGCAGATCTTGATTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAAGTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTGGGCTCCAATCGGAGCCTTTTTATg                                                                                                                                     <  2:658229/250‑1 (MQ=255)                                                                                                aaaaaTCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACt                                                                                                                                  >  2:654263/1‑251 (MQ=255)                                                                                                       aaaaTATATGGCACAGCCGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTACCGCATTTGGCTCCATTTGGAGCCTTTTTATGACtatttatt                                                                                                                           >  2:361174/1‑251 (MQ=255)                                                                                                                         gATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGAATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTa                                                                                                           >  2:779554/1‑249 (MQ=255)                                                                                                                          aTCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTa                                                                                                                                                                                                         <  2:814309/154‑1 (MQ=255)                                                                                                                          aTCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTa                                                                                                                                                                                                         >  1:814309/1‑154 (MQ=255)                                                                                                                              tGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCt                                                                                                      <  1:654263/249‑1 (MQ=255)                                                                                                                              tGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTAGTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGTCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGcc                                                                                                       <  2:351956/248‑1 (MQ=255)                                                                                                                                  tttGTTTATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGttttt                                                                                                >  1:556860/1‑251 (MQ=255)                                                                                                                                  tttGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGtttt                                                                                                 <  1:779554/250‑1 (MQ=255)                                                                                                                                     gttgATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTaa                                                                                              <  2:397965/250‑1 (MQ=255)                                                                                                                                                   aaTGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCg                                                                                <  1:463657/250‑1 (MQ=255)                                                                                                                                                                aataTAAATTATTCGCAATCTAAGCTAGAAATTTTTGGTTCTTATATCATGGATATATTAAAGTAAATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGTATGAAAATTTATTCTTGTAACGCATTTGACTCCAATTGGAGCCTTTTTATGACTATTTATTACAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGACGCggcaactc                                                                  >  2:806593/4‑245 (MQ=255)                                                                                                                                                                  aaTAACTTATTTGCAATCTCAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTCTTCTTGTAACGCATTTGACTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCAt                                                                >  1:448422/1‑251 (MQ=255)                                                                                                                                                                     aaCTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGTTATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGTCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGAGAAGCGGACACACATAcc                                                             >  2:436913/1‑251 (MQ=255)                                                                                                                                                                          aTTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCt                                                         <  1:479982/250‑1 (MQ=255)                                                                                                                                                                              gCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTaaa                                                    <  2:192919/251‑1 (MQ=255)                                                                                                                                                                                     aaGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACCTTAAGAATGATGTTAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTGTTTATGACTATTTATTAAAATATCCTAAGGTAAGCCTTTGTTGAAAGCAGGCTTGTGCGTGAAGCGGACACTCATACCTGCTTTAAAggg                                                 <  1:702555/247‑3 (MQ=255)                                                                                                                                                                                                    ggTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCa                              <  2:41244/251‑1 (MQ=255)                                                                                                                                                                                                      ttcttcTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCAt                             >  1:562300/1‑250 (MQ=255)                                                                                                                                                                                                      ttcttcTATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTATATAAAATTAATTAATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGAATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTACGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTACAGTGGCGTTTTTGTCGTTACCCAt                             >  2:702896/1‑250 (MQ=255)                                                                                                                                                                                                       tcttctATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGAGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATAc                           >  2:761219/1‑251 (MQ=255)                                                                                                                                                                                                        cttctATTATGGATATATTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATAcc                          >  1:435884/1‑251 (MQ=255)                                                                                                                                                                                                          tctATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGACTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATAc                           <  2:305080/248‑1 (MQ=255)                                                                                                                                                                                                                     tatatTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACATGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCtgtt             >  2:384555/1‑251 (MQ=255)                                                                                                                                                                                                                        atTAAAGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTgttgtt          >  2:268603/1‑251 (MQ=255)                                                                                                                                                                                                                            aaGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                 >  2:625956/1‑208 (MQ=255)                                                                                                                                                                                                                            aaGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                 <  1:625956/208‑1 (MQ=255)                                                                                                                                                                                                                            aaGTATATATTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAAc      >  1:817493/1‑251 (MQ=255)                                                                                                                                                                                                                               tatataTTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTATAAGAATGATGTGAATGATAATTTATTCTTGTAACGCATTTGTCTCCAATTGGAGACTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTAGAGCGTGAAGCGGACACTCATACCTGCTTTAAATTGACGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTACCTGc   >  2:184494/1‑251 (MQ=255)                                                                                                                                                                                                                                atatatTTGGTTTAAGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGcc  <  1:379271/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             CAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCC  >  NC_000913/2774710‑2775183

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