Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,349,856	T→C	E84E (GAA→GAG)	yciW ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,349,856	0	T	C	98.2%	209.5 / NA	55	E84E (GAA→GAG)	yciW	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/29);  total (25/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 = 1.00e+00

GCGGGCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTG  >  NC_000913/1349709‑1350003                                                                                                                                                    |                                                                                                                                                    gCGGGCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCa                                                                                                                                                  <  1:937225/151‑1 (MQ=255)  cGGGCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCAt                                                                                                                                                 <  1:124499/151‑1 (MQ=255)      cATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGc                                                                                                                                             <  2:885019/151‑1 (MQ=255)          ttGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTg                                                                                                                                         >  2:471416/1‑151 (MQ=255)           tGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTg                                                                                                                                         >  1:233171/1‑150 (MQ=255)             cGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGa                                                                                                                                        >  2:534800/1‑149 (MQ=255)               tGATTTGTGACAACCGATTATTTCTTACCCTGGTGGAATCAGGCCCGGGGAGCAACGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAg                                                                                                                                      <  2:854286/149‑1 (MQ=255)                         cAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACgg                                                                                                                          >  2:598698/1‑151 (MQ=255)                          aaCCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACggg                                                                                                                         >  1:845727/1‑151 (MQ=255)                             cGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTc                                                                                                                       >  2:529300/1‑150 (MQ=255)                               atgatTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACg                                                                                                                    >  1:589905/1‑151 (MQ=255)                                  atTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGt                                                                                                                  >  1:334481/1‑150 (MQ=255)                                  atTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGt                                                                                                                  >  1:837579/1‑150 (MQ=255)                                   ttACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTa                                                                                                                 >  1:181450/1‑150 (MQ=255)                                    tACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTaa                                                                                                                <  2:381931/150‑1 (MQ=255)                                     aCTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTGCATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCGCGGTAAGc                                                                                                              <  1:257757/151‑1 (MQ=255)                                      cTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCt                                                                                                             >  2:1127112/1‑151 (MQ=255)                                           ccTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCa                                                                                                        <  2:1013100/151‑1 (MQ=255)                                                      cAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATc                                                                                             >  1:282948/1‑151 (MQ=255)                                                          cccGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACt                                                                                         >  1:45245/1‑151 (MQ=255)                                                               ggAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAg                                                                                    <  2:342845/151‑1 (MQ=255)                                                                  gcggcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACGTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCAACTGCAAGaa                                                                                  <  1:623314/150‑1 (MQ=255)                                                                  gcggcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAac                                                                                 <  1:122741/151‑1 (MQ=255)                                                                         aaGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCa                                                                          <  1:760702/151‑1 (MQ=255)                                                                          aGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACgcgc                                                                           >  1:396932/1‑149 (MQ=255)                                                                             cgcgCCGAATAGTAGTCACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAggg                                                                       <  2:1074316/150‑1 (MQ=255)                                                                               cgcCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTcc                                                                    <  1:102216/151‑1 (MQ=255)                                                                                  cGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGc                                                                 <  2:548171/151‑1 (MQ=255)                                                                                  cGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGc                                                                 >  2:744865/1‑151 (MQ=255)                                                                                   gAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGcc                                                                <  2:950400/151‑1 (MQ=255)                                                                                    aaTAGCAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGAGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGcc                                                                >  2:256312/1‑150 (MQ=255)                                                                                     aTAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCAt                                                              <  1:110690/151‑1 (MQ=255)                                                                                      tagtagTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATg                                                             >  2:848254/1‑151 (MQ=255)                                                                                        gtagtTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGaa                                                           <  1:956765/151‑1 (MQ=255)                                                                                         tagtTACATAACCGCTGGTCGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGGGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGACATGaa                                                           >  1:1093438/1‑150 (MQ=255)                                                                                                aTAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAggg                                                   <  1:519132/151‑1 (MQ=255)                                                                                                     cGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCa                                              <  1:689045/151‑1 (MQ=255)                                                                                                          gACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGtt                                         <  2:853959/151‑1 (MQ=255)                                                                                                          gACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGT‑CAGCCATGAAAGGAAGGGGGCCAGCGttt                                        >  2:992812/1‑151 (MQ=255)                                                                                                             gccgcAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGATTCt                                      >  1:446213/1‑151 (MQ=255)                                                                                                               cgcAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCttt                                    <  2:592245/151‑1 (MQ=255)                                                                                                                gcAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCtttt                                   >  2:1171006/1‑151 (MQ=255)                                                                                                                 cAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCttttt                                  >  2:111427/1‑151 (MQ=255)                                                                                                                  aaaCCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTc                                 <  2:141812/151‑1 (MQ=255)                                                                                                                  aaaCCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTc                                 <  2:336547/151‑1 (MQ=255)                                                                                                                        gAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAgg                             <  1:91592/149‑1 (MQ=255)                                                                                                                          gCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCa                         <  1:154972/151‑1 (MQ=255)                                                                                                                           cAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCAc                        <  2:809935/151‑1 (MQ=255)                                                                                                                            aaCCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACg                       <  2:835005/151‑1 (MQ=255)                                                                                                                              ccGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTc                     >  1:1082250/1‑151 (MQ=255)                                                                                                                                            gACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGAt       >  2:814122/1‑151 (MQ=255)                                                                                                                                             aCAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATc      <  1:22832/151‑1 (MQ=255)                                                                                                                                              cAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATc      <  1:522749/150‑1 (MQ=255)                                                                                                                                               aaGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGc    >  2:1006507/1‑151 (MQ=255)                                                                                                                                                aGCGTTCATAAACGGTGAAGGTCTGTGAACGGGTCACGGTAAGCCGGTCAGGGAACAAATCCGCTACAAGAACACGCGCAGGGTCCAGCCATGGAAGGAAGGGGGCCACCGTTTCTTTTTCAAGGATCACGACGAGTAGAAAGGGATCGCt   <  1:325519/151‑1 (MQ=255)                                                                                                                                                 gcgcTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTg  <  1:438075/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCGGGCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTG  >  NC_000913/1349709‑1350003

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