Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG45.6% 82.9 / 19.8 57N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (5/26);  new base G (26/0);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 2.63e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGC  >  NC_000913/2677899‑2678186
                                                                                                                                              |                                                                                                                                                    
tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGcc                                                                                                                                              <  1:733423/149‑1 (MQ=255)
     gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCCGAGGGTTTTCGCCCCTGTc                                                                                                                                        >  1:15718/1‑150 (MQ=255)
      ccGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAg                                                                                                                                      <  2:547852/151‑1 (MQ=255)
       cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGGCCGGGGGAGAGGGGTTTCGCCACTCTCAGc                                                                                                                                     >  2:325417/1‑151 (MQ=255)
       cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCAGGGGGTTTTCGcccc                                                                                                                                            >  2:267182/1‑142 (MQ=255)
        gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAGCGTCACCGCCCAGCGCCTGTAAGCAAGGTGGGGCATCCTGCGCCTTTTGCAGGTAATCAATCAGGCCGTTGCGGCGGGGCCACTGACCGAGGCATAGGTTTTCCGCCACTGTCAg                                                                                                                                      <  2:858648/149‑1 (MQ=255)
         ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGGTTTTCGCCCCCGTCAGccc                                                                                                                                   >  2:61469/1‑151 (MQ=255)
            ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGCGGCCGAGGGTTTTCGCCCCTGTCAGCCCtt                                                                                                                                 >  2:868529/1‑150 (MQ=255)
                cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTccacc                                                                                                                            <  1:471925/151‑1 (MQ=255)
                 gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTCACCCCGGCAGAGGGTTTCCCCCACTGTCAGCCCTcccccc                                                                                                                            >  1:368473/1‑150 (MQ=255)
                   tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAg                                                                                                                          <  1:447288/150‑1 (MQ=255)
                    gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCCCCAGAc                                                                                                                        >  1:455909/1‑151 (MQ=255)
                         aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAAc                                                                                                                   <  2:304405/151‑1 (MQ=255)
                          aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGGAGAGGGTTTCCGCCCCCGTCCGCCCCTCCCCCCGACtt                                                                                                                      >  1:676441/1‑147 (MQ=255)
                           gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAACGACCGAGGCAGGGGGTTTTCGCCACCGTCAGCCCTTTCACCAGACttta                                                                                                                    >  1:781140/1‑146 (MQ=255)
                              gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCCCCACTGTCCGCCCCTCCCCCCGACTTAACTcc                                                                                                                >  1:221936/1‑149 (MQ=255)
                              gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGGAAAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCCGACtttaaccc                                                                                                                >  2:539919/1‑143 (MQ=255)
                                  aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCCACCGGGCGGGGGGTTTTTGCCCCTGTCAGCCCCTCCCCCCGACTTCACTCCTGAt                                                                                                            >  2:191234/1‑149 (MQ=255)
                                           gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGGCCGGGGGGGAGGGTTTCCGCCACTGTCAGCCCCTTCACCCGACTTAAATCCCGATaaaaccc                                                                                                     >  2:171352/1‑143 (MQ=255)
                                             cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAc                                                                                              <  1:299190/151‑1 (MQ=255)
                                                cGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAcccc                                                                                           <  2:392284/151‑1 (MQ=255)
                                                  cccAGCGCCTGTAAGCAACGTTGGTCATCCTGCCCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCa                                                                                         <  2:205846/151‑1 (MQ=255)
                                                  cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAccccc                                                                                          <  1:841659/150‑1 (MQ=255)
                                                      gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGtt                                                                                      <  1:212645/150‑1 (MQ=255)
                                                      gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                     <  2:249984/151‑1 (MQ=255)
                                                        gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGGAGAGGGTTTTCGCCACCGTCAGCCCCTCCACCCGACTTAACTCCTGATAAAACGCCCGCCCCCCCCGtt                                                                                      >  1:123243/1‑148 (MQ=255)
                                                        gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGAACGAGGGAGAGGGTTTCCTCCACCGTCATCCCCTCCCCCAGACTTTACTCCCGCTAAACCCCGCGCGccccccct                                                                                       >  2:83970/1‑144 (MQ=255)
                                                         ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCCCCCCTGTCAGCCCCTCCACCAAACTTTACTCCTGATAAACCGCCCGCccccccccct                                                                                      >  1:769418/1‑143 (MQ=255)
                                                              aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGGGGGAGAGGGTTTCCGCCCCCGTCCGCCCCTCCCCCACACTTAAAACCTCATAAAACCCGCGCCCCCCCCGCTCAGCgg                                                                                >  2:265413/1‑148 (MQ=255)
                                                               aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGGGGGAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCACACTTTACTCCCGATAAACCGCGCGCCacccc                                                                                          >  2:715282/1‑137 (MQ=255)
                                                                 cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGGTTTTCGCCACCGTCAGCCCTTCCCCCCGACTTAAATCCTGATAAACCCCGCGCACCCCCCGTTCCGCGGCAc                                                                             >  1:65570/1‑148 (MQ=255)
                                                                   aCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGGCCGGGGGAGAGGTTTTTCGCCCCCGTCAGGCCTTCCACCCGGCTTTATACCTTATTAACCCCGCGCACCCCCAGCTCAGCGGccccg                                                                           >  2:440430/1‑144 (MQ=255)
                                                                     gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGGCGCAGGTTTTCGCCCCCGTTCACCCTTCCACCCCACCTAACACCTTAAAAAACACGCCCCACCCCCACTTctccgga                                                                               >  2:586042/1‑137 (MQ=255)
                                                                       tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCCCCCCTTTCAGCCCCTCCCCCAGACTTAACTCCTGATAAACCCCGCGCACCCCCCCTTCAGCGGCCCGGCg                                                                         >  2:371111/1‑146 (MQ=255)
                                                                       tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGAGGGTTTTCCCCCCCGTCAGCCCTTCCACCAGACCTAACTCCTGATAAACCCCGCGCACCCCCCGTTCAGCGGCACGGCGAGTc                                                                     >  1:45667/1‑150 (MQ=255)
                                                                                 gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc                                                          <  2:578783/151‑1 (MQ=255)
                                                                                     cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTCCGCCACTGTCAGCCCTTCCACCAGACCTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCGTAGGtt                                                           >  1:133115/1‑146 (MQ=255)
                                                                                      aTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACCGTCAGCCCTTCCACCAGACTTAATTCCTGATAAACCGCACGAACCCCCAGTTCAGCGGCACGGCGAGTCAACGTAGCTTCGTCAc                                                     <  1:779481/151‑1 (MQ=255)
                                                                                             ggTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCCAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCCGACTTTACTCCCGATAAAACCCGCGCACCCCCCGTTCAGCGGCCCGGCGAGTCAGCGGAGCTTCGTCACCCTcc                                                >  2:188036/1‑149 (MQ=255)
                                                                                               tAATCAATCATGCCGTTGCGGCGGGGCCACTCACCGAGGGAGAGGGTTTTCCCCCCTGTCAGCCCCTCCACCCAACTTAACTCCCCATAAAACCCGCGCCCCCCCCGTTCAGCGGCACGGCGGGTCAGCGGAGCTTCGTCACCtctccct                                             >  2:528729/1‑144 (MQ=255)
                                                                                                 atcaatcaTGCCGTTGCGGCGGGGCCCCTCTCCCGGGCAGAGGGTTTTCCCCCCTGTCAGCCCTTCCACCAGAATTTACTCCTTATAAAAAGCGCGCACCCCCCCGTCACCGGCACAGCGAGGCAGCGGAGCCTCGGCACCCTTCCGTcgc                                          >  1:440677/1‑151 (MQ=255)
                                                                                                    aatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCACGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGt                                        <  2:838849/150‑1 (MQ=255)
                                                                                                        aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCCCGCGCACCCCCCGTTCAGCGGCCCGGCGAGGCAGCGTAGCTTCGTCACCCTCCAGTCGCGTCTCGc                                   >  1:333854/1‑151 (MQ=255)
                                                                                                          gCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGcc                                  <  1:325469/150‑1 (MQ=255)
                                                                                                          gCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCa                                 <  1:835588/151‑1 (MQ=255)
                                                                                                                  ggcggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                         <  1:392235/151‑1 (MQ=255)
                                                                                                                     ggggCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                       <  1:699433/150‑1 (MQ=255)
                                                                                                                     ggggCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                      <  2:63008/151‑1 (MQ=255)
                                                                                                                          ccACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCg                  <  2:369648/151‑1 (MQ=255)
                                                                                                                               gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACAGCTAt              <  1:543317/150‑1 (MQ=255)
                                                                                                                                   gAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg         >  1:684658/1‑151 (MQ=255)
                                                                                                                                      ggCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc        <  2:4907/150‑1 (MQ=255)
                                                                                                                                       gCAGCGGTTGTCCCACACTGTCAGCCCTTCCACCAGACTTAACTCCTGATACACCGCGCGCACCCCCAGTTCAGCGGCACGGGGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTGGCCAATCCAGATATCACCGCTATCCGGGCGt      <  2:205572/151‑1 (MQ=255)
                                                                                                                                       gCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt      <  2:551851/151‑1 (MQ=255)
                                                                                                                                         agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc    >  1:714439/1‑151 (MQ=255)
                                                                                                                                          gagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg   >  2:645268/1‑151 (MQ=255)
                                                                                                                                            ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGc  <  2:794989/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGC  >  NC_000913/2677899‑2678186

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