Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,165,046	A→G	T470A (ACC→GCC)	btuB →	cobalamin/cobinamide outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,165,046	0	A	G	100.0%	226.6 / NA	65	T470A (ACC→GCC)	btuB	cobalamin/cobinamide outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/31);  total (34/31)

CTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACACCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACGCC  >  NC_000913/4164910‑4165179                                                                                                                                         |                                                                                                                                      cTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGccc                                                                                                                                    >  1:445091/1‑140 (MQ=255) cTACGGAAATCAGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGccc                                                                                                                                    <  1:1087719/140‑1 (MQ=255)   aCGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGccc                                                                                                                                    >  1:805374/1‑138 (MQ=255)   aCGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGccc                                                                                                                                    <  2:805374/138‑1 (MQ=255)   aCGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTg                                                                                                                                  >  1:1112743/1‑140 (MQ=255)       aaaTCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAat                                                                                                                              >  2:398060/1‑140 (MQ=255)           ccGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTa                                                                                                                          <  2:445091/140‑1 (MQ=255)            cGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTAc                                                                                                                         <  2:136867/140‑1 (MQ=255)              aaTCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACaa                                                                                                                       >  2:483575/1‑140 (MQ=255)                    gACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAgg                                                                                                                 >  1:1069893/1‑140 (MQ=255)                           agaAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAgcg                                                                                                          <  2:510190/140‑1 (MQ=255)                           agaAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAgcg                                                                                                          <  2:810676/140‑1 (MQ=255)                           agaAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAgcg                                                                                                          <  2:811751/140‑1 (MQ=255)                           agaAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAgcg                                                                                                          <  1:744953/140‑1 (MQ=255)                                     aCAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTTACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAgg                                                                                                >  2:873279/1‑140 (MQ=255)                                     aCAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAgg                                                                                                <  2:1069893/140‑1 (MQ=255)                                     aCAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAgg                                                                                                >  1:909866/1‑140 (MQ=255)                                                 cgcgTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGa                                                                                    >  2:640904/1‑140 (MQ=255)                                                 cgcgTTTGAAGGCTTAACCGCTGGGGTGAACTGACGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGa                                                                                    >  2:342023/1‑140 (MQ=255)                                                   cgTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCTATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGAcc                                                                                  >  1:154362/1‑140 (MQ=255)                                                      ttGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATCATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGcc                                                                               <  2:1025212/140‑1 (MQ=255)                                                           ggCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCa                                                                              >  2:545837/1‑136 (MQ=255)                                                           ggCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCa                                                                              <  1:545837/136‑1 (MQ=255)                                                               tAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGAt                                                                      <  2:909866/140‑1 (MQ=255)                                                                aaCCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATa                                                                     >  1:904128/1‑140 (MQ=255)                                                                 aCCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATAc                                                                    <  1:483575/140‑1 (MQ=255)                                                                     cTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGa                                                                >  1:673251/1‑140 (MQ=255)                                                                     cTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGa                                                                >  2:426140/1‑140 (MQ=255)                                                                      tGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGAc                                                               >  1:1012674/1‑140 (MQ=255)                                                                       ggggTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGAcc                                                              >  1:89547/1‑140 (MQ=255)                                                                          gTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACt                                                           >  1:406626/1‑140 (MQ=255)                                                                              aCTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACg                                                       <  2:1012674/140‑1 (MQ=255)                                                                               cTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGc                                                      <  2:904128/140‑1 (MQ=255)                                                                                tGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCa                                                     <  1:648488/140‑1 (MQ=255)                                                                                 ggCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCAt                                                    <  1:426140/140‑1 (MQ=255)                                                                                  gCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATa                                                   >  1:405390/1‑140 (MQ=255)                                                                                  gCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTAAGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGc                                                      >  1:357384/1‑137 (MQ=255)                                                                                  gCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTAAGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGc                                                      <  2:357384/137‑1 (MQ=255)                                                                                   cGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCa                                                             <  1:151053/129‑1 (MQ=255)                                                                                   cGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCa                                                             >  2:151053/1‑129 (MQ=255)                                                                                   cGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATAc                                                  <  1:59318/140‑1 (MQ=255)                                                                                         tCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACt                                                           <  2:189216/125‑1 (MQ=255)                                                                                         tCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACt                                                           >  1:189216/1‑125 (MQ=255)                                                                                          ccGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGt                                           <  1:640904/140‑1 (MQ=255)                                                                                             gATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTAt                                        <  2:406626/140‑1 (MQ=255)                                                                                               tatCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGa                                      <  2:488022/140‑1 (MQ=255)                                                                                                    tAACGATGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATg                                 <  1:32844/140‑1 (MQ=255)                                                                                                          tGTCAGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATg                           >  2:308871/1‑140 (MQ=255)                                                                                                              aGTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATgcgcg                       >  2:762817/1‑140 (MQ=255)                                                                                                               gTGACTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATgcgcgc                      <  1:762817/140‑1 (MQ=255)                                                                                                                   cTTGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATg                  <  2:673251/140‑1 (MQ=255)                                                                                                                    ttGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGc                 >  1:334023/1‑140 (MQ=255)                                                                                                                    ttGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGc                 >  2:763648/1‑140 (MQ=255)                                                                                                                     tGATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCg                >  2:846303/1‑140 (MQ=255)                                                                                                                      gATCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGa               >  1:1124235/1‑140 (MQ=255)                                                                                                                        tCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACTGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGAtt             >  2:546489/1‑140 (MQ=255)                                                                                                                        tCGATTATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGAtt             >  1:830932/1‑140 (MQ=255)                                                                                                                            ttATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATgcg                         >  2:794462/1‑124 (MQ=255)                                                                                                                            ttATGATGATCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATgcg                         <  1:794462/124‑1 (MQ=255)                                                                                                                                 atgatCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCGTACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACg    >  1:672976/1‑140 (MQ=255)                                                                                                                                 atgatCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACg    >  2:531828/1‑140 (MQ=255)                                                                                                                                  tgatCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCGTACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACGc   <  2:672976/140‑1 (MQ=255)                                                                                                                                  tgatCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACGc   <  1:398060/140‑1 (MQ=255)                                                                                                                                  tgatCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACGc   >  2:117291/1‑140 (MQ=255)                                                                                                                                   gatCACGCCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACGcc  <  1:308871/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      CTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCGTTTGAAGGCTTAACCGCTGGGGTGAACTGGCGTATTTCCGGATATCGTAACGATGTCAGTGACTTGATCGATTATGATGATCACACCCTGAAATATTACAACGAAGGGAAAGCGCGGATTAAGGGCGTCGAGGCGACCGCCAATTTTGATACCGGACCACTGACGCATACTGTGAGTTATGATTATGTCGATGCGCGCAATGCGATTACCGACACGCC  >  NC_000913/4164910‑4165179

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