Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,032 (G)6→7 coding (1002/1047 nt) wzzE → Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,0261.G97.1% 238.4 / ‑0.7 69coding (996/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (34/33);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

ACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTG  >  NC_000913/3969899‑3970148
                                                                                                                                |                                                                                                                           
acacCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑gggggggctga                                                                                                                 <  1:1406126/139‑6 (MQ=255)
    ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgg                                                                                                              >  2:1316447/1‑139 (MQ=255)
    ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgg                                                                                                              >  2:463068/1‑139 (MQ=255)
    ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgg                                                                                                              >  1:132133/1‑139 (MQ=255)
              ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGgtgt                                                                                                     >  1:1100051/1‑139 (MQ=255)
              ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTc                                                                                                    >  2:1119748/1‑139 (MQ=255)
              ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTc                                                                                                    >  2:647281/1‑139 (MQ=255)
              ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCAAGTCGTGCCTTCCTGATGATTATCCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTc                                                                                                    >  2:962471/1‑139 (MQ=255)
                tCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGc                                                                                                  <  2:75713/139‑1 (MQ=255)
                            tCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                      <  1:142789/139‑1 (MQ=255)
                               gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGc                                                                                                           >  1:1209518/1‑115 (MQ=255)
                               gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGc                                                                                                           <  2:1209518/115‑1 (MQ=255)
                               gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                   >  1:721646/1‑139 (MQ=255)
                               gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                   >  1:280076/1‑139 (MQ=255)
                                 gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                 <  1:962471/139‑1 (MQ=255)
                                  ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                <  1:1316447/139‑1 (MQ=255)
                                    ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              <  2:721646/139‑1 (MQ=255)
                                    ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              <  2:280076/139‑1 (MQ=255)
                                     tCAGACCTATTGCTATTTGCGTCCGCAGGAAGAACCGGTAAACCGCGATAGCCCACGTGGTGTCTTCCAAATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              <  2:1100051/139‑1 (MQ=255)
                                             aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑GGGGGGGCTATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                  >  2:1265094/1‑125 (MQ=255)
                                             aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑GGGGGGGCTATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                  <  1:1265094/125‑1 (MQ=255)
                                                      tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                            >  2:468102/1‑139 (MQ=255)
                                                      tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                            >  2:721841/1‑139 (MQ=255)
                                                        cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                          <  2:400472/139‑1 (MQ=255)
                                                            cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtg                                                                                                      <  1:1143630/91‑1 (MQ=255)
                                                            cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtg                                                                                                      >  2:1143630/1‑91 (MQ=255)
                                                               cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCa                                                   <  1:721841/139‑1 (MQ=255)
                                                                  aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgggg                                                                                                            <  1:156574/79‑1 (MQ=255)
                                                                  aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgggg                                                                                                            >  2:156574/1‑79 (MQ=255)
                                                                  aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgggg                                                                                                            >  2:852362/1‑79 (MQ=255)
                                                                  aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgggg                                                                                                            <  1:852362/79‑1 (MQ=255)
                                                                       ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTg                                                         <  2:850951/125‑1 (MQ=255)
                                                                       ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTg                                                         >  1:850951/1‑125 (MQ=255)
                                                                          gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgc                                        >  2:237959/1‑139 (MQ=255)
                                                                          gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgc                                        >  1:986664/1‑139 (MQ=255)
                                                                          gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgc                                        <  1:647281/139‑1 (MQ=255)
                                                                           tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              <  1:1248936/100‑1 (MQ=255)
                                                                           tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              >  2:1248936/1‑100 (MQ=255)
                                                                            aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                      <  1:237959/139‑1 (MQ=255)
                                                                              aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGTTGTCGTATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                    <  2:132133/139‑1 (MQ=255)
                                                                              aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                    <  1:468102/139‑1 (MQ=255)
                                                                                cgcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                    <  1:1068618/137‑1 (MQ=255)
                                                                                cgcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                    >  2:1068618/1‑137 (MQ=255)
                                                                                    aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtg                                                                                                      >  2:1183323/1‑67 (MQ=255)
                                                                                    aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtg                                                                                                      <  1:1183323/67‑1 (MQ=255)
                                                                                    aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                              <  1:1119748/139‑1 (MQ=255)
                                                                                     tAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGa                             >  1:214067/1‑139 (MQ=255)
                                                                                             gtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATc                     >  1:979576/1‑139 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGAt                                                                                                                 >  1:1353360/1‑43 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGAt                                                                                                                 <  2:1353360/43‑1 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                           <  2:1246405/113‑1 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                           >  1:1246405/1‑113 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgc                                     <  2:283255/119‑1 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgc                                     >  1:283255/1‑119 (MQ=255)
                                                                                                 tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                 <  1:463068/139‑1 (MQ=255)
                                                                                                           atgatTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                <  1:1016494/66‑1 (MQ=255)
                                                                                                           atgatTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                >  2:1016494/1‑66 (MQ=255)
                                                                                                                tATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg  >  2:728496/1‑139 (MQ=255)
                                                                                                                  tgcggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCg                                                                                                   <  2:1107345/37‑1 (MQ=255)
                                                                                                                  tgcggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCg                                                                                                   >  1:1107345/4‑40 (MQ=255)
                                                                                                                     ggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATTGATGTGAAAGTactgac    >  2:1440845/1‑132 (MQ=255)
                                                                                                                     ggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac    <  1:1440845/132‑1 (MQ=255)
                                                                                                                      gggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                              >  2:346996/1‑105 (MQ=255)
                                                                                                                      gggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                              <  1:346996/105‑1 (MQ=255)
                                                                                                                       ggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt                                <  1:1246869/102‑1 (MQ=255)
                                                                                                                       ggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt                                >  2:1246869/1‑102 (MQ=255)
                                                                                                                        gCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa                                                                            >  2:71828/1‑57 (MQ=255)
                                                                                                                        gCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa                                                                            <  1:71828/57‑1 (MQ=255)
                                                                                                                          aTTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                    >  2:271469/1‑63 (MQ=255)
                                                                                                                          aTTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                    <  1:271469/63‑1 (MQ=255)
                                                                                                                                |                                                                                                                           
ACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTG  >  NC_000913/3969899‑3970148

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