Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,161	G→A	9.0%	noncoding (76/171 nt) Q6* (CAA→TAA)	rttR ← tpr ←	rtT sRNA, processed from tyrT transcript protamine‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,161	0	G	A	9.0%	120.4 / 4.0	58	noncoding (76/171 nt) Q6* (CAA→TAA)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand):  ref base G (29/20);  new base A (3/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

GCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTG  >  NC_000913/1287013‑1287309                                                                                                                                                     |                                                                                                                                                     gCAGAATATTTGATTGCGGATTCGCTTGAGAGTGCAGGGACTTTTGAAAGTGAGGGGGGTGGGGGAATGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                                                                     <  1:1031584/150‑1 (MQ=255)             aTTGCGGATTCGTTTGAGACTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                                         <  2:641083/150‑1 (MQ=255)              ttGCGGATTCGGTTGAGATTTCCGGGTCTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                                        <  2:1819704/150‑1 (MQ=11)                              gaattCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                           >  1:285790/4‑131 (MQ=0)                              gaattCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                           <  2:285790/128‑1 (MQ=0)                                gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCGTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGCGCGTTGGTCGGAGCTTCTCATCCTTGCCCGccccccct                                                                                                                      >  1:211695/1‑143 (MQ=14)                                gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                      >  2:705312/1‑150 (MQ=255)                                        gcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                                         <  2:283782/122‑1 (MQ=37)                                        gcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                                         >  1:283782/2‑123 (MQ=37)                                         cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                                  <  1:1309560/145‑1 (MQ=37)                                         cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGACGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                                  >  2:1309560/1‑145 (MQ=21)                                                  gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                                                                    <  1:929024/150‑1 (MQ=38)                                                   tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCg                                                                                                       <  2:965604/146‑1 (MQ=38)                                                   tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCg                                                                                                       >  1:965604/1‑146 (MQ=21)                                                          ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTATCGGTGGCGCGCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGGATGGGGAGAATTTTTTATTGCGGGTtgggtggg                                                                                            >  1:661600/1‑150 (MQ=11)                                                                ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCCAAGCTTCTAATCCATCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGtta                                                                                      >  2:820826/1‑149 (MQ=255)                                                                    ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                                                       <  2:1297803/97‑1 (MQ=33)                                                                    ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCAt                                                                                                                                       >  1:1297803/1‑97 (MQ=18)                                                                         aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCGAGTCGAAGCTTGGTCGAAGCATCTCATCCTTCCCCGGATGGGCAGAATTTTTAATTGCGGATTGGTTGGGAAGTTCAGGGACgtt                                                                             >  1:145048/1‑150 (MQ=17)                                                                          cTCAGCGCTGCGCGCTTCGTCCTTAGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGTGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                                                            <  1:129083/150‑1 (MQ=255)                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                                                      >  2:747766/1‑92 (MQ=33)                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                                                      <  1:747766/92‑1 (MQ=33)                                                                            cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAGTGCGGAGTCGTTGGGAAGTTCAGGGACTTTTGa                                                                          >  1:323622/1‑150 (MQ=255)                                                                                gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGt                                                                      <  1:705312/150‑1 (MQ=255)                                                                                  tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTATTTGCGGATTCGTTGGGAAGTTCTGGGGCTTTTGAAAGTgt                                                                    >  1:482099/1‑149 (MQ=255)                                                                                  ggcgcgcTTGTCCGTTCGGGTGGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                    <  1:634287/149‑1 (MQ=255)                                                                                   gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAt                                                                   <  1:1825882/150‑1 (MQ=255)                                                                                   gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCAAACGCTGGCGCTCGAGTCGAACCGTGGTCGAAACTTATCAAACATACCCCGAAGTGCAGAATATTTAATTTCGGATTTGTTGGGGAATTGAGGGCCATTTGaaaattat                                                                   >  2:1516884/1‑145 (MQ=9)                                                                                        cTTCGCCCTTCGGGTCGTTGCCTGCAGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                                                              <  1:1772946/150‑1 (MQ=255)                                                                                              ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                                                                   <  1:335640/123‑1 (MQ=255)                                                                                              ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                                                                   >  2:335640/1‑123 (MQ=255)                                                                                              ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGGCGCGCGAGTCGAACCTTGGTAGAAACTTCTCAACCTTACCCCCAAGGGGAGAAAATTTTATTTCGGAATGGGTTGGAAGGTAAGGGGCATTTTAAATTTATGGTgggggggg                                                        >  2:676237/1‑150 (MQ=9)                                                                                                    gagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                  >  2:1759811/4‑150 (MQ=255)                                                                                                      gcGTTGCCGCAGGCAACGCTCTCTTGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCAACCTTCCCCGCATCGGCAGAAAATTAAAAAGCGGAAACAATGGAAAATAAAGGGACTTTTGAAAGTGATGGTGGAGGGGGAAGGaatg                                                >  2:1451225/2‑147 (MQ=9)                                                                                                            ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGTTGGTGGGGGAAGGATTCGAAcat                                          >  1:1198106/1‑148 (MQ=255)                                                                                                            ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAcgt                                          >  1:1339575/1‑148 (MQ=255)                                                                                                            ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTCGAACCt                                          >  1:1981805/1‑150 (MQ=255)                                                                                                              tGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc                                        <  2:1702699/150‑1 (MQ=255)                                                                                                               gCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCg                                       <  1:896086/150‑1 (MQ=255)                                                                                                                      cGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAgaaga                                >  2:1911484/1‑146 (MQ=255)                                                                                                                      cGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGa                                >  2:102957/1‑150 (MQ=255)                                                                                                                        ctctctCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATg                              <  1:102957/150‑1 (MQ=255)                                                                                                                         tctctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTAAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                             >  1:39683/1‑150 (MQ=255)                                                                                                                            ctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                          >  1:18316/1‑150 (MQ=255)                                                                                                                              cGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa                        >  1:1809503/1‑150 (MQ=255)                                                                                                                                      ctctAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGTAAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACa                <  1:1664996/150‑1 (MQ=255)                                                                                                                                           gccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                    >  1:355407/3‑77 (MQ=255)                                                                                                                                           gccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                    <  2:355407/75‑1 (MQ=255)                                                                                                                                            tCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGc          <  1:1516884/150‑1 (MQ=255)                                                                                                                                            ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                                  >  2:1415658/2‑78 (MQ=255)                                                                                                                                            ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                                  <  1:1415658/77‑1 (MQ=255)                                                                                                                                             cGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                                                               <  1:808547/64‑1 (MQ=40)                                                                                                                                             cGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                                                               >  2:808547/1‑64 (MQ=40)                                                                                                                                             cGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGca         >  1:675418/1‑149 (MQ=255)                                                                                                                                              gAACCTTTGTCGAAGCTTCTCTTCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGTAAGTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACATTCGAAGTCGATGACGGCAGATTTACAGTCTGCTc        <  1:977126/150‑1 (MQ=18)                                                                                                                                               aaCCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTcc       <  1:1968153/150‑1 (MQ=255)                                                                                                                                                 ccTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTCATGGTGGTGGGGGAAGGATTCGAACATTCTAAGTCGATGACCGCTGATTTACATTCTGCTCCCt     >  2:974427/1‑150 (MQ=21)                                                                                                                                                 ccTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAACTTCAGGGACTTTTGAAAGTGAAGCTGGAGGGGGAAGGATTCTACCAATCTAACTCCATGACGCCAGATTTCCAGTCTTCTCCCt     >  2:968012/1‑150 (MQ=17)                                                                                                                                                    tGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTg  >  1:451432/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTG  >  NC_000913/1287013‑1287309

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