Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	216,568	C→T	A146T (GCC→ACC)	yaeF ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,568	0	C	T	95.3%	270.8 / 1.6	85	A146T (GCC→ACC)	yaeF	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (47/34);  total (49/36)
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.95e-01

CGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGA  >  NC_000913/216476‑216661                                                                                             |                                                                                              cGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgcg                                                                                        >  1:227131/1‑100 (MQ=255)   gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                      <  2:908031/100‑1 (MQ=255)   gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                      >  2:152535/1‑100 (MQ=255)       cAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGata                                                                                  >  2:1307471/1‑100 (MQ=255)         gTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATc                                                                                >  1:1954108/1‑100 (MQ=255)         gTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATc                                                                                >  2:1615538/1‑100 (MQ=255)            aGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgt                                                                             <  1:860678/100‑1 (MQ=255)             gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                            <  2:1820993/100‑1 (MQ=255)             gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                            <  2:479732/100‑1 (MQ=255)                caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                         >  1:1105171/1‑100 (MQ=255)                caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                         >  2:1761740/1‑100 (MQ=255)                caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                         >  2:685453/1‑100 (MQ=255)                caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGAACTTCTTCGCCAACACggccatcaat                                                                               >  1:1396268/1‑86 (MQ=255)                   cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                      >  2:1282209/1‑100 (MQ=255)                   cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                      >  1:5581/1‑100 (MQ=255)                     ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                    >  1:516521/1‑100 (MQ=255)                     ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                    >  2:1921873/1‑100 (MQ=255)                         aGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCggg                                                                >  1:753062/1‑100 (MQ=255)                          gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                               >  2:589121/1‑100 (MQ=255)                          gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                               <  2:646713/100‑1 (MQ=255)                           tCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGt                                                              <  1:1795269/100‑1 (MQ=255)                            ctccATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                             <  2:1284617/98‑1 (MQ=255)                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                             <  1:152536/100‑1 (MQ=255)                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATCTTATTGGTAAACGCGGTGATCTCTGTCGCTTGTTGCGGGGTa                                                             >  2:1916223/1‑100 (MQ=255)                             aCCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTaa                                                            <  2:1176715/100‑1 (MQ=255)                              ccATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAg                                                           >  1:1284623/1‑100 (MQ=255)                              ccATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAg                                                           >  1:304845/1‑100 (MQ=255)                               cATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGa                                                          >  1:95057/1‑100 (MQ=255)                                aTAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGCGATATCTGTGGCTTGTTGCGGGGTCAGa                                                          >  2:1904890/1‑99 (MQ=255)                                  aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                       >  1:529510/1‑100 (MQ=255)                                  aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                       <  1:564975/100‑1 (MQ=255)                                  aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                       >  1:1254826/1‑100 (MQ=255)                                  aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                       >  1:251329/1‑100 (MQ=255)                                   aaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                      >  2:717446/1‑100 (MQ=255)                                   aaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                      >  1:908047/1‑100 (MQ=255)                                     gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                    >  1:2022376/1‑100 (MQ=255)                                         aTAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTc                                                >  1:1101774/1‑100 (MQ=255)                                          tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                               <  1:2015846/100‑1 (MQ=255)                                           aaaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                              >  1:1025529/1‑100 (MQ=255)                                           aaaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                              >  1:720949/1‑100 (MQ=255)                                            aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                             <  1:876035/100‑1 (MQ=255)                                            aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                             <  2:1101769/100‑1 (MQ=255)                                             aTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTaa                                            >  1:1108537/1‑100 (MQ=255)                                              ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                           <  1:1921891/100‑1 (MQ=255)                                              ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                           <  2:1363803/100‑1 (MQ=255)                                                cGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGAc                                         <  2:1121739/100‑1 (MQ=255)                                                cGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGAc                                         >  2:1679105/1‑100 (MQ=255)                                                cGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGAc                                         <  2:235896/100‑1 (MQ=255)                                                  gatctGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGa                                       <  2:882482/96‑1 (MQ=255)                                                   cAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaa                                      >  1:1744172/1‑100 (MQ=255)                                                    aaTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaa                                     >  1:1970646/1‑100 (MQ=255)                                                     aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTAGTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaa                                     <  2:503327/99‑1 (MQ=255)                                                       gCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaagatc                                  >  1:882485/1‑96 (MQ=255)                                                       gCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGc                                  <  1:1307479/100‑1 (MQ=255)                                                       gCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGc                                  >  2:1832463/1‑100 (MQ=255)                                                        ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                                 <  2:737368/100‑1 (MQ=255)                                                        ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                                 >  2:1503624/1‑100 (MQ=255)                                                            gATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATc                             >  2:188697/1‑100 (MQ=255)                                                             aTAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCa                            <  1:1530140/100‑1 (MQ=255)                                                             aTAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCa                            >  2:1587259/1‑100 (MQ=255)                                                              tAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCAc                           >  2:129977/1‑100 (MQ=255)                                                                gTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTa                         >  2:962979/1‑100 (MQ=255)                                                                gTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTa                         >  1:852169/1‑100 (MQ=255)                                                                gTTATAAACGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTa                         >  1:1080297/1‑100 (MQ=255)                                                                 ttATAACCGCTATCTTTGATTTTATTGGTCAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTAt                        <  2:1232464/100‑1 (MQ=255)                                                                 ttATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTAt                        <  2:5581/100‑1 (MQ=255)                                                                    taACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                     <  2:2023026/100‑1 (MQ=255)                                                                      aCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTc                   >  1:392634/1‑100 (MQ=255)                                                                      aCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGAGCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTc                   <  2:1254820/100‑1 (MQ=255)                                                                       ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                  >  1:1257616/1‑100 (MQ=255)                                                                       ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATACCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                  <  2:1253806/100‑1 (MQ=255)                                                                        cGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCAt                 >  1:1119170/1‑100 (MQ=255)                                                                        cGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCAt                 <  2:1642307/100‑1 (MQ=255)                                                                          cTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCg               <  1:717448/100‑1 (MQ=255)                                                                          cTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCg               <  2:2022358/100‑1 (MQ=255)                                                                              cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt           >  2:150543/1‑100 (MQ=255)                                                                                ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt         >  2:311580/1‑100 (MQ=255)                                                                                ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt         >  1:574550/1‑100 (MQ=255)                                                                                  gATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTa       >  2:561944/1‑100 (MQ=255)                                                                                     tttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAgg    <  2:852166/100‑1 (MQ=255)                                                                                     tttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAgg    <  2:1970628/100‑1 (MQ=255)                                                                                     tttATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAgg    <  1:974926/100‑1 (MQ=255)                                                                                      ttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGTGTAAGATCCGGGACTCGTATGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAggg   <  2:1119165/100‑1 (MQ=255)                                                                                       tATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGa  <  1:1904908/100‑1 (MQ=255)                                                                                       tATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGa  <  1:1679118/100‑1 (MQ=255)                                                                                             |                                                                                              CGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGA  >  NC_000913/216476‑216661

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