Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,454,626	A→T	L48L (CTA→CTT)	USA300HOU_2315 →	hypothetical membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,454,626	0	A	T	100.0%	201.5 / NA	60	L48L (CTA→CTT)	USA300HOU_2315	hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (33/27);  total (33/27)

ATGAAGAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTCTAATTATT  >  CP000730/2454483‑2454770                                                                                                                                                |                                                                                                                                                  aTGAAGAACAATAAGATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAc                                                                                                                                               <  2:67273/148‑1 (MQ=255)  tGAAGAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAcc                                                                                                                                              >  1:630639/1‑148 (MQ=255)   gaagaaCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCa                                                                                                                                             <  1:285696/148‑1 (MQ=255)    aagaaCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGTACCAt                                                                                                                                            >  1:169049/1‑148 (MQ=255)    aagaaCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCAt                                                                                                                                            <  1:307483/148‑1 (MQ=255)    aagaaCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCACTATGCTCAGATATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCAt                                                                                                                                            <  1:147484/148‑1 (MQ=255)           aTAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgc                                                                                                                                     <  1:580757/148‑1 (MQ=255)            tAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgct                                                                                                                                    <  1:225028/148‑1 (MQ=255)               aaTTTCTGGTTTTCAATGGGCAATGACGATTTTAGGAGCGTTTTTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATGTGAAGTTACAGATCTTGCACCATTAATTGCTGCa                                                                                                                                 >  1:150119/1‑148 (MQ=255)                     tGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCAGTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑Att                                                                                                                          >  2:260940/1‑148 (MQ=255)                      ggTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑Attt                                                                                                                         >  1:440581/1‑148 (MQ=255)                      ggTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTATCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTAGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑Attt                                                                                                                         >  1:405907/1‑148 (MQ=255)                      ggTTTTCAATGGGCAATGACGATTTTCGGCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATACTTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑Attt                                                                                                                         <  1:275197/148‑1 (MQ=255)                            cAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATAc                                                                                                                   >  2:550515/1‑148 (MQ=255)                             aaTGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACt                                                                                                                  >  2:244129/1‑148 (MQ=255)                                gggCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGt                                                                                                               <  1:586681/148‑1 (MQ=255)                                        aCGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATAGTCGTTTTCAAat                                                                                                       >  2:153521/1‑148 (MQ=255)                                         cgctttctgTCTTCTTCGTCACTACCATGGCGTTATCCATTATGCTCTGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTTTATTTGAAGTTACAGATCTTGCACCATTACTTGCTGCAATC‑ATTTGTTTACTCGTTTTCAAgta                                                                                                      <  2:469693/140‑1 (MQ=255)                                         cGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACGATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAata                                                                                                      >  2:320564/1‑148 (MQ=255)                                            ttttCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATaa                                                                                                   <  1:90439/148‑1 (MQ=255)                                            ttttCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATA‑ATTTGTATACTCGTTTTCAAATATaa                                                                                                   >  2:11418/1‑148 (MQ=255)                                                  tcttctTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGt                                                                                             <  1:445642/148‑1 (MQ=255)                                                   cttcttTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTc                                                                                            >  1:390762/1‑148 (MQ=255)                                                   cttcctGGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTc                                                                                            <  2:480151/148‑1 (MQ=255)                                                    ttcttTGGCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTcc                                                                                           <  1:388864/148‑1 (MQ=255)                                                          gTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCTCCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTg                                                                                     >  2:224880/1‑148 (MQ=255)                                                          gTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTg                                                                                     >  2:529486/1‑148 (MQ=255)                                                               tACAATGGCGTTATCCATTATGCTCAGAGGTTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGACCTTGCACCATTAATTGCTGCAAGC‑ATTTGGATACTCGTTTTCAAATATAAAAAGGTCCAACTTgagggt                                                                                >  1:335093/1‑143 (MQ=255)                                                               tACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGt                                                                                >  2:93044/1‑148 (MQ=255)                                                                      tcGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTCACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAATAGGTCCAACTTGCAGGTTTAAAAt                                                                         <  1:311625/147‑1 (MQ=255)                                                                            tCCATTAAGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATAGGTG‑ATTCTTTTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCaa                                                                   >  1:458465/1‑148 (MQ=255)                                                                                 tATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGc                                                              <  1:529557/148‑1 (MQ=255)                                                                                   tGCTCAGAGATTTCCAGTCTATAATTGGTCTGCCTCGCTTTATATTTGAAGTTACAGATCTTGCACCATAAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCt                                                            >  2:446909/1‑148 (MQ=255)                                                                                   tGCTCAGAGATTTCCAGACTATAATTGGTATCAAACACTTTATATTTGACGTTACAGATCTTGAAACATTAATTGCTGCAATC‑ATTTCTATACTCGTTTTCAAATATAAAAAGGTCCAAATTGCAGGTTTAAAATTCTCAATCaactt                                                            >  2:311119/1‑144 (MQ=255)                                                                                      tCAGACATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTTCAGGTTTAAAATTCTCAATCAGCCTGaa                                                         <  1:37323/148‑1 (MQ=255)                                                                                       cAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTATTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGaaa                                                        <  2:548023/148‑1 (MQ=255)                                                                                                   gTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGt                                            <  1:476832/148‑1 (MQ=255)                                                                                                   gTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGt                                            <  2:390762/148‑1 (MQ=255)                                                                                                    tCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTc                                           <  1:260940/148‑1 (MQ=255)                                                                                                        taATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTAtt                                       >  1:393443/1‑148 (MQ=255)                                                                                                            tGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTa                                   >  1:483652/1‑148 (MQ=255)                                                                                                              gtgtCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCTAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGc                                 <  1:416723/148‑1 (MQ=255)                                                                                                               tgtCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCt                                >  2:525243/1‑148 (MQ=255)                                                                                                                gtCAAACAATTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCtt                               <  1:208217/148‑1 (MQ=255)                                                                                                                   aaaCACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCACTC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTaa                            >  2:423239/1‑148 (MQ=255)                                                                                                                      cacTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAAttt                         >  1:189818/1‑148 (MQ=255)                                                                                                                       acTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAAtttt                        <  1:25684/148‑1 (MQ=255)                                                                                                                        cTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTATCTTTAATTTTa                       <  2:169049/148‑1 (MQ=255)                                                                                                                         tttATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTAc                      <  2:440581/148‑1 (MQ=255)                                                                                                                             tatTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTCGCTTTAATTTTACCttt                  >  2:560354/1‑148 (MQ=255)                                                                                                                             tatTTGAAGTTACAGATCTTGCACCATTAATTGATGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGCACGTCTATTGCTAGTCCTAAttttccatat                  >  2:560760/1‑142 (MQ=255)                                                                                                                              atTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTa                 <  1:287982/148‑1 (MQ=255)                                                                                                                              atTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTa                 >  2:198234/1‑148 (MQ=255)                                                                                                                              atTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTa                 <  1:79099/148‑1 (MQ=255)                                                                                                                                ttGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTaagtcagacgataat               >  2:128127/1‑135 (MQ=255)                                                                                                                                ttGAAGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAAGTTTACCTTTAat               >  2:452893/1‑148 (MQ=255)                                                                                                                                   aaGTTACAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAattat            >  1:263002/1‑148 (MQ=255)                                                                                                                                         cAGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTCTAat      >  1:244112/1‑148 (MQ=255)                                                                                                                                          aGATCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCATGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATACTAatt     >  2:61746/1‑148 (MQ=255)                                                                                                                                             tCTTGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATCTACAAAGGAACAAATTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTCTAatgttt  >  1:530163/1‑144 (MQ=255)                                                                                                                                                |                                                                                                                                                  ATGAAGAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATC‑ATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTCTAATTATT  >  CP000730/2454483‑2454770

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