Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,628	A→C	N109H (AAC→CAC)	fimC →	periplasmic chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,628	0	A	C	100.0%	230.9 / NA	60	N109H (AAC→CAC)	fimC	periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/33);  total (27/33)

AATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCG  >  NC_000913/4544484‑4544774                                                                                                                                                 |                                                                                                                                                   aatGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCaa                                                                                                                                               <  1:677482/150‑1 (MQ=255) aatGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAAt                                                                                                                                              <  2:919836/151‑1 (MQ=255) aatGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAAt                                                                                                                                              <  1:331938/151‑1 (MQ=255)       gaAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCcac                                                                                                                                        >  1:467173/1‑151 (MQ=255)        aaaaTAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCcac                                                                                                                                        <  1:235382/150‑1 (MQ=255)         aaaTAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCcac                                                                                                                                        >  1:908112/1‑149 (MQ=255)           aTAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAgg                                                                                                                                     >  1:615612/1‑150 (MQ=255)             aGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGAcc                                                                                                                                  >  2:56946/1‑151 (MQ=255)              gTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCg                                                                                                                                 <  2:242910/151‑1 (MQ=255)                aCCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCggg                                                                                                                               >  1:828618/1‑151 (MQ=255)                   tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGaa                                                                                                                             >  2:312715/1‑150 (MQ=255)                    aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGaa                                                                                                                             <  2:550198/149‑1 (MQ=255)                    aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAg                                                                                                                           <  2:542182/151‑1 (MQ=255)                      ttAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGt                                                                                                                          <  1:17262/150‑1 (MQ=255)                       tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGtt                                                                                                                         <  1:237615/150‑1 (MQ=255)                              aTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                                                                                                                 >  1:773587/1‑151 (MQ=255)                                cATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                                                                                                               <  2:105725/151‑1 (MQ=255)                                   gggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGaa                                                                                                            <  2:170880/151‑1 (MQ=255)                                                 gATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc                                                                                              <  2:472691/151‑1 (MQ=255)                                                     gtgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc                                                                                          <  2:804609/151‑1 (MQ=255)                                                          aaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATgg                                                                                     >  1:755367/1‑151 (MQ=255)                                                                    gTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCaaa                                                                           <  2:564827/151‑1 (MQ=255)                                                                      tttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCaaaa                                                                          <  2:645959/150‑1 (MQ=255)                                                                                 gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGaa                                                               >  1:744947/1‑150 (MQ=255)                                                                                 gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAAt                                                              <  2:134427/151‑1 (MQ=255)                                                                                    tcctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACg                                                           >  1:31147/1‑151 (MQ=255)                                                                                     cctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATAc                                                            <  2:337038/149‑1 (MQ=255)                                                                                     cctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACg                                                           <  2:256025/150‑1 (MQ=255)                                                                                     cctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGc                                                          >  2:170921/1‑151 (MQ=255)                                                                                        ctGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTa                                                        >  1:430882/1‑150 (MQ=255)                                                                                        ctGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTAc                                                       >  1:62875/1‑151 (MQ=255)                                                                                          gTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTa                                                        <  2:283622/148‑1 (MQ=255)                                                                                          gTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTAc                                                       >  1:71315/1‑149 (MQ=255)                                                                                          gTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACa                                                      <  2:416612/150‑1 (MQ=255)                                                                                          gTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAg                                                     <  1:31706/151‑1 (MQ=255)                                                                                           tttGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGc                                                    >  2:48012/1‑151 (MQ=255)                                                                                              gCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTc                                                  <  1:168043/150‑1 (MQ=255)                                                                                                gATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGc                                                >  1:603578/1‑150 (MQ=255)                                                                                                  tGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAAt                                             >  2:69105/1‑151 (MQ=255)                                                                                                    aaGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTa                                           <  2:828600/151‑1 (MQ=255)                                                                                                       ggAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATc                                         >  2:705747/1‑150 (MQ=255)                                                                                                        gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAg                                       <  2:619125/151‑1 (MQ=255)                                                                                                         aaaaaaaGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGc                                      <  1:868217/151‑1 (MQ=255)                                                                                                          aaaaaaGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGcc                                     >  1:78883/1‑151 (MQ=255)                                                                                                               agagAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCAtt                                <  2:522742/151‑1 (MQ=255)                                                                                                                   aaTACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAAc                            <  1:96969/151‑1 (MQ=255)                                                                                                                       ccTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGt                         >  2:649757/1‑150 (MQ=255)                                                                                                                       ccTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTa                        <  1:925800/151‑1 (MQ=255)                                                                                                                           aCGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTa                     <  2:428788/150‑1 (MQ=255)                                                                                                                            cGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATc                   <  2:607598/151‑1 (MQ=255)                                                                                                                              tATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGc                 >  2:805132/1‑151 (MQ=255)                                                                                                                               aTTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGcc                <  1:514849/151‑1 (MQ=255)                                                                                                                                 tCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCg              >  1:392709/1‑151 (MQ=255)                                                                                                                                  cTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCgg             >  2:854071/1‑151 (MQ=255)                                                                                                                                    tGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCg              <  1:555710/148‑1 (MQ=255)                                                                                                                                      aTGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTaa         <  2:31147/151‑1 (MQ=255)                                                                                                                                      aTGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTa          >  1:557330/1‑150 (MQ=255)                                                                                                                                         caacaaATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAAtt      >  2:685289/1‑151 (MQ=255)                                                                                                                                         caacaaATCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAAtt      >  1:100536/1‑151 (MQ=255)                                                                                                                                             aaaTCACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCg  >  1:188392/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCG  >  NC_000913/4544484‑4544774

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