Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,095,255 T→A 5.9% K222M (AAG→ATG)  USA300HOU_RS10605 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,095,2550TA 5.9% 178.5 / 3.0 61K222M (AAG→ATG) USA300HOU_RS10605hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/22);  new base A (2/2);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TGTTTTGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGTTTTTCTTTTTG  >  USA300TCH1516_ALE/2095116‑2095389
                                                                                                                                           |                                                                                                                                      
tGTTTTGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCtt                                                                                                                                       >  2:634229/1‑141 (MQ=255)
   tttGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTcc                                                                                                                                    >  1:162897/1‑141 (MQ=255)
    ttGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTTATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTTACTTTGCGCTTCGTCCTCTGTCACAACTAGCTCTAAACGCTTTCCt                                                                                                                                   >  2:776/1‑141 (MQ=255)
    ttGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCt                                                                                                                                   >  2:92541/1‑141 (MQ=255)
    ttGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCt                                                                                                                                   >  1:295679/1‑141 (MQ=255)
           aaGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCAt                                                                                                                            <  2:50361/141‑1 (MQ=255)
           aaGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCATTCCTTTGTCAt                                                                                                                            >  2:398449/1‑141 (MQ=255)
           aaGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCATTCCTTTGTCAt                                                                                                                            >  2:398395/1‑141 (MQ=255)
                cTTAATTGTGTTTCTTCCATATTTTGATCATCTGATTGTGGTTCATTTATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCg                                                                                                                       <  2:158054/141‑1 (MQ=255)
                cTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCg                                                                                                                       <  1:287570/141‑1 (MQ=255)
                  tAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCggg                                                                                                                     <  1:324828/141‑1 (MQ=255)
                   aaTCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGt                                                                                                                    >  1:129511/1‑141 (MQ=255)
                      cGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTccc                                                                                                                 >  1:33189/1‑141 (MQ=255)
                                   tatTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCCCAGCAATTAATg                                                                                                    >  1:325560/1‑141 (MQ=255)
                                   tatTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATg                                                                                                    >  2:239549/1‑141 (MQ=255)
                                   tatTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGCCCCACAGCACTTAATg                                                                                                    >  2:101146/1‑141 (MQ=255)
                                     ttttGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGcg                                                                                                  >  1:207558/1‑140 (MQ=255)
                                         gATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGCATTGACTTTGCGCTTCTTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGttt                                                                                              <  2:74695/141‑1 (MQ=255)
                                             acatGATTGTGGTTCATATATCCCCCAAATATAGCGCGTAGGTAGGTTGAATTGACTTTGCGCTTCCTCCTCTGTCACACCTACCTCTAAACGCTTTCCTCTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGAt                                                                                          <  2:332209/138‑1 (MQ=255)
                                                       ggTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACata                                                                                <  2:33189/141‑1 (MQ=255)
                                                        gTTCCTATATCCCCCAAATATACCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACatat                                                                               <  1:31227/141‑1 (MQ=255)
                                                        gTTCATTTCTCCCCGAAATGTCGCTCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACatat                                                                               <  2:173495/141‑1 (MQ=255)
                                                        gTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGActtct                                                                               >  2:274055/1‑137 (MQ=255)
                                                            acttcgcccccAAATATCGCGCATACGTAGCTTGAATTGCATTGCCGCTCCGTCGTCTGTAACAACTAGCTCTACACGCTGTCCTTTGTTATTTTACGGCATCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAAt                                                                           <  2:15770/135‑1 (MQ=255)
                                                             tatatCCCCCAAAGATAGCGCNTNGGTNGGTTGAATTGACTTTGCGCTTCGTCATCTGTCCCAACTAGCTCTAAACGCTGTCCTTTGTCATTTTCGGGTCCCCCAGCAATTAATGCTGTTCTGATTTCTGACAGATCAATc                                                                          <  1:438145/141‑1 (MQ=255)
                                                              atatCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCt                                                                         >  1:682860/1‑141 (MQ=255)
                                                               tatCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCtt                                                                        >  2:473106/1‑141 (MQ=255)
                                                                     ccAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGt                                                                  <  1:179986/141‑1 (MQ=255)
                                                                         atatAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTCTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGCCATATCAATCTTCCTAGTTAAc                                                              >  2:552493/1‑141 (MQ=255)
                                                                           ataGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAAcgc                                                            >  2:457581/1‑139 (MQ=255)
                                                                            taGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCAGCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGAGTTCTGACATATCAATCTTCCTAGTTAACcca                                                           <  2:139962/141‑1 (MQ=255)
                                                                              gcgcATAGGTTGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACccaac                                                         <  2:483687/141‑1 (MQ=255)
                                                                               cgcATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACccaacc                                                        <  2:595160/141‑1 (MQ=255)
                                                                                gcATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGACATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACccaacca                                                       <  2:162897/141‑1 (MQ=255)
                                                                                  aTAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCaaa                                                     >  2:691356/1‑141 (MQ=255)
                                                                                      gtaggtTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCtt                                                 >  2:429523/1‑141 (MQ=255)
                                                                                      gtaggtTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCtt                                                 >  2:275582/1‑141 (MQ=255)
                                                                                             gAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCaga                                          >  1:64867/1‑139 (MQ=255)
                                                                                              aaTTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGAGTTCTGGCATATCAATCTTCCTAGTTAACCCAACCCAATCATTACCATAc                                         >  1:255059/1‑141 (MQ=255)
                                                                                              aaTTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGGTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATAc                                         >  1:632413/1‑141 (MQ=255)
                                                                                               aTTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCATTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACt                                        <  1:398449/141‑1 (MQ=255)
                                                                                               aTTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCATTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACt                                        <  1:398395/141‑1 (MQ=255)
                                                                                               aTTGACTTTGCGCTTCGTCATCTGTCACAAATAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGCTTTCTGACATATCAAACTTCCTAGATAACCAACCCAAATATTTTCCATACt                                        >  2:23056/1‑141 (MQ=255)
                                                                                                ttGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTCTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTc                                       <  1:18105/141‑1 (MQ=255)
                                                                                                 tGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCa                                      >  2:545617/1‑141 (MQ=255)
                                                                                                   aCTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTACCACAACCATACCTCTACCCTACTCAAt                                    >  2:344804/1‑141 (MQ=255)
                                                                                                    cTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTACCCCAACCAAATCTTTACCATACTCAAtt                                   >  1:495999/1‑141 (MQ=255)
                                                                                                     tttGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAAATACACCAAATCTTTACCATACTCAAttt                                  >  2:142324/1‑141 (MQ=255)
                                                                                                      ttGCGCTTCGTCATCTGTCACAACTAGCTCGAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTc                                 <  2:124175/141‑1 (MQ=255)
                                                                                                          gcTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATCCTCAATTTCtttc                             >  1:31769/1‑140 (MQ=255)
                                                                                                              cGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCtt                         >  1:133925/1‑141 (MQ=255)
                                                                                                                 cATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGa                      >  2:488234/1‑141 (MQ=255)
                                                                                                                   tCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATCCTCAATTTCTTTACCTTTGAAt                    >  2:290950/1‑141 (MQ=255)
                                                                                                                      gTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTGGAATaga                 >  2:245947/1‑139 (MQ=255)
                                                                                                                         acaACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTg              >  1:620380/1‑141 (MQ=255)
                                                                                                                            aCTAGCTCTAAACGCTTTCCCTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGttt           >  1:423079/1‑141 (MQ=255)
                                                                                                                             cTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGtttt          >  2:233032/1‑141 (MQ=255)
                                                                                                                               tgCTTTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGGTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCGACCAAATCTTTACCATACTCAATTTCTTCACCTTTGAATAAGCTGTTTTTc        <  2:494049/140‑1 (MQ=255)
                                                                                                                               aGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTGAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGTTTTTc        <  1:20967/141‑1 (MQ=255)
                                                                                                                                  tctAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGTTTTTCttt     <  2:133925/141‑1 (MQ=255)
                                                                                                                                   ctAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGACTCCGCTGTCTTTCtttt    >  2:106392/1‑141 (MQ=255)
                                                                                                                                     aaaCGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGTTTTTCTTTTTg  <  2:682860/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
TGTTTTGGCTAAAGAACTTAATCGTGTTTCATTCATATTTTGATCATCTGATTGTGGTTCATATATCCCCCAAATATAGCGCATAGGTAGGTTGAATTGACTTTGCGCTTCGTCATCTGTCACAACTAGCTCTAAACGCTTTCCTTTGTCATTTTCGGGTCCCACAGCAATTAATGCTGTTTTGATTTCTGACATATCAATCTTCCTAGTTAACCCAACCAAATCTTTACCATACTCAATTTCTTTACCTTTGAATAAGCTGTTTTTCTTTTTG  >  USA300TCH1516_ALE/2095116‑2095389

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