Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,446,372	T→C	E58E (GAA→GAG)	rpsH ←	30S ribosomal subunit protein S8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,446,372	0	T	C	100.0%	219.2 / NA	61	E58E (GAA→GAG)	rpsH	30S ribosomal subunit protein S8
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/22);  total (39/22)

AACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACGGATA  >  NC_000913/3446237‑3446507                                                                                                                                        |                                                                                                                                        aacTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAg                                                                                                                                     >  2:1083886/1‑140 (MQ=255)    tGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCt                                                                                                                                  >  2:824787/1‑140 (MQ=255)        aTACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGgt                                                                                                                              >  1:280411/1‑140 (MQ=255)              aGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGgtg                                                                                                                             >  2:28254/1‑135 (MQ=255)              aGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGgtg                                                                                                                             <  1:28254/135‑1 (MQ=255)                ccccGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCtt                                                                                                                      <  1:674467/139‑1 (MQ=255)                aCCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCtt                                                                                                                      >  1:588070/1‑140 (MQ=255)                     ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGgtg                                                                                                                             <  2:430675/128‑1 (MQ=255)                     ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGgtg                                                                                                                             >  1:430675/1‑128 (MQ=255)                     ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACt                                                                                                                 >  2:117116/1‑140 (MQ=255)                     ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACt                                                                                                                 >  2:531393/1‑140 (MQ=255)                     ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACt                                                                                                                 >  1:165694/1‑140 (MQ=255)                              tCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCt                                                                                                        >  2:1165576/1‑140 (MQ=255)                              tCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGACTTCAACTTTAAAATCt                                                                                                        <  1:1131779/140‑1 (MQ=255)                                                    ataGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCa                                                                                  >  2:138998/1‑140 (MQ=255)                                                     taGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAg                                                                                 <  1:824787/140‑1 (MQ=255)                                                        aTGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCAc                                                                              >  1:636705/1‑140 (MQ=255)                                                          gcgcAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAAc                                                                                              <  1:647001/122‑1 (MQ=255)                                                          gcgcAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAAc                                                                                              >  2:647001/1‑122 (MQ=255)                                                          gcgcAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGt                                                                            >  1:191980/1‑140 (MQ=255)                                                                    gggCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGc                                                                  <  1:433373/140‑1 (MQ=255)                                                                     ggcggcTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGcc                                                                 >  2:375534/1‑140 (MQ=255)                                                                          cTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACttt                                                            >  1:265304/1‑140 (MQ=255)                                                                          cTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACttt                                                            >  2:189086/1‑140 (MQ=255)                                                                          cTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACttt                                                            >  1:1097314/1‑140 (MQ=255)                                                                          cTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACttt                                                            >  2:587159/1‑140 (MQ=255)                                                                          cTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACttt                                                            <  2:636705/140‑1 (MQ=255)                                                                          cTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACttt                                                            >  2:1142722/1‑140 (MQ=255)                                                                               acGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCt                                                       <  1:117116/140‑1 (MQ=255)                                                                                   tGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTgga                                                   <  2:1097314/140‑1 (MQ=255)                                                                                   tGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTgga                                                   >  1:541755/1‑140 (MQ=255)                                                                                       tGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGATGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGaa                                               >  2:809013/1‑140 (MQ=255)                                                                                               acaacaGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGt                                       >  2:365263/1‑140 (MQ=255)                                                                                               acaacaGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGt                                       >  2:1093355/1‑140 (MQ=255)                                                                                               acaacaGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGt                                       <  2:588070/140‑1 (MQ=255)                                                                                                   caGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGAcc                                   <  1:1093355/140‑1 (MQ=255)                                                                                                      cTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCGGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGt                                       >  2:204783/1‑133 (MQ=255)                                                                                                      cTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCGGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGt                                       <  1:204783/133‑1 (MQ=255)                                                                                                      cTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGTTGACCGCa                                >  2:425770/1‑140 (MQ=255)                                                                                                      cTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCa                                >  2:964297/1‑140 (MQ=255)                                                                                                      cTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCa                                >  2:121873/1‑140 (MQ=255)                                                                                                         tGCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCt                             >  2:427646/1‑140 (MQ=255)                                                                                                          gCCCTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCtt                            >  2:1077132/1‑140 (MQ=255)                                                                                                           cccTGGAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCttt                           >  1:947576/1‑140 (MQ=255)                                                                                                                gAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTcg                      >  1:802829/1‑140 (MQ=255)                                                                                                                gAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTcg                      >  1:352600/1‑140 (MQ=255)                                                                                                                gAAATACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTcg                      <  1:1142722/140‑1 (MQ=255)                                                                                                                    tACTTCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGc                  <  1:587159/140‑1 (MQ=255)                                                                                                                       ttCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTg               <  2:165694/140‑1 (MQ=255)                                                                                                                        tCAGAGTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGa              <  1:541896/140‑1 (MQ=255)                                                                                                                          agagTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGAcc            >  1:1093493/1‑140 (MQ=255)                                                                                                                          agagTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGAcc            <  1:365263/140‑1 (MQ=255)                                                                                                                             gTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGATCGCGGCCTGACCGtt         <  1:427646/140‑1 (MQ=255)                                                                                                                             gTAAGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAACTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGtt         <  2:265304/140‑1 (MQ=255)                                                                                                                                aGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACg      >  1:93461/1‑140 (MQ=255)                                                                                                                                aGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACg      >  2:1143284/1‑140 (MQ=255)                                                                                                                                aGTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACg      >  1:1113189/1‑140 (MQ=255)                                                                                                                                 gTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGAcc            <  1:1048207/133‑1 (MQ=255)                                                                                                                                 gTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGAcc            >  2:1048207/1‑133 (MQ=255)                                                                                                                                 gTTCCAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACgg     >  1:922844/1‑140 (MQ=255)                                                                                                                                    ccAGCTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACGGATa  <  1:189086/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        AACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTGCCACTTTCAGCTTGGAGGAAGGCATGGTGACCGCAGCTTTGTTCGCGGCCTGACCGTTACGGATA  >  NC_000913/3446237‑3446507

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