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.G98.3% 206.7 / ‑5.6 59V332V (GTC→GTG) wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (27/31);  minor base A (0/1);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA  >  NC_000913/3969888‑3970157
                                                                                                                                           |                                                                                                                                    
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                       >  2:1393547/1‑139 (MQ=255)
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                       >  1:241092/1‑139 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggg                                                                                                                              <  1:987759/139‑2 (MQ=255)
               ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCgg                                                                                                                       >  1:1604765/1‑139 (MQ=255)
                   aaTGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCt                                                                                                                   >  1:356246/1‑139 (MQ=255)
                        tGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtgt                                                                                                              <  1:1504243/139‑1 (MQ=255)
                            ccAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCa                                                                                                          <  1:871608/139‑1 (MQ=255)
                             cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCAGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGc                                                                                                                    <  1:335016/128‑1 (MQ=255)
                             cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCAGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGc                                                                                                                    >  2:335016/1‑128 (MQ=255)
                                cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTa                                                                                                       >  1:835464/1‑138 (MQ=255)
                                cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTa                                                                                                       <  2:835464/138‑1 (MQ=255)
                                     gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                                                 >  1:746952/1‑139 (MQ=255)
                                        ccGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                               <  1:1393547/139‑1 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGACGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                            >  2:226938/1‑139 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgt                                                                                                                <  2:368121/118‑1 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgt                                                                                                                >  1:368121/1‑118 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          <  2:1604765/139‑1 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          <  1:698767/139‑1 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          <  2:1271746/139‑1 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCCGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCA‑GGGGGCCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                          <  1:1133086/139‑1 (MQ=255)
                                             ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         <  1:162054/139‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggg                                                                                                                              <  2:272268/89‑2 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggg                                                                                                                              >  1:272268/1‑88 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑g                                                                                                                                    <  2:1582836/83‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑g                                                                                                                                    >  1:1582836/1‑83 (MQ=255)
                                                            cTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                       <  2:79815/126‑1 (MQ=255)
                                                            cTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                       >  1:79815/1‑126 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAgc                                                               >  1:1089611/1‑139 (MQ=255)
                                                                         ccGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑gggggggc                                                                                                                             <  2:107385/74‑3 (MQ=255)
                                                                         ccGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑gggggggc                                                                                                                             >  1:107385/1‑72 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAg                                                         >  1:854104/1‑139 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAg                                                         >  1:170767/1‑139 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                    >  2:35183/1‑139 (MQ=255)
                                                                                    ggTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcg                                                  >  1:683834/1‑139 (MQ=255)
                                                                                      tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         <  1:867046/98‑1 (MQ=255)
                                                                                      tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         >  2:867046/1‑98 (MQ=255)
                                                                                               aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                               <  2:858179/99‑1 (MQ=255)
                                                                                               aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                               >  1:858179/1‑99 (MQ=255)
                                                                                               aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                                       >  1:125193/1‑139 (MQ=255)
                                                                                                      aCGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa                                <  2:241092/139‑1 (MQ=255)
                                                                                                        gtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa                                >  2:298985/1‑137 (MQ=255)
                                                                                                        gtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa                                <  1:298985/137‑1 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         >  2:732673/1‑75 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         <  1:732673/75‑1 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGTTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                        >  1:103766/1‑139 (MQ=255)
                                                                                                                    tgatgaTTATGCGGGGCATTGCCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                  <  2:1089611/139‑1 (MQ=255)
                                                                                                                     gatgatTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTac                 <  2:746952/139‑1 (MQ=255)
                                                                                                                        gatTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga              <  2:125193/139‑1 (MQ=255)
                                                                                                                        gatTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga              <  1:226938/139‑1 (MQ=255)
                                                                                                                         atTATGCGGGGCATTGTCG‑GGGGGGTTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac             <  2:103766/139‑1 (MQ=255)
                                                                                                                         atTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac             <  2:170767/139‑1 (MQ=255)
                                                                                                                         atTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac             <  2:854104/139‑1 (MQ=255)
                                                                                                                          ttATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact            <  1:777263/139‑1 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtgt                                                                                                              <  1:778136/40‑1 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGgtgt                                                                                                              >  2:778136/1‑40 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg           >  2:195778/1‑139 (MQ=255)
                                                                                                                             tgcggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         >  1:1204999/4‑59 (MQ=255)
                                                                                                                             tgcggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                         <  2:1204999/56‑1 (MQ=255)
                                                                                                                             tgcggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                             <  2:6783/68‑1 (MQ=255)
                                                                                                                             tgcggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                             >  1:6783/4‑71 (MQ=255)
                                                                                                                                ggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                             <  1:950337/68‑1 (MQ=255)
                                                                                                                                ggggCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                             >  2:950337/1‑68 (MQ=255)
                                                                                                                                   gCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAATTAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                          >  2:1118107/1‑100 (MQ=255)
                                                                                                                                   gCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                          <  1:1118107/100‑1 (MQ=255)
                                                                                                                                   gCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga              <  1:831236/128‑1 (MQ=255)
                                                                                                                                   gCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga              >  2:831236/1‑128 (MQ=255)
                                                                                                                                     aTTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa  >  1:1472562/1‑138 (MQ=255)
                                                                                                                                     aTTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa  <  2:1472562/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA  >  NC_000913/3969888‑3970157

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