Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP028101	880,783	A→G	K1885K (AAA→AAG)	C7Y58_04335 →	autotransporter domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP028101	880,783	0	A	G	100.0%	120.2 / NA	39	K1885K (AAA→AAG)	C7Y58_04335	autotransporter domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/15);  total (24/15)

ATCAACAGGAATAATCTTAGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAA‑CATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGGA  >  CP028101/880537‑880885                                                                                                                                                                                                                                                       |                                                                                                        aTCAACAGGAATAATCTTAGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGata                                                                                                     >  1:334726/1‑251 (MQ=255)        ggAATAATCTTAGCTAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGCACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACATAGGTATATGTCATTATACATTTAAGTTCAAGGATATAGGaaa                                                                                              >  2:775252/1‑251 (MQ=17)         gAATAATCTTAGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGa                                                                                                <  1:548481/248‑1 (MQ=255)              atCTTAGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATACTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACATGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCaaa                                                                                        >  1:804444/1‑251 (MQ=255)                cTTAGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCaaaaa                                                                                      <  2:490726/251‑2 (MQ=255)                                    atCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAg                                                                 >  2:461782/1‑251 (MQ=255)                                    atCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAg                                                                 >  1:370958/1‑251 (MQ=255)                                     tCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGGAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTCGATGGTACACAGGTATAGTTCATAATACATTTAACTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAgg                                                                >  1:739370/1‑251 (MQ=14)                                            gaTGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGatt                                                         >  2:504018/1‑251 (MQ=255)                                               gAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGattatt                                                      >  2:231010/1‑251 (MQ=255)                                                 aTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATCGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTaa                                                    >  1:798096/1‑251 (MQ=255)                                                 aTGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTaa                                                    <  1:286194/251‑1 (MQ=255)                                                   ggAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTaaaa                                                  >  1:883191/1‑251 (MQ=255)                                                   ggAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGATGCAAGCAAAGGTAGGATTATTaaaa                                                  >  2:220478/1‑251 (MQ=255)                                                    gagaACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTaaaa                                                  <  1:461782/250‑1 (MQ=255)                                                      gaACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGTTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGACTATTAAACTCa                                               >  2:607398/1‑251 (MQ=14)                                                       aaCAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATc                                                <  2:101559/249‑1 (MQ=255)                                                         cAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGAttatta                                                     <  1:203703/242‑1 (MQ=255)                                                         cAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGAttatta                                                     >  2:203703/1‑242 (MQ=255)                                                           gCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGCATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGGATAGTTAATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAAACA‑ATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGca                                          >  2:532968/1‑250 (MQ=11)                                                           gCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTACGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGc                                           >  1:828486/1‑250 (MQ=255)                                                            cTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCa                                         >  2:796339/1‑251 (MQ=255)                                                            cTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCa                                         >  1:543362/1‑251 (MQ=255)                                                                        cAAATAACATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAAc                             <  1:835803/251‑1 (MQ=255)                                                                           ataaaataaaAGCATTTGGAACTAGAGGAGAATACAAAACAGATACTGCAGGAGTTATAGACTACTAGAATCATGCTTATGGAGTAGCTTATGTACATGAAATTGAAGATATTAAACTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAAt                          <  1:777737/251‑1 (MQ=11)                                                                           ataaaataaaAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAAt                          <  2:899525/251‑1 (MQ=255)                                                                           ataaaataaaAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAAt                          <  2:370958/251‑1 (MQ=255)                                                                           ataaaataaaAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAAt                          <  1:818225/251‑1 (MQ=255)                                                                                ataaaaGCATTTGGAACAAGAGGAGAATACAAAACAGATACAGCAGGAGATATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGTAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGAAAGCAAAGGTAGGTTTATTAAAATCAGCGCAATTTGATGATAACATTAGCtt                     >  2:744043/1‑251 (MQ=11)                                                                                    aaGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTGATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGTGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAAt                 <  1:775252/251‑1 (MQ=255)                                                                                      gCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTg               >  2:569215/1‑251 (MQ=255)                                                                                         tttGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAACAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGAc            <  1:638722/251‑1 (MQ=255)                                                                                           tGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACaa          >  1:526017/1‑251 (MQ=16)                                                                                             gAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAata        >  1:83957/1‑251 (MQ=16)                                                                                              aaCTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAatat       >  2:406275/1‑251 (MQ=12)                                                                                              aaCTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAatat       >  1:65392/1‑251 (MQ=12)                                                                                                cTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCa     >  2:419122/1‑251 (MQ=9)                                                                                                    aggagAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGga  <  1:867870/250‑1 (MQ=2)                                                                                                    aggagAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAGGATATAGGAAATTCAAAAGAA‑CAAATGTTGCAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGga  <  1:220478/250‑1 (MQ=2)                                                                                                                                                                                                                                                       |                                                                                                        ATCAACAGGAATAATCTTAGATAAAGAATTTAAATATCTAAGAGATGAATGGAGAACAGCTTCTAAAGATTCAAATAAAATAAAAGCATTTGGAACTAGAGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAA‑CATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTGGACAATATCAGGA  >  CP028101/880537‑880885

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