Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	54,480	T→A	I75F (ATC→TTC)	surA ←	peptidyl‑prolyl cis‑trans isomerase SurA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	54,480	0	T	A	100.0%	220.4 / NA	62	I75F (ATC→TTC)	surA	peptidyl‑prolyl cis‑trans isomerase SurA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (29/33);  total (29/33)

AAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGAC  >  NC_000913/54345‑54612                                                                                                                                        |                                                                                                                                     aaGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACaaa                                                                                                                                  <  2:411477/140‑1 (MQ=255)   gCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACg                                                                                                                                <  2:33038/140‑1 (MQ=255)   gCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACg                                                                                                                                >  1:837323/1‑140 (MQ=255)                       aTCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCa                                                                                                            <  2:705175/140‑1 (MQ=255)                         ccAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAgc                                                                                                          >  2:315725/1‑140 (MQ=255)                              gTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAgcgtcgc                                                                                                     >  2:541383/1‑140 (MQ=255)                               tCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAgcgtcgcg                                                                                                    <  2:1101406/140‑1 (MQ=255)                               tCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAgcgtcgcg                                                                                                    <  2:429150/140‑1 (MQ=255)                                 aTGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgtc                                                                                                  >  1:498596/1‑140 (MQ=255)                                 aTGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgtc                                                                                                  >  1:778073/1‑140 (MQ=255)                                  tgttgtTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtca                                                                                                 >  2:1149289/1‑140 (MQ=255)                                   gttgttCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAATTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcat                                                                                                <  2:80486/140‑1 (MQ=255)                                      gttCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAg                                                                                             <  2:674459/140‑1 (MQ=255)                                         cTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGaa                                                                                          >  2:641910/1‑140 (MQ=255)                                         cTGTTTCACGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGaa                                                                                          >  2:516027/1‑140 (MQ=255)                                            tttCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt                                                                                       >  1:641564/1‑140 (MQ=255)                                              tcacgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgc                                                                                     <  1:516027/137‑1 (MQ=255)                                                 cgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctgc                                                                                  <  1:541383/140‑1 (MQ=255)                                                 cgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctgc                                                                                  <  1:282569/140‑1 (MQ=255)                                                 cgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctgc                                                                                  <  1:641910/140‑1 (MQ=255)                                                 cgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctgc                                                                                  <  1:849749/140‑1 (MQ=255)                                                  gATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcc                                                                                 <  1:1149289/140‑1 (MQ=255)                                                              cGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGagcag                                                                     >  2:514152/1‑140 (MQ=255)                                                              cGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGagcag                                                                     >  1:368241/1‑140 (MQ=255)                                                                ccTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCg                                                                   >  2:716845/1‑140 (MQ=255)                                                                ccTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCg                                                                   >  2:715812/1‑140 (MQ=255)                                                                    aTCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCa                                                               >  1:2782/1‑140 (MQ=255)                                                                         gctgctCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGtttt                                                          <  2:778073/140‑1 (MQ=255)                                                                          ctgctCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTa                                                         >  2:1044253/1‑140 (MQ=255)                                                                            gctCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTAcc                                                       <  2:961918/140‑1 (MQ=255)                                                                              tCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGa                                                     <  1:156102/140‑1 (MQ=255)                                                                              tCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGa                                                     >  2:946909/1‑140 (MQ=255)                                                                              tCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGa                                                     <  2:4369/140‑1 (MQ=255)                                                                               cATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGAc                                                    <  2:2782/140‑1 (MQ=255)                                                                               cATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGAc                                                    >  1:570361/1‑140 (MQ=255)                                                                               cATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGAc                                                    >  1:615047/1‑140 (MQ=255)                                                                                  cGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGc                                                 >  2:1112703/1‑140 (MQ=255)                                                                                          ttCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATcc                                         <  2:830762/140‑1 (MQ=255)                                                                                             aCTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATc                                      <  2:615047/140‑1 (MQ=255)                                                                                                    ttttCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCg                               <  2:702014/140‑1 (MQ=255)                                                                                                    ttttCTGCCCCATCTGCAGGATGATTAGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCg                               <  1:812454/140‑1 (MQ=255)                                                                                                     tttCTGCCCCATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGc                              <  1:1112703/140‑1 (MQ=255)                                                                                                           ccccATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCg                               >  1:557850/1‑133 (MQ=255)                                                                                                           ccccATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCg                               <  2:557850/133‑1 (MQ=255)                                                                                                              cATCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCa                     <  2:952240/140‑1 (MQ=255)                                                                                                               aTCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCac                    <  2:498596/140‑1 (MQ=255)                                                                                                                tCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacg                   >  2:583550/1‑140 (MQ=255)                                                                                                                tCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacg                   >  1:1044714/1‑140 (MQ=255)                                                                                                                tCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacg                   >  1:724629/1‑140 (MQ=255)                                                                                                                tCTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacg                   >  2:1038775/1‑140 (MQ=255)                                                                                                                 cTGCAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacga                  <  1:475085/140‑1 (MQ=255)                                                                                                                    cAGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGc               <  2:1044714/140‑1 (MQ=255)                                                                                                                     aGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacgacg                >  1:720296/1‑138 (MQ=255)                                                                                                                     aGGATGATTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCacgacg                <  2:720296/138‑1 (MQ=255)                                                                                                                       gatgatTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGt            >  2:635954/1‑140 (MQ=255)                                                                                                                          gatTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTAt         >  2:686610/1‑140 (MQ=255)                                                                                                                          gatTTGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTAt         >  1:87422/1‑140 (MQ=255)                                                                                                                            tttGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTg       <  1:912986/140‑1 (MQ=255)                                                                                                                              tGATCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTgac     <  2:837323/140‑1 (MQ=255)                                                                                                                                aTCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTgacga   <  1:635954/140‑1 (MQ=255)                                                                                                                                 tCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTgacgac  >  2:483210/1‑140 (MQ=255)                                                                                                                                 tCCATGAACAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTgacgac  <  1:1044253/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     AAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGAC  >  NC_000913/54345‑54612

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