Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,471,105 G→A pseudogene (291/291 nt)
P138L (CCT→CTT) 
tfaS →
yfdK ←
pseudogene, CPS‑53 (KpLE1) prophage, tail fiber assembly protein fragment,Phage or Prophage Related
CPS‑53 (KpLE1) prophage, conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,471,1050GA100.0% 226.2 / NA 70pseudogene (291/291 nt)
P138L (CCT→CTT) 
tfaS
yfdK
pseudogene, CPS‑53 (KpLE1) prophage, tail fiber assembly protein fragment,Phage or Prophage Related
CPS‑53 (KpLE1) prophage, conserved protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/37);  total (33/37)

TACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACT  >  NC_000913/2470972‑2471238
                                                                                                                                     |                                                                                                                                     
tACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCAt                                                                                                                                  <  2:1105828/139‑1 (MQ=255)
         ggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCt                                                                                                                         >  2:810830/1‑139 (MQ=255)
         ggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCt                                                                                                                         >  1:272494/1‑139 (MQ=255)
                    aTGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGa                                                                                                              >  2:168374/1‑139 (MQ=255)
                     tGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                    >  2:561319/1‑132 (MQ=255)
                     tGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                    <  1:561319/132‑1 (MQ=255)
                             gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     <  2:234883/139‑1 (MQ=255)
                                aCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAg                                                                                                  <  2:877398/139‑1 (MQ=255)
                                    cccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTcc                                                                                              >  1:1047802/1‑139 (MQ=255)
                                     ccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCa                                                                                             <  1:810830/139‑1 (MQ=255)
                                        tcTTAACACTGTGCTGGATTACCTTGACGCGATTACGGCAACAGATACCAACACCGCGACGGATGTCATCTGTACTGAACTGCCGGAGGACTAAGCGATTCAATATCTGGCGCACTGGAACTATCGACCAACCCCAgcg                                                                                          >  1:160909/2‑137 (MQ=14)
                                                   tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTgg                                                                                                                       >  2:517219/1‑99 (MQ=18)
                                                   tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTgg                                                                                                                       <  1:517219/99‑1 (MQ=18)
                                                   tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                               >  1:1043516/1‑139 (MQ=255)
                                                   tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                               >  1:872109/1‑139 (MQ=255)
                                                   tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGACCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                               >  2:717539/1‑139 (MQ=255)
                                                    gCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     <  2:695517/116‑1 (MQ=255)
                                                    gCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     >  1:695517/1‑116 (MQ=255)
                                                     cTGCATTACATTGACGCGGTGACGGTAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATc                                                                             <  1:717539/139‑1 (MQ=255)
                                                          ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  2:631848/139‑1 (MQ=255)
                                                          ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  2:1047802/139‑1 (MQ=255)
                                                          ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  1:686292/139‑1 (MQ=255)
                                                             cATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCcaca                                                                     >  2:816524/1‑139 (MQ=255)
                                                                  aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATta                                                                >  1:310711/1‑139 (MQ=255)
                                                                  aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATta                                                                >  1:463377/1‑139 (MQ=255)
                                                                       gTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTg                                                           >  1:359488/1‑139 (MQ=255)
                                                                       gTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTg                                                           >  2:706335/1‑139 (MQ=255)
                                                                          aCGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTgcgc                                                        <  2:65118/139‑1 (MQ=255)
                                                                          aCGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCAAGATAATCCAGCCACAAATTATATTgcgc                                                        <  2:1078100/139‑1 (MQ=255)
                                                                          aCGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAG‑‑TCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCa                                                      >  1:871463/1‑139 (MQ=255)
                                                                              cAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGt                                                    >  1:783899/1‑139 (MQ=255)
                                                                                aCAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     <  2:781393/88‑1 (MQ=255)
                                                                                aCAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     >  1:781393/1‑88 (MQ=255)
                                                                                  aGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCt                                                <  1:816524/139‑1 (MQ=255)
                                                                                                cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        >  1:593538/1‑101 (MQ=255)
                                                                                                cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  2:593538/101‑1 (MQ=255)
                                                                                                cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                  >  2:457105/1‑139 (MQ=255)
                                                                                                cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                  >  2:720050/1‑139 (MQ=255)
                                                                                                  ccGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                            >  1:204306/1‑79 (MQ=255)
                                                                                                  ccGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                            <  2:204306/79‑1 (MQ=255)
                                                                                                   cGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               <  2:474036/59‑1 (MQ=255)
                                                                                                   cGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               >  1:474036/1‑59 (MQ=255)
                                                                                                              tGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAG‑‑TCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                  <  2:871463/139‑1 (MQ=255)
                                                                                                                gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTTCCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATCATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                  <  1:706335/139‑1 (MQ=255)
                                                                                                                gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                  >  1:634308/1‑139 (MQ=255)
                                                                                                                gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                  >  2:615656/1‑139 (MQ=255)
                                                                                                                gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                  <  1:457105/139‑1 (MQ=255)
                                                                                                                gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                  <  1:168374/139‑1 (MQ=255)
                                                                                                                 ccTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               <  2:123163/45‑1 (MQ=255)
                                                                                                                 ccTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               >  1:123163/1‑45 (MQ=255)
                                                                                                                    gAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc              <  2:1036202/139‑1 (MQ=255)
                                                                                                                    gAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc              <  2:463377/139‑1 (MQ=255)
                                                                                                                     aaCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        >  1:322810/1‑80 (MQ=255)
                                                                                                                     aaCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  2:322810/80‑1 (MQ=255)
                                                                                                                     aaCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCttt                                         <  2:258989/111‑1 (MQ=255)
                                                                                                                     aaCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCttt                                         >  1:258989/1‑111 (MQ=255)
                                                                                                                      aCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGccc            <  2:634308/139‑1 (MQ=255)
                                                                                                                      aCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGccc            <  2:1043516/139‑1 (MQ=255)
                                                                                                                       cTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCa           <  2:783899/139‑1 (MQ=255)
                                                                                                                          ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     <  2:565503/46‑1 (MQ=255)
                                                                                                                          ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     >  1:565503/1‑46 (MQ=255)
                                                                                                                              aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa    <  1:945563/139‑1 (MQ=255)
                                                                                                                              aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa    <  1:720050/139‑1 (MQ=255)
                                                                                                                               ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGt                                                                                       >  1:134318/1‑55 (MQ=255)
                                                                                                                               ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGt                                                                                       <  2:134318/55‑1 (MQ=255)
                                                                                                                               ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGcc                                                       >  2:192835/1‑87 (MQ=255)
                                                                                                                               ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGccc            >  2:35010/1‑130 (MQ=255)
                                                                                                                               ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGccc            <  1:35010/130‑1 (MQ=255)
                                                                                                                               ggagTAAGCCATTCAATATATGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGcc                                                       <  1:192835/87‑1 (MQ=255)
                                                                                                                                gagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACt  <  2:310711/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACT  >  NC_000913/2470972‑2471238

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