Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,542,775	T→G	intergenic (+142/‑340)	fimE → / → fimA	regulator for fimA/type 1 fimbriae major subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,775	0	T	G	100.0%	244.2 / NA	63	intergenic (+142/‑340)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/33);  total (30/33)

ACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAA  >  NC_000913/4542638‑4542911                                                                                                                                          |                                                                                                                                         aCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGagat                                                                                                                                        >  2:1066000/1‑140 (MQ=255) aCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGagat                                                                                                                                        >  1:300103/1‑140 (MQ=255)                    aaaaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTa                                                                                                                     >  2:957942/1‑140 (MQ=255)                             aataaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAAttt                                                                                                            <  1:18974/140‑1 (MQ=255)                                    aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCa                                                                                                     >  2:981617/1‑140 (MQ=255)                                     cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCAc                                                                                                    <  1:1073211/140‑1 (MQ=255)                                     cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCAc                                                                                                    <  2:300103/140‑1 (MQ=255)                                       aaaCATTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                                                  <  2:563774/138‑1 (MQ=255)                                       aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                                                  <  2:209615/140‑1 (MQ=255)                                        gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc                                                                                                   <  1:237375/138‑1 (MQ=255)                                        gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc                                                                                                   >  2:237375/1‑138 (MQ=255)                                        gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTc                                                                                                 >  2:894976/1‑140 (MQ=255)                                            aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCt                                                                                             <  2:502480/140‑1 (MQ=255)                                             ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                                                  >  1:127209/1‑134 (MQ=255)                                             ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt                                                                                                  <  2:127209/134‑1 (MQ=255)                                             ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCta                                                                                            >  2:948839/1‑140 (MQ=255)                                             ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCt                                                                                             <  2:503323/139‑1 (MQ=255)                                             ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCt                                                                                             >  1:503323/1‑139 (MQ=255)                                                  gCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGt                                                                                       <  2:778197/140‑1 (MQ=255)                                                  gCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGt                                                                                       >  2:984554/1‑140 (MQ=255)                                                           tCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGc                                                                                              <  1:867948/124‑1 (MQ=255)                                                           tCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGc                                                                                              >  2:867948/1‑124 (MQ=255)                                                           tCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAAc                                                                              >  1:923035/1‑140 (MQ=255)                                                           tCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAAc                                                                              >  1:115054/1‑140 (MQ=255)                                                            ccATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACg                                                                             >  2:172322/1‑140 (MQ=255)                                                             cATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGt                                                                            >  1:768364/1‑140 (MQ=255)                                                             cATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGt                                                                            >  1:602142/1‑140 (MQ=255)                                                               tatCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGttt                                                                          <  1:261166/140‑1 (MQ=255)                                                                atCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTc                                                                         >  2:846543/1‑140 (MQ=255)                                                                 tCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCt                                                                        <  2:937544/140‑1 (MQ=255)                                                                  cATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCtg                                                                       <  1:301152/140‑1 (MQ=255)                                                                     aaataaGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCt                                                                                 >  2:713365/1‑127 (MQ=255)                                                                     aaataaGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCt                                                                                 <  1:713365/127‑1 (MQ=255)                                                                     aaataaGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTgg                                                                    >  1:1079557/1‑140 (MQ=255)                                                                        taaGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTc                                                                 <  2:197914/140‑1 (MQ=255)                                                                         aaGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCg                                                                >  2:86114/1‑140 (MQ=255)                                                                          aGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGa                                                               >  1:950822/1‑140 (MQ=255)                                                                           gTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGAc                                                              <  2:602142/140‑1 (MQ=255)                                                                            ttACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACg                                                             >  1:856514/1‑140 (MQ=255)                                                                            ttACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACg                                                             >  2:426688/1‑140 (MQ=255)                                                                               cGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCAt                                                          <  2:115054/140‑1 (MQ=255)                                                                               cGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCAt                                                          <  1:749354/140‑1 (MQ=255)                                                                               cGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCAt                                                          >  2:362847/1‑140 (MQ=255)                                                                               cGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCAt                                                          >  2:555459/1‑140 (MQ=255)                                                                                            aaGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCAttctt                                             <  1:172322/140‑1 (MQ=255)                                                                                            aaGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCAttctt                                             <  2:950822/140‑1 (MQ=255)                                                                                            aaGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCAttctt                                             <  1:948839/140‑1 (MQ=255)                                                                                             aGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCAttcttc                                            <  2:923035/140‑1 (MQ=255)                                                                                              gCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATtcttct                                           <  1:1066000/140‑1 (MQ=255)                                                                                              gCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATtcttct                                           >  2:1028096/1‑140 (MQ=255)                                                                                              gCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATtcttct                                           >  1:590428/1‑140 (MQ=255)                                                                                                  aaaaaTATTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTcc                                       <  2:337665/140‑1 (MQ=255)                                                                                                        atTAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAac                                 >  2:759006/1‑140 (MQ=255)                                                                                                         ttAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAacc                                >  2:872379/1‑140 (MQ=255)                                                                                                          tAAAAAACGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAacca                               >  1:407745/1‑140 (MQ=255)                                                                                                            aaaaaCGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAaccacc                             <  1:555459/140‑1 (MQ=255)                                                                                                            aaaaaCGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAaccacc                             <  1:894976/140‑1 (MQ=255)                                                                                                            aaaaaCGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAaccacc                             >  1:1147491/1‑140 (MQ=255)                                                                                                                aCGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCAt                         <  1:984554/140‑1 (MQ=255)                                                                                                                 cGACAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATg                        <  2:407745/140‑1 (MQ=255)                                                                                                                   aCAAAAAGCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCa                      <  1:1028096/140‑1 (MQ=255)                                                                                                                          gCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa               <  2:590428/140‑1 (MQ=255)                                                                                                                          gCATCTAACTGTTTGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa               <  1:981617/140‑1 (MQ=255)                                                                                                                                      ttGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaata   <  1:821883/140‑1 (MQ=255)                                                                                                                                       tGAGATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTataaataa  >  1:994109/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                         ACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAA  >  NC_000913/4542638‑4542911

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