Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,196,508	C→A	S338Y (TCT→TAT)	ftsZ →	Cell division protein FtsZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,196,508	0	C	A	100.0%	194.7 / NA	61	S338Y (TCT→TAT)	ftsZ	Cell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (36/25);  total (36/25)

CGTTAATATGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTTCAA  >  USA300TCH1516_ALE/1196362‑1196651                                                                                                                                                   |                                                                                                                                                cGTTAATATGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGATCAAGCGTAAATACTtata                                                                                                                                               <  1:807546/149‑1 (MQ=255)       tatGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACCACCCAATAGCGCATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATg                                                                                                                                         <  1:488204/149‑1 (MQ=255)       tatGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCaagg                                                                                                                                         >  1:628590/1‑147 (MQ=255)       tatGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGAACTGGATTCGGAACAAGCGTAAATACTTATAGCAATg                                                                                                                                         >  1:95561/1‑149 (MQ=255)          gATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAAGACTTATAGCAATGCaa                                                                                                                                      >  2:670660/1‑149 (MQ=255)             tttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACtt                                                                                                                                   <  1:524192/149‑1 (MQ=255)              ttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACtta                                                                                                                                  >  1:80128/1‑148 (MQ=255)                 ggTACAGTTATTAATCCTGACTTACAAGATGAGATTGTTGTAAAAGTTATTGCAACTGGTTTTGAGGACAACCCACCATCACATGGTCGTAAATCTGGTAGCAGTGGATTCGGAACAAGCGTAAATACTTATAGCAACCGAACTTCTaa                                                                                                                               >  1:161821/1‑149 (MQ=255)                   tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAg                                                                                                                             <  1:978716/149‑1 (MQ=255)                       gTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATACGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGa                                                                                                                         >  1:208218/1‑149 (MQ=255)                        ttattaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                                                                                        >  2:669815/1‑149 (MQ=255)                        ttattaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                                                                                        >  1:616184/1‑149 (MQ=255)                             aaTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAGTACTTATAGCAATGCAACTTCTAAAGATGAATCAtt                                                                                                                   >  1:992570/1‑149 (MQ=255)                              aTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTc                                                                                                                  >  1:617762/1‑149 (MQ=255)                              aTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTc                                                                                                                  <  1:328931/149‑1 (MQ=255)                                ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                                >  1:345585/1‑149 (MQ=255)                                  tGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                                                                              <  1:685870/149‑1 (MQ=255)                                   gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACATACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                                                                             >  2:1015288/1‑149 (MQ=255)                                    aaTTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTGGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCa                                                                                                            >  1:857146/1‑149 (MQ=255)                                    aaTTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCa                                                                                                            >  1:631040/1‑149 (MQ=255)                                         cAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAtaa                                                                                                       >  2:275837/1‑147 (MQ=255)                                               tagaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATAATTCACTTCAACTTCATCaaa                                                                                                 >  2:553082/2‑149 (MQ=255)                                                agaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAAt                                                                                                >  1:108797/1‑149 (MQ=255)                                                   ttgttgTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGca                                                                                             >  2:496885/1‑149 (MQ=255)                                                     gttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGcaca                                                                                           <  2:866417/149‑1 (MQ=255)                                                      ttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACaa                                                                                          >  1:945068/1‑149 (MQ=255)                                                      ttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACaa                                                                                          >  2:1031053/1‑149 (MQ=255)                                                         aaaCAGTTATGGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCa                                                                                       <  2:1091141/148‑1 (MQ=255)                                                            cAGTTCGTGCAACTGGTTTTGCGGACGGTGCAACAGCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGGAACt                                                                                    <  2:779499/149‑1 (MQ=255)                                                            cAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAACTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCCCTTCAAATTCCTCCAATGCACAAGCaaat                                                                                    >  1:841106/1‑147 (MQ=255)                                                            cAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACt                                                                                    >  2:272470/1‑149 (MQ=255)                                                            cAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACt                                                                                    >  1:841105/1‑149 (MQ=255)                                                                tactGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATa                                                                                <  1:705973/146‑1 (MQ=255)                                                                  ttGTAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACGTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgt                                                                              <  2:673657/149‑1 (MQ=255)                                                                  ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCCAATGCACAAGCAACTGATAga                                                                              >  2:908320/1‑148 (MQ=255)                                                                   tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                             <  2:177842/149‑1 (MQ=255)                                                                   tGCAACTGGTTTTGAAGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                             <  2:720948/149‑1 (MQ=255)                                                                   ggCAACGGGTTTTGATGACAAACCAACAGCACATGGGCGGAAGGCGGGGAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTGCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                             <  1:756950/148‑1 (MQ=255)                                                                     cAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTa                                                                           <  1:421045/149‑1 (MQ=255)                                                                     aaaCGGGTTTTGATGACAAACCAACAACACATGGTCGTAAATCTGGTAGCCCTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAATTCACCGATATCAAATGCACAAGCAACTGATAGTGTa                                                                           <  1:269702/148‑1 (MQ=255)                                                                          ggTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCGTCAAATGCACAAGCAACTGATAGTGTAAGTGa                                                                      <  1:712243/149‑1 (MQ=255)                                                                           ggtttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGaa                                                                     <  1:802662/147‑1 (MQ=255)                                                                            ttttGATGACAAACCAACATCACCAGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAGGTGaaa                                                                    >  1:222514/1‑149 (MQ=255)                                                                            ttttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAACGCACAAGCAACTGATAGTGTAAGTGaaa                                                                    <  1:679457/149‑1 (MQ=255)                                                                              ttGATGAAAAACCAACATCACATGGTCGTAAATCCGGTAGCACCGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGa                                                                  <  1:277402/149‑1 (MQ=255)                                                                                       aaCCAACATCACATGGTCGTAACACTGGTAGCTGGGGATTCGGCCGGTCCTGGTGTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                         <  1:806024/149‑1 (MQ=255)                                                                                               tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAaaga                                                 <  1:351969/149‑1 (MQ=255)                                                                                                  caTGGTCGTAAATCTGGTAGCACTGGATTCGGAACACGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                              >  2:565583/1‑149 (MQ=255)                                                                                                   aTGGTCGTAAATCTGGTAGCACTTGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgat                                             <  2:382758/149‑1 (MQ=255)                                                                                                     ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACCTACAACTAAAGAAgatga                                           <  1:807535/149‑1 (MQ=255)                                                                                                     ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatga                                           >  1:621543/1‑149 (MQ=255)                                                                                                         gTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATAtt                                       >  2:578002/1‑149 (MQ=255)                                                                                                          tAAATCTGGTAGCACTGGATTCGGAACAAGCGTAACTTCTGATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTc                                      >  1:1101031/1‑149 (MQ=255)                                                                                                                   tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAATAAGATGATAATCCTAGCTTcc                             >  2:962218/1‑148 (MQ=255)                                                                                                                    aGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCAt                            <  2:124909/149‑1 (MQ=255)                                                                                                                           gATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAAt                     >  1:288186/1‑149 (MQ=255)                                                                                                                              tCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATaga                  >  1:216469/1‑149 (MQ=255)                                                                                                                                   aCAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCGAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTGGAAATAGagaaga             >  1:618459/1‑149 (MQ=255)                                                                                                                                          tAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGACGATGCTATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGt      <  2:1101031/149‑1 (MQ=255)                                                                                                                                             aTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTTCa   >  1:73793/1‑149 (MQ=255)                                                                                                                                              tACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAGGTAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTGCaa  >  1:31432/1‑149 (MQ=255)                                                                                                                                              tACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAATACGTTCaa  >  1:503480/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                CGTTAATATGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTTCAA  >  USA300TCH1516_ALE/1196362‑1196651

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