Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,161,737	C→A	intergenic (‑186/‑148)	fhuE ← / → hinT	ferric coprogen/ferric rhodotorulic acid outer membrane transporter/purine nucleoside phosphoramidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,161,737	0	C	A	96.8%	211.4 / 0.2	62	intergenic (‑186/‑148)	fhuE/hinT	ferric coprogen/ferric rhodotorulic acid outer membrane transporter/purine nucleoside phosphoramidase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (31/29);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGAGCGATTATATTTAGTGTGCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGCGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGGCGCAACGCATGGCTATTTGAAAA  >  NC_000913/1161603‑1161873                                                                                                                                       |                                                                                                                                         tgaGCGATTATATTTAGTGTGCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGatat                                                                                                                                     >  1:784973/1‑140 (MQ=255)                  tgtgCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGa                                                                                                                    >  2:505876/1‑140 (MQ=255)                  tgtgCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGa                                                                                                                    >  2:1473847/1‑140 (MQ=255)                  tgtgCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGa                                                                                                                    >  2:1083125/1‑140 (MQ=255)                  tgtgCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGa                                                                                                                    >  1:751624/1‑140 (MQ=255)                  tgtgCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGa                                                                                                                    >  1:480826/1‑140 (MQ=255)                  tgtgCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGa                                                                                                                    >  1:394220/1‑140 (MQ=255)                     gCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGAttt                                                                                                                 >  2:256706/1‑140 (MQ=255)                     gCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGAttt                                                                                                                 >  2:790293/1‑140 (MQ=255)                     gCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGAttt                                                                                                                 >  1:472544/1‑140 (MQ=255)                      cGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGAtttt                                                                                                                >  1:1394667/1‑140 (MQ=255)                       gAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGAttttt                                                                                                               <  2:529729/140‑1 (MQ=255)                                    tgtAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATc                                                                                                  <  1:573459/140‑1 (MQ=255)                                       aaaTGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTg                                                                                               >  1:853211/1‑140 (MQ=255)                                           gCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGtt                                                                                           <  1:100317/140‑1 (MQ=255)                                           gCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGtt                                                                                           <  2:1394667/140‑1 (MQ=255)                                                tGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGt                                                                                      <  1:834162/140‑1 (MQ=255)                                                        ataCGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAagag                                                                              <  2:597322/140‑1 (MQ=255)                                                        ataCGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAagag                                                                              <  1:504022/140‑1 (MQ=255)                                                        ataCGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAagag                                                                              <  2:805128/140‑1 (MQ=255)                                                              aTTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGtt                                                                        >  1:47153/1‑140 (MQ=255)                                                                 cATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCGGGTTTaa                                                                     >  2:1248479/1‑140 (MQ=255)                                                                 cATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTaa                                                                     >  2:1277867/1‑140 (MQ=255)                                                                 cATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTaa                                                                     >  1:861920/1‑140 (MQ=255)                                                                   tatTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCGGGTTTAACt                                                                   <  1:1248479/140‑1 (MQ=255)                                                                   tatTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACt                                                                   <  2:780166/140‑1 (MQ=255)                                                                   tatTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACt                                                                   <  1:1277867/140‑1 (MQ=255)                                                                   tatTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGAAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACt                                                                   <  1:1344779/140‑1 (MQ=255)                                                                      ttGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGa                                                                <  1:997176/140‑1 (MQ=255)                                                                      ttGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGa                                                                <  2:853211/140‑1 (MQ=255)                                                                       tGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAc                                                               <  1:1129912/140‑1 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              <  2:1261682/140‑1 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              <  2:960204/140‑1 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              <  1:1152107/140‑1 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              <  2:751624/140‑1 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              >  1:960204/1‑140 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              <  1:898300/140‑1 (MQ=255)                                                                        gTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGAcc                                                              >  1:720127/1‑140 (MQ=255)                                                                            tGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATAc                                                          >  1:564022/1‑140 (MQ=255)                                                                            tGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATAc                                                          >  2:1468448/1‑140 (MQ=255)                                                                              tAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATActc                                                        <  2:1266949/140‑1 (MQ=255)                                                                              tAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATActc                                                        <  1:316983/140‑1 (MQ=255)                                                                              tAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATActc                                                        >  2:316983/1‑140 (MQ=255)                                                                                aaCCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCAt                                                            >  2:547050/1‑134 (MQ=255)                                                                                aaCCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCAt                                                            <  1:547050/134‑1 (MQ=255)                                                                                      aCAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGTAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCg                                                >  2:1390389/1‑140 (MQ=255)                                                                                       cAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGc                                               >  1:1092621/1‑140 (MQ=255)                                                                                            aaTGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGCGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAgtgt                                          >  1:1405998/1‑140 (MQ=255)                                                                                            aaTGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAgtgt                                          >  2:96183/1‑140 (MQ=255)                                                                                                tCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGgcg                                      <  1:1473847/140‑1 (MQ=255)                                                                                                tCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGgcg                                      <  1:1083125/140‑1 (MQ=255)                                                                                                  ttttCAGCGCATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGcgcg                                    >  1:1269679/1‑140 (MQ=255)                                                                                                        gcgcATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAg                              >  1:139968/1‑140 (MQ=255)                                                                                                          gcATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCa                            <  2:480826/140‑1 (MQ=255)                                                                                                          gcATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCa                            >  2:1369623/1‑140 (MQ=255)                                                                                                          gcATTCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCa                            <  1:1369623/140‑1 (MQ=255)                                                                                                              tCGCAGGCAGAAATGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGgc                        >  2:1046226/1‑140 (MQ=255)                                                                                                               cGCAGGCAGAAATGGGAATAATGGCGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGgcg                       <  2:1405998/140‑1 (MQ=255)                                                                                                                         aaTGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGGCGCAACGCATgg             <  2:720127/140‑1 (MQ=255)                                                                                                                          aTGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGgcg                       <  2:28381/129‑1 (MQ=255)                                                                                                                          aTGGGAATAATGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGgcg                       >  1:28381/1‑129 (MQ=255)                                                                                                                                    tGGAGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGGCGCAACGCATGGCTATTTGaaaa  >  1:1249742/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         TGAGCGATTATATTTAGTGTGCGAATAATTTTGTTTGTAAATGCAAATGAGAAATATACGCATTCATATTTGTCATGTAAACCAAACAGAGAATGTCTTTTCAGCGCATTCGCAGGCAGAAATGGGAATAATGGCGATATATACGGCAGCAAAACGATTTTTTGCTTAAGAAATCGTGAGTTAAGGTTGAAAGAGCAGGTTTAACTCGACCATACTCTATACTCGCAGTGTGGCGCGGCGTAGCATGGCGCAACGCATGGCTATTTGAAAA  >  NC_000913/1161603‑1161873

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