Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	19.8%	S27S (AGT→AGC)	tpn ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,471	0	A	G	19.8%	92.5 / 16.0	64	S27S (AGT→AGC)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (9/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCG  >  CP000730/36332‑36608                                                                                                                                            |                                                                                                                                          gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCCGATATTTCGCGATaa                                                                                                                                           >  1:296488/1‑140 (MQ=11) gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                           >  2:112231/1‑140 (MQ=25)  ccttttGTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                          <  2:164249/137‑1 (MQ=255)  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAAc                                                                                                                                          >  2:9725/1‑140 (MQ=255)  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                          >  2:253726/1‑140 (MQ=25)  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGAGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGACAAc                                                                                                                                          >  2:111098/1‑140 (MQ=255)  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAAc                                                                                                                                          >  2:90814/1‑140 (MQ=255)    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTCTATACAGTCCTTTTTCACCTCGTTCCATTAATATGTCAAATATATCACGATAACTc                                                                                                                                        >  1:231109/1‑140 (MQ=11)     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                       <  2:279583/140‑1 (MQ=18)     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                       <  1:69881/140‑1 (MQ=21)     tttGTAATAAGCTTTTTTAGGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATGCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                       >  1:192272/1‑140 (MQ=11)      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  1:96719/140‑1 (MQ=11)      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCaa                                                                                                                                      >  1:302352/1‑140 (MQ=21)      ttGGAATCAGCTTTTTTATGCTTTTTCTTCCAACTTTGATGCAAAATGTGGCCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:66765/140‑1 (MQ=11)       tGTACTAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCCCAAAAGCTCAAt                                                                                                                                     >  2:276793/1‑140 (MQ=11)        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATCAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    <  2:306129/140‑1 (MQ=21)        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCGACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  2:35310/1‑140 (MQ=11)        gTAATAAGCTTTTTTATGACTATTATACCAAATACGAGGTAAAATTTGGGCATATTGTTGACCCCAACAGTAGGGCTGTTAAGGGGCGGCCTTGACACCACTTTCCCGTTATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    <  2:27888/140‑1 (MQ=2)         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAACCCAACGGTAAACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAGCTATATCATGATAACTCaacac                                                                                                                                   >  2:50910/1‑137 (MQ=11)         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGTTATATCACCATAGCTCAATgc                                                                                                                                   >  2:41429/1‑139 (MQ=11)         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATGTCACGATAACTCAATgc                                                                                                                                   >  1:262844/1‑139 (MQ=16)         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCACTgc                                                                                                                                   >  1:229535/1‑139 (MQ=16)          aataaGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGAATATTCTTGAACCCAACAGTCGTTCGTTGAATAATGAACGTTTACACCACGTTCTATTAATATATCACATATATCAAGATAGCGCAAtcca                                                                                                                                  >  2:295049/1‑137 (MQ=11)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCGTATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                  <  2:218039/140‑3 (MQ=12)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGACCGTTTACACCACGTTCTCTTAATATTTCCGATATATCACGATAGCTCAATgca                                                                                                                                  >  1:46820/1‑138 (MQ=11)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  >  2:125489/1‑138 (MQ=11)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  <  2:279768/140‑3 (MQ=11)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACCATAGCTCAATgca                                                                                                                                  >  1:13704/1‑138 (MQ=11)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTGTACTCCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  <  2:176336/140‑3 (MQ=11)          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                  <  2:244742/140‑3 (MQ=17)          aataaGCTTTTTTATGCTTTTTCTTCCAAATGTGATATAAAATTGGAGCATAGTCTTCACAGCAACGATTGACCGTTGAATGATGAACGTTTACACCACGATCCCTTAATATTTCAAATATATCACAATAACTCAAtaca                                                                                                                                  >  2:286840/1‑137 (MQ=11)           ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCACTgctt                                                                                                                                 >  2:87544/1‑137 (MQ=255)           ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATACCTCAATGtat                                                                                                                                 >  1:60004/1‑140 (MQ=255)             aaGCTTTTTTATGCTTTTTCTTCCAACTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTatat                                                                                                                               >  1:306940/1‑140 (MQ=255)                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAACGTATATCTTAGAtagt                                                                                                                     <  1:158796/140‑1 (MQ=25)                        tGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                    <  2:295903/140‑1 (MQ=35)                        tGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                    <  1:199742/140‑1 (MQ=35)                         gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAACTTAGAtagtag                                                                                                                   >  1:4352/1‑140 (MQ=25)                         gCTTGTTCTTCCCATTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagtag                                                                                                                   <  2:133656/140‑1 (MQ=255)                              ttcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAc                                                                                                              >  2:12640/1‑140 (MQ=35)                                  tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATAAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAATGAATATCTTAGATAGTAGCCAACGGCt                                                                                                          >  2:66217/1‑140 (MQ=12)                                         ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAgaga                                                                                                   >  1:3088/1‑140 (MQ=25)                                         ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGa                                                                                                   <  2:186893/140‑1 (MQ=35)                                                agaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATc                                                                                            <  1:144650/138‑1 (MQ=18)                                                     tGGGGGATATTCTGGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                       <  2:23722/140‑1 (MQ=17)                                                        ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                    <  2:294509/140‑1 (MQ=32)                                                        ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                    >  2:61128/1‑140 (MQ=32)                                                        ggCATATTCTTAAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                    <  1:91274/140‑1 (MQ=18)                                                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                        >  1:292517/1‑140 (MQ=32)                                                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                        >  1:245500/1‑140 (MQ=33)                                                                       ccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                     <  2:275889/140‑1 (MQ=33)                                                                          aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAAt                                                                  <  2:295702/140‑1 (MQ=33)                                                                          aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAATTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAAt                                                                  >  2:205027/1‑140 (MQ=18)                                                                                              gaACGTTTACACCACGTTCTGTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAGTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCt                                              >  1:207628/1‑140 (MQ=255)                                                                                              gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCt                                              >  2:275793/1‑140 (MQ=32)                                                                                                  gTTTACACCTCGTTTCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                          <  2:135760/140‑1 (MQ=17)                                                                                                     tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                       >  1:149302/1‑140 (MQ=32)                                                                                                     tACAACACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                       <  2:86624/140‑1 (MQ=17)                                                                                                       caccacGTTCCCTTAATCTTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATACTTGTTAAATTGTTTATATCTGAAATAGCTGCTACAGAAGACTCCTTTTTGTTaa                                     >  1:296678/1‑140 (MQ=11)                                                                                                              ttCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtata                              <  2:4352/140‑1 (MQ=32)                                                                                                               tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                             <  1:132764/140‑1 (MQ=33)                                                                                                                      atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                      <  2:75347/140‑1 (MQ=35)                                                                                                                      atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                      <  1:12640/140‑1 (MQ=34)                                                                                                                                     cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg       <  1:207486/140‑1 (MQ=35)                                                                                                                                      aCGATAGCTCAATGCATTTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGa      <  1:9725/140‑1 (MQ=255)                                                                                                                                          tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATGAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCg  <  1:61128/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          GCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCG  >  CP000730/36332‑36608

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