Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,100,871	A→G	S272P (TCG→CCG)	csgG ←	curli secretion channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,100,871	0	A	G	100.0%	224.1 / NA	62	S272P (TCG→CCG)	csgG	curli secretion channel
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/35);  total (27/35)

GCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAA  >  NC_000913/1100734‑1101005                                                                                                                                          |                                                                                                                                       gcgcGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCa                                                                                                                                      <  1:1153201/140‑1 (MQ=255)   gcgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCata                                                                                                                                    <  1:241673/140‑1 (MQ=255)    cgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCatat                                                                                                                                   >  1:862510/1‑140 (MQ=255)    cgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCatat                                                                                                                                   >  1:1127544/1‑140 (MQ=255)    cgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCatat                                                                                                                                   <  1:956057/140‑1 (MQ=255)           cAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTa                                                                                                                            <  1:1333417/140‑1 (MQ=255)           cAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTa                                                                                                                            >  1:932802/1‑140 (MQ=255)              aGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTAttt                                                                                                                         <  2:1022768/140‑1 (MQ=255)                 cGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCAc                                                                                                                      <  1:802398/140‑1 (MQ=255)                 cGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCAc                                                                                                                      <  1:1139727/140‑1 (MQ=255)                   gcgcACCGTGAGCAGGAAATCGAAAATAAAGCGCGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCa                                                                                                                    >  2:560996/1‑140 (MQ=255)                   gcgcACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCa                                                                                                                    >  2:1028587/1‑140 (MQ=255)                   gcgcACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCa                                                                                                                    >  2:516865/1‑140 (MQ=255)                                   aaaTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCg                                                                                                    >  1:1224184/1‑140 (MQ=255)                                   aaaTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCg                                                                                                    <  1:1289930/140‑1 (MQ=255)                                   aaaTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCg                                                                                                    <  1:516865/140‑1 (MQ=255)                                     aTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGtt                                                                                                  >  2:159401/1‑140 (MQ=255)                                     aTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGtt                                                                                                  >  1:790828/1‑140 (MQ=255)                                     aTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGtt                                                                                                  <  2:786175/140‑1 (MQ=255)                                              aaaGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAt                                                                                         <  2:369131/140‑1 (MQ=255)                                                      aGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCaaa                                                                                 <  1:872095/140‑1 (MQ=255)                                                           cACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCca                                                                            <  1:1028587/140‑1 (MQ=255)                                                             cTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc                                                                             <  1:1034694/137‑1 (MQ=255)                                                             cTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc                                                                             >  2:1034694/1‑137 (MQ=255)                                                             cTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca                                                                          >  1:181083/1‑140 (MQ=255)                                                             cTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca                                                                          >  2:966003/1‑140 (MQ=255)                                                                 cGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAcc                                                                      <  2:469256/140‑1 (MQ=255)                                                                  gCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCa                                                                     <  1:844663/140‑1 (MQ=255)                                                                  gCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCa                                                                     <  2:181083/140‑1 (MQ=255)                                                                  gCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCa                                                                     <  2:1032093/140‑1 (MQ=255)                                                                   cTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCAc                                                                    <  1:292444/140‑1 (MQ=255)                                                                     ttCGCTTAAACAGTAAAATGTCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACgg                                                                  >  2:889725/1‑140 (MQ=255)                                                                       cGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTc                                                                >  1:281869/1‑140 (MQ=255)                                                                        gCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCg                                                               >  1:798888/1‑140 (MQ=255)                                                                        gCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCg                                                               >  1:143044/1‑140 (MQ=255)                                                                        gCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCg                                                               >  1:901296/1‑140 (MQ=255)                                                                           tAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATa                                                            <  2:798888/140‑1 (MQ=255)                                                                            aaaCAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATAc                                                           <  1:676539/140‑1 (MQ=255)                                                                            aaaCAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATAc                                                           >  1:8488/1‑140 (MQ=255)                                                                              aCAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACca                                                         >  1:178951/1‑140 (MQ=255)                                                                                aGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACcatc                                                       >  1:257800/1‑140 (MQ=255)                                                                                 gTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACcatca                                                      <  2:1127544/140‑1 (MQ=255)                                                                                 gTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCAACAGACCACGGTCGATACcatca                                                      <  1:1190258/140‑1 (MQ=255)                                                                                       tGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATc                                                <  2:977215/140‑1 (MQ=255)                                                                                        gCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCa                                               <  2:862510/140‑1 (MQ=255)                                                                                            gatgatAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGaa                                           >  1:431722/1‑140 (MQ=255)                                                                                               gatAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAAt                                        <  1:1003226/140‑1 (MQ=255)                                                                                                  aaTTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAgg                                             >  2:91158/1‑132 (MQ=255)                                                                                                  aaTTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAgg                                             <  1:91158/132‑1 (MQ=255)                                                                                                   aTTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGAcc                                    <  2:901296/140‑1 (MQ=255)                                                                                                   aTTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGAcc                                    >  1:467753/1‑140 (MQ=255)                                                                                                       cGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTg                                >  1:383797/1‑140 (MQ=255)                                                                                                           ttttttATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTTATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTc                            >  2:174852/1‑140 (MQ=255)                                                                                                                tATCTGTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAg                       <  2:257800/140‑1 (MQ=255)                                                                                                                     gTCAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGAc                  >  2:795189/1‑140 (MQ=255)                                                                                                                       cAGGATTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACAt                <  1:159401/140‑1 (MQ=255)                                                                                                                           aTTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAgg            <  2:1224184/140‑1 (MQ=255)                                                                                                                           aTTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAgg            <  1:966003/140‑1 (MQ=255)                                                                                                                           aTTCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAgg            <  1:560996/140‑1 (MQ=255)                                                                                                                             tCCGGTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGca          <  2:431722/140‑1 (MQ=255)                                                                                                                                ggTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAgg            <  2:16664/135‑1 (MQ=255)                                                                                                                                ggTGGAACCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAgg            >  1:16664/1‑135 (MQ=255)                                                                                                                                     aaCCGGCATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATaa  >  2:63362/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       GCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCACGGTCGATACCATCATTAATCAGGAAAATGACCCCTGTTTCGATAGCCGACATCAGGCACAGCATAA  >  NC_000913/1100734‑1101005

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