Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,428,700	A→G	D9D (GAT→GAC)	insH‑5 ←	Rac prophage; IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,428,700	0	A	G	100.0%	240.3 / NA	62	D9D (GAT→GAC)	insH‑5	Rac prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/34);  total (28/34)

AGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACG  >  NC_000913/1428564‑1428835                                                                                                                                         |                                                                                                                                        aGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCg                                                                                                                                      >  2:661468/1‑140 (MQ=255)    gacGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCg                                                                                                                                   <  2:939252/138‑1 (MQ=255)    gCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCg                                                                                                                                   <  1:757782/140‑1 (MQ=255)        tCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAgg                                                                                                                               >  2:1073636/1‑140 (MQ=255)        ccGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCg                                                                                                                                   <  2:606647/135‑1 (MQ=14)        ccGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCg                                                                                                                                   >  1:606647/2‑136 (MQ=14)         cGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCg                                                                                                                                   <  2:995481/135‑1 (MQ=255)         cGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCg                                                                                                                                   >  1:995481/1‑135 (MQ=255)                tGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGa                                                                                                                       <  1:312418/140‑1 (MQ=255)                     gCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGACt                                                                                                                  <  1:405486/140‑1 (MQ=255)                       cTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCa                                                                                                                <  2:469263/140‑1 (MQ=255)                       cTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCa                                                                                                                <  2:198183/140‑1 (MQ=255)                           ggATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatgat                                                                                                            <  1:618673/140‑1 (MQ=255)                            gATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                                                           <  1:1085546/140‑1 (MQ=255)                            gATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                                                           <  2:861402/140‑1 (MQ=255)                            gATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                                                           <  1:159323/140‑1 (MQ=255)                                aaaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTccc                                                                                                       <  2:167399/140‑1 (MQ=255)                                 aaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCct                                                                                                      >  1:230031/1‑140 (MQ=255)                                 aaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCct                                                                                                      >  2:237709/1‑140 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                                                           <  1:219924/130‑1 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                                                           >  2:219924/1‑130 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatga                                                                                                             <  2:786181/128‑1 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCatga                                                                                                             >  1:786181/1‑128 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTc                                                                                                         >  2:595185/1‑132 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTc                                                                                                         <  1:595185/132‑1 (MQ=255)                                      cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTggg                                                                                                 >  2:112394/1‑140 (MQ=255)                                           tGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATgcg                                                                                            <  2:95887/140‑1 (MQ=255)                                               aGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTcc                                                                                        >  2:61883/1‑140 (MQ=255)                                                  ggTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGa                                                                                     <  1:658655/140‑1 (MQ=255)                                                  ggTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGa                                                                                     <  1:839335/140‑1 (MQ=255)                                                   gTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGAt                                                                                    <  1:237709/140‑1 (MQ=255)                                                    tCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATg                                                                                   >  2:382780/1‑140 (MQ=255)                                                         ttttGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAatat                                                                              >  1:314764/4‑140 (MQ=14)                                                           cTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATatga                                                                            >  1:894778/1‑140 (MQ=255)                                                           cTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATatga                                                                            <  1:382780/140‑1 (MQ=255)                                                            tGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATatgat                                                                           <  1:661468/140‑1 (MQ=255)                                                              ccATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATAtgatga                                                                         >  1:707472/1‑140 (MQ=255)                                                              ccATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATAtgatga                                                                         >  2:126020/1‑140 (MQ=255)                                                                 tGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGAtct                                                                      >  1:1121921/1‑140 (MQ=255)                                                                       aaTCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATc                                                                <  1:61883/140‑1 (MQ=255)                                                                            gCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGaa                                                           <  2:34106/140‑1 (MQ=255)                                                                             cTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                                                          >  1:907906/1‑140 (MQ=255)                                                                             cTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                                                          >  1:234361/1‑140 (MQ=255)                                                                             cTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                                                          >  1:722757/1‑140 (MQ=255)                                                                             cTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                                                          >  2:990114/1‑140 (MQ=255)                                                                             cTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                                                          >  1:579339/1‑140 (MQ=255)                                                                             cTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                                                          <  1:1073636/140‑1 (MQ=255)                                                                                 aTGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGt                                                      <  2:907906/140‑1 (MQ=255)                                                                                 aTGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGt                                                      >  2:126420/1‑140 (MQ=255)                                                                                  tGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGtt                                                     <  2:230031/140‑1 (MQ=255)                                                                                  tGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGtt                                                     >  1:580074/1‑140 (MQ=255)                                                                                  tGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGtt                                                     >  1:427274/1‑140 (MQ=255)                                                                                      ggAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCa                                                 <  1:195048/140‑1 (MQ=255)                                                                                        agagGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCAcc                                               <  1:990114/140‑1 (MQ=255)                                                                                            gAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGtt                                                     >  2:383053/1‑130 (MQ=255)                                                                                            gAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGtt                                                     <  1:383053/130‑1 (MQ=255)                                                                                              aaaTCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCa                                         <  2:894778/140‑1 (MQ=255)                                                                                                 tctcTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAgg                                      <  2:234361/140‑1 (MQ=255)                                                                                                     ttttCGGGTCTGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGaa                                  <  1:126420/140‑1 (MQ=255)                                                                                                               tGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCacgac                        >  2:219850/1‑140 (MQ=255)                                                                                                                    gcgcTTAGTGCTGAATTCACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCt                   >  1:762541/1‑140 (MQ=255)                                                                                                                                     cACTGTCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACg  <  1:178230/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        AGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACG  >  NC_000913/1428564‑1428835

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