Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,303,774	Δ1 bp	coding (327/1008 nt)	yhbW →	putative luciferase‑like monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,303,771	0	G	.	100.0%	252.3 / NA	59	coding (324/1008 nt)	yhbW	putative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (31/28);  total (31/28)

ATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTTGGGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCC  >  NC_000913/3303633‑3303899                                                                                                                                           |                                                                                                                                 aTTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTTgg                                                                                                                                 >  1:608748/1‑140 (MQ=255)  ttGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTTggg                                                                                                                                <  2:606743/140‑1 (MQ=255)     aTCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTTgggctg                                                                                                                             >  2:63337/1‑137 (MQ=255)      tCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTTgggctgg                                                                                                                            >  1:21641/1‑136 (MQ=255)         gCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTc                                                                                                                        >  2:98025/1‑140 (MQ=255)         gCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTc                                                                                                                        <  2:724689/140‑1 (MQ=255)         gCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTc                                                                                                                        <  1:933119/140‑1 (MQ=255)               tGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTc                                                                                                                  <  1:692245/140‑1 (MQ=255)                       aaTACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGt                                                                                                          <  2:52718/140‑1 (MQ=255)                            caccacGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCa                                                                                                     <  2:806165/140‑1 (MQ=255)                             accacGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCaa                                                                                                    <  1:888210/140‑1 (MQ=255)                              ccacGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAAc                                                                                                   >  2:339738/1‑140 (MQ=255)                                 cGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGa                                                                                                >  2:653590/1‑140 (MQ=255)                                 cGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGa                                                                                                >  2:824714/1‑140 (MQ=255)                                       aTCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAatga                                                                                          <  2:608748/140‑1 (MQ=255)                                           ggggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAatga                                                                                          >  2:898937/1‑136 (MQ=255)                                           ggggTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAatga                                                                                          <  1:898937/136‑1 (MQ=255)                                              gTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCt                                                                                   <  1:105129/140‑1 (MQ=255)                                              gTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCt                                                                                   <  2:470190/140‑1 (MQ=255)                                              gTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCt                                                                                   <  1:865503/140‑1 (MQ=255)                                              gTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCt                                                                                   <  2:773746/140‑1 (MQ=255)                                                cTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTAc                                                                                 <  1:248715/140‑1 (MQ=255)                                                cTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTAc                                                                                 <  2:250339/140‑1 (MQ=255)                                                 tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCt                                                                                   >  1:280205/1‑137 (MQ=255)                                                 tGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCt                                                                                   <  2:280205/137‑1 (MQ=255)                                                      gcgTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACgtcgtc                                                                           <  1:81909/140‑1 (MQ=255)                                                             gTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTAc                                                                                 >  2:25076/1‑127 (MQ=255)                                                             gTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTAc                                                                                 <  1:25076/127‑1 (MQ=255)                                                                      ccACTCACCGTTGGTCATTGCAGAACAGTTCGACACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGa                                                                     <  2:499287/130‑1 (MQ=255)                                                                      ccACTCACCGTTGGTCATTGCAGAACAGTTCGACACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGa                                                                     >  1:499287/1‑130 (MQ=255)                                                                          tCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGAt                                                       >  1:498358/1‑140 (MQ=255)                                                                           cACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATa                                                      >  1:382923/1‑140 (MQ=255)                                                                            aCCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATaa                                                     >  1:589360/1‑140 (MQ=255)                                                                                ttGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTc                                                 >  1:980828/1‑140 (MQ=255)                                                                                 tGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTcc                                                >  1:429235/1‑140 (MQ=255)                                                                                        tGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGa                                         >  1:701105/1‑140 (MQ=255)                                                                                        tGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGa                                         >  1:793268/1‑140 (MQ=255)                                                                                                  ttCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAg                               >  2:102373/1‑140 (MQ=255)                                                                                                  ttCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAg                               >  1:806418/1‑140 (MQ=255)                                                                                                  ttCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAg                               >  1:350/1‑140 (MQ=255)                                                                                                  ttCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAg                               >  2:57112/1‑140 (MQ=255)                                                                                                     ggCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCtg                            <  2:429235/140‑1 (MQ=255)                                                                                                     ggCACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCtg                            <  2:793268/140‑1 (MQ=255)                                                                                                       cACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCgg                                 <  1:359564/133‑1 (MQ=255)                                                                                                       cACGCTTAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCgg                                 >  2:359564/1‑133 (MQ=255)                                                                                                            ttAATACACTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACt                     >  1:409301/1‑140 (MQ=255)                                                                                                                   acTCTATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGa              <  1:824714/140‑1 (MQ=255)                                                                                                                     tctATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACg            <  1:63337/140‑1 (MQ=255)                                                                                                                      ctATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGc           >  2:853377/1‑140 (MQ=255)                                                                                                                      ctATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGc           >  2:508399/1‑140 (MQ=255)                                                                                                                      ctATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGc           >  1:227610/1‑140 (MQ=255)                                                                                                                       tATCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGcc          <  2:21641/140‑1 (MQ=255)                                                                                                                        aTCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACg            >  2:68600/1‑137 (MQ=255)                                                                                                                        aTCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACg            <  1:68600/137‑1 (MQ=255)                                                                                                                        aTCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGccc         >  1:798427/1‑140 (MQ=255)                                                                                                                        aTCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGccc         >  1:751561/1‑140 (MQ=255)                                                                                                                         tCCGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCcg        <  1:818870/140‑1 (MQ=255)                                                                                                                          ccGGGGCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCcgc       >  1:575353/1‑140 (MQ=255)                                                                                                                             gggCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTCTGACGCCCGCGAt    >  1:366874/1‑140 (MQ=255)                                                                                                                              ggCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATc   >  1:466790/1‑140 (MQ=255)                                                                                                                               gCGAATCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATcc  <  2:227610/140‑1 (MQ=255)                                                                                                                                   aTCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGa              >  2:166707/1‑124 (MQ=255)                                                                                                                                   aTCGATTT‑GGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGa              <  1:166707/124‑1 (MQ=255)                                                                                                                                           |                                                                                                                                 ATTGATCGGCTATCTGGCGGCGAATACCACCACGCTGCATCTGGGGTCTGGCGGCGTGATGTTGCCTAACCACTCACCGTTGGTCATTGCAGAACAGTTCGGCACGCTTAATACACTCTATCCGGGGCGAATCGATTTGGGGCTGGGTCGTGCTCCGGGTAGTGACCAACGGACAATGATGGCGCTACGTCGTCATATGAGCGGCGATATTGATAATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCC  >  NC_000913/3303633‑3303899

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