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% 227.7 / NA 71pseudogene (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 (32/39);  total (32/39)

ATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAC  >  NC_000913/2470970‑2471237
                                                                                                                                       |                                                                                                                                    
aTTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGcc                                                                                                                                   <  2:554636/139‑1 (MQ=255)
          gggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATc                                                                                                                         >  2:145019/1‑139 (MQ=255)
           ggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCt                                                                                                                        >  1:835059/1‑139 (MQ=255)
             cGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTgg                                                                                                                      <  1:880620/139‑1 (MQ=255)
                     gATGCAGGCAGAGTCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                              <  2:926755/139‑1 (MQ=255)
                     gATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                              <  1:835683/139‑1 (MQ=255)
                     gATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                              >  1:57327/1‑139 (MQ=255)
                            gcagAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATc                                                                                                       >  1:714098/1‑139 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    <  1:570398/139‑1 (MQ=255)
                                agaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    >  2:654896/1‑138 (MQ=255)
                                agaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    <  1:654896/138‑1 (MQ=255)
                                      cccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    >  2:99396/1‑132 (MQ=255)
                                      cccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    <  1:99396/132‑1 (MQ=255)
                                      cccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGACGCACTGGAAGTATCGACCAGCTcc                                                                                             >  2:312597/1‑139 (MQ=255)
                                       ccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCa                                                                                            <  2:1024756/139‑1 (MQ=255)
                                        cGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                           <  2:57327/139‑1 (MQ=255)
                                        cGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                           <  2:714098/139‑1 (MQ=255)
                                             tAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGt                                                                                      <  1:145019/139‑1 (MQ=255)
                                              aaCACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTc                                                                                     >  2:783811/1‑139 (MQ=255)
                                               acacTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    >  1:30676/1‑123 (MQ=255)
                                               acacTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    <  2:30676/123‑1 (MQ=255)
                                                    gtgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATa                                                                               >  1:846266/1‑139 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                              >  1:701398/1‑139 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                              >  1:538320/1‑139 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                              >  1:324437/1‑139 (MQ=255)
                                                      gCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAAt                                                                             >  2:674112/1‑139 (MQ=255)
                                                         ggATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAg                                                                         <  2:835059/139‑1 (MQ=255)
                                                          gATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCAGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATACTCCAg                                                                         >  1:668368/1‑139 (MQ=255)
                                                           aTTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGc                                                                        <  2:668368/139‑1 (MQ=255)
                                                            ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                       <  2:976179/139‑1 (MQ=255)
                                                            ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGACGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                       <  1:312597/139‑1 (MQ=255)
                                                               cATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCcaca                                                                    >  2:589461/1‑139 (MQ=255)
                                                                    aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCa                                                                                                                              >  1:478787/1‑76 (MQ=21)
                                                                    aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCa                                                                                                                              <  2:478787/76‑1 (MQ=21)
                                                                     cgcgGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                              <  2:878389/91‑1 (MQ=18)
                                                                     cgcgGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                              >  1:878389/1‑91 (MQ=18)
                                                                       cgGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTatat                                                            >  1:427948/1‑139 (MQ=255)
                                                                       cgGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTatat                                                            <  2:846266/139‑1 (MQ=255)
                                                                          tGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    <  1:256397/96‑1 (MQ=255)
                                                                          tGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    >  2:256397/1‑96 (MQ=255)
                                                                              ggCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCa                                                     >  1:746653/1‑139 (MQ=255)
                                                                                   cAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTcc                                                >  2:796313/1‑139 (MQ=255)
                                                                                   cAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATAGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTcc                                                >  2:941233/1‑139 (MQ=255)
                                                                                    aGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCt                                               <  2:746653/139‑1 (MQ=255)
                                                                                           aGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    <  2:852189/79‑1 (MQ=255)
                                                                                           aGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                    >  1:852189/1‑79 (MQ=255)
                                                                                             cACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                           <  2:428176/86‑1 (MQ=255)
                                                                                             cACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                           >  1:428176/1‑86 (MQ=255)
                                                                                              aCCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAg                                     >  1:249085/1‑139 (MQ=255)
                                                                                               ccGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAga                                    <  2:324437/139‑1 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                   >  1:311145/1‑57 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                   <  2:311145/57‑1 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                           >  1:680917/1‑81 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                       >  1:691327/1‑101 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                       <  2:691327/101‑1 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCAAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                 >  1:642789/1‑139 (MQ=255)
                                                                                                  cgcAGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                           <  2:680917/81‑1 (MQ=255)
                                                                                                        aTGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATCGCGCCAGTTCCTCACCTTTCAGACGACCAATa                           <  2:639736/139‑1 (MQ=255)
                                                                                                                 ggCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCATAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTa                  <  1:783811/139‑1 (MQ=255)
                                                                                                                      gAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc             <  1:674112/139‑1 (MQ=255)
                                                                                                                      gAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc             <  2:642789/139‑1 (MQ=255)
                                                                                                                        aCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGccc           <  2:249085/139‑1 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc             <  1:472231/133‑1 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc             >  2:472231/1‑133 (MQ=255)
                                                                                                                               gaggagTAAGCCATTCAATATCTGGCGCATTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCAtttttt    <  1:941233/139‑1 (MQ=255)
                                                                                                                               gaggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGttt    <  1:796313/139‑1 (MQ=255)
                                                                                                                                aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa   <  2:427948/139‑1 (MQ=255)
                                                                                                                                aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa   <  2:701398/139‑1 (MQ=255)
                                                                                                                                 ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTgcgc                                                       >  2:653438/1‑86 (MQ=255)
                                                                                                                                 ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTgcgc                                                       <  1:653438/86‑1 (MQ=255)
                                                                                                                                 ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAc  <  1:589461/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
ATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAC  >  NC_000913/2470970‑2471237

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