Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE707,8790AT23.8% 128.7 / 46.6 63L11I (TTA→ATA) znuC_1High‑affinity zinc uptake system ATP‑binding protein ZnuC
Reads supporting (aligned to +/- strand):  ref base A (28/20);  new base T (8/7);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTAT  >  USA300TCH1516_ALE/707742‑708017
                                                                                                                                         |                                                                                                                                          
aCCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTaaaaa                                                                                                                                         <  1:663052/141‑1 (MQ=255)
 ccGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAAc                                                                                                                                        >  2:293774/1‑141 (MQ=255)
 ccGGTAGCATTAAATTCACCAATTAAAGAATTGATAAGGGAAGATTTACCAGCAGCATTCGGGCCCATGATACTAAGTATTTCGGCGCTTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAAc                                                                                                                                        >  2:401048/1‑141 (MQ=255)
    gTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGa                                                                                                                                     >  1:396454/1‑141 (MQ=255)
    gTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGa                                                                                                                                     <  1:285376/141‑1 (MQ=255)
          ttAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTaaa                                                                                                                               >  1:911614/1‑141 (MQ=255)
                tCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCtttt                                                                                                                         <  2:420591/141‑1 (MQ=255)
                  aCCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGt                                                                                                                       >  2:199382/1‑141 (MQ=255)
                    cAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGttt                                                                                                                     <  1:740209/141‑1 (MQ=255)
                         aaaGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAAc                                                                                                                >  1:235278/1‑141 (MQ=255)
                         aaaGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAAc                                                                                                                >  1:74851/1‑141 (MQ=255)
                           aGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACGGATTTAAATCTTTTGTTTCTAACaa                                                                                                              >  2:257885/1‑141 (MQ=255)
                               ttGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGt                                                                                                          <  2:132098/141‑1 (MQ=255)
                                tGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGtt                                                                                                         >  2:312713/1‑141 (MQ=255)
                                   tAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtat                                                                                                      >  1:390100/1‑141 (MQ=255)
                                   tAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtat                                                                                                      >  2:155368/1‑141 (MQ=255)
                                         aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcctcc                                                                                                <  1:731775/141‑1 (MQ=255)
                                         aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctcc                                                                                                <  2:341785/141‑1 (MQ=255)
                                          aGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcctcct                                                                                               >  1:679755/1‑141 (MQ=255)
                                          aGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctcct                                                                                               >  1:622114/1‑141 (MQ=255)
                                          aGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTGTCTAACAAGCGTTTATAcctcct                                                                                               >  2:24179/1‑141 (MQ=255)
                                              ttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAAtt                                                                                           <  1:176372/141‑1 (MQ=255)
                                              ttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAAtt                                                                                           <  1:229812/141‑1 (MQ=255)
                                                    gCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGt                                                                                     <  2:253284/141‑1 (MQ=255)
                                                         aTTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAgg                                                                                >  2:342064/1‑141 (MQ=255)
                                                             gggCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACGTGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTaa                                                                            >  1:554427/1‑141 (MQ=255)
                                                                   aTGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCtaat                                                                      >  1:185317/1‑141 (MQ=255)
                                                                        aCCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATATTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTCACCTaattaatt                                                                 >  2:187557/1‑141 (MQ=255)
                                                                              tattTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTataa                                                           <  1:713141/141‑1 (MQ=255)
                                                                                 ttCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTAtaataa                                                        >  1:722016/1‑141 (MQ=255)
                                                                                 ttCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTAtaataa                                                        >  2:172450/1‑141 (MQ=255)
                                                                                   cgccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAAc                                                      <  1:464742/141‑1 (MQ=255)
                                                                                    gccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACATACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACt                                                     >  1:471038/1‑141 (MQ=255)
                                                                                     ccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTg                                                    <  1:549836/141‑1 (MQ=255)
                                                                                       ccgTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAANAGTTTAGGCTAACCTAATTAATTGTATAATAAACTggg                                                  <  1:161443/140‑3 (MQ=255)
                                                                                            cTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAatat                                             <  1:88934/141‑1 (MQ=255)
                                                                                             tGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAACAACAGATTTACATATTTTGTTTCTAACAAGAGCACACACCTCCTCATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATAtt                                            <  1:492558/141‑1 (MQ=255)
                                                                                              ggTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATAttt                                           >  1:100354/1‑141 (MQ=255)
                                                                                                   cGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCAt                                      >  2:140418/1‑141 (MQ=255)
                                                                                                   cGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCAt                                      >  2:140417/1‑141 (MQ=255)
                                                                                                    gATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATg                                     <  1:342539/141‑1 (MQ=255)
                                                                                                    gATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATg                                     <  1:563962/141‑1 (MQ=255)
                                                                                                        aGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTAGTTTCTAACAAGCGTTTATACCGCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCaa                                 >  1:193288/1‑141 (MQ=255)
                                                                                                         ggAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                <  1:90405/141‑1 (MQ=255)
                                                                                                         ggAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                <  1:456284/141‑1 (MQ=255)
                                                                                                           aaaTGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTa                              <  1:782276/141‑1 (MQ=255)
                                                                                                                  tttAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                       >  2:324858/1‑141 (MQ=255)
                                                                                                                  tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                       >  2:244345/1‑141 (MQ=255)
                                                                                                                  tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                       >  1:318094/1‑141 (MQ=255)
                                                                                                                  tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                       >  1:23538/1‑141 (MQ=255)
                                                                                                                   ttAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                      <  1:236369/141‑1 (MQ=255)
                                                                                                                   ttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                      >  1:147266/1‑141 (MQ=255)
                                                                                                                   ttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTATAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATAGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                      >  1:432073/1‑141 (MQ=255)
                                                                                                                         aCATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataa                >  1:679083/1‑141 (MQ=255)
                                                                                                                            tGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTTCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCACGTCAAGTAAATTCGTGatataatat             <  1:110527/141‑1 (MQ=255)
                                                                                                                              cTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAg           >  1:597677/1‑141 (MQ=255)
                                                                                                                                 attaCCTAAAAACAGATTTAAATCGTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACa        <  1:88666/141‑1 (MQ=255)
                                                                                                                                  ttaCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACaa       <  2:294902/141‑1 (MQ=255)
                                                                                                                                   taCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAat      <  1:703930/141‑1 (MQ=255)
                                                                                                                                    aCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatg     >  2:225696/1‑141 (MQ=255)
                                                                                                                                    aCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATATATCATGTCAAGTAAATTCGTGATATAATATAGACAatg     >  1:257240/1‑141 (MQ=255)
                                                                                                                                     ccTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgt    >  2:319593/1‑141 (MQ=255)
                                                                                                                                     ccTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgt    <  1:431421/141‑1 (MQ=255)
                                                                                                                                       tAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgtat  <  2:257982/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
ACCGGTAGCATTAAATTCACCAATTAAAGACTTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTAT  >  USA300TCH1516_ALE/707742‑708017

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