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% 247.4 / NA 77pseudogene (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 (35/42);  total (35/42)

ATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTA  >  NC_000913/2470970‑2471236
                                                                                                                                       |                                                                                                                                   
aTTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGcc                                                                                                                                  <  2:1108570/139‑1 (MQ=255)
  tACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCAt                                                                                                                                <  2:1304901/139‑1 (MQ=255)
         tGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCAtt                                                                                                                               >  2:677553/1‑133 (MQ=255)
         tGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCAtt                                                                                                                               <  1:677553/133‑1 (MQ=255)
          gggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATc                                                                                                                        >  2:977828/1‑139 (MQ=255)
           ggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCt                                                                                                                       >  1:93803/1‑139 (MQ=255)
           ggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCt                                                                                                                       >  1:1495561/1‑139 (MQ=255)
            gCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTg                                                                                                                      <  2:1495561/139‑1 (MQ=255)
               gAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcg                                                                                                                   <  1:128475/139‑1 (MQ=255)
               gAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcg                                                                                                                   <  2:1184048/139‑1 (MQ=255)
                aaGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                  <  1:692661/139‑1 (MQ=255)
                 aGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                  >  2:427078/1‑138 (MQ=255)
                 aGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                  <  1:427078/138‑1 (MQ=255)
                 aGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCa                                                                                                                 >  1:1281428/1‑139 (MQ=255)
                     gATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                             >  1:214840/1‑139 (MQ=255)
                            gcagAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATc                                                                                                      >  2:1456231/1‑139 (MQ=255)
                             cagAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCg                                                                                                     >  2:636314/1‑139 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  2:30513/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  1:1145308/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  2:1281428/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  1:654865/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  1:636314/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  2:766751/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  1:1565092/139‑1 (MQ=255)
                               gagaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  2:93803/139‑1 (MQ=255)
                                agaCCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAcc                                                                                                  >  2:506195/1‑139 (MQ=255)
                                        cGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                          <  1:1225153/139‑1 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCa                                                                                                                             <  1:1599/91‑1 (MQ=21)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCa                                                                                                                             >  2:1599/1‑91 (MQ=18)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                  <  1:1338908/102‑1 (MQ=18)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcgc                                                                                                                  >  2:1338908/1‑102 (MQ=17)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  1:1263402/117‑1 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   >  2:1263402/1‑117 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                             >  1:1373463/1‑139 (MQ=255)
                                                     tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                             >  1:1036696/1‑139 (MQ=255)
                                                      gCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   >  2:1176363/1‑116 (MQ=255)
                                                      gCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                   <  1:1176363/116‑1 (MQ=255)
                                                            ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      <  2:208587/139‑1 (MQ=255)
                                                             tACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCGGTGCGTCCAGATAATCCAGCca                                                                     >  2:104813/1‑139 (MQ=255)
                                                               cATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCcaca                                                                   >  1:312465/1‑139 (MQ=255)
                                                                    aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      >  1:170128/1‑131 (MQ=255)
                                                                    aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      <  2:170128/131‑1 (MQ=255)
                                                                    aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATta                                                              >  2:83578/1‑139 (MQ=255)
                                                                    aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATta                                                              >  1:51213/1‑139 (MQ=255)
                                                                                   cAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      >  1:450090/1‑116 (MQ=255)
                                                                                   cAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      <  2:450090/116‑1 (MQ=255)
                                                                                           aGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCttt                                       <  1:506195/139‑1 (MQ=255)
                                                                                             cACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCa                                     >  2:219519/1‑139 (MQ=255)
                                                                                                 gcgcCGGATGTCATCTGGCCTGAACTACCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacg                                 >  2:739343/1‑139 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                             >  2:614306/1‑62 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                             <  1:614306/62‑1 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                          >  1:277247/1‑81 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                          <  2:277247/81‑1 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                >  1:482660/1‑139 (MQ=255)
                                                                                                  cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                >  2:184407/1‑139 (MQ=255)
                                                                                                        aTGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCGCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATa                          >  2:712850/1‑139 (MQ=255)
                                                                                                            cATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCg                      <  1:184407/139‑1 (MQ=255)
                                                                                                            cATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCg                      <  1:379469/139‑1 (MQ=255)
                                                                                                                  gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCGGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                <  1:104813/139‑1 (MQ=255)
                                                                                                                  gCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                <  1:1456231/139‑1 (MQ=255)
                                                                                                                    cTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGc                     <  2:1135649/132‑1 (MQ=255)
                                                                                                                    cTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGc                     >  1:1135649/1‑132 (MQ=255)
                                                                                                                      gAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc            <  1:83578/139‑1 (MQ=255)
                                                                                                                      gAACTGCCGGAGGAGTAAGCCATTCAATATCAGGCGCACTGGAAGTATCGACCAGCTCTAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc            <  2:1036696/139‑1 (MQ=255)
                                                                                                                      gAACTACCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc            <  1:739343/139‑1 (MQ=255)
                                                                                                                        aCTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGccc          <  2:312465/139‑1 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      <  2:1484397/75‑1 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                      >  1:1484397/1‑75 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc            <  2:475643/133‑1 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc            >  1:475643/1‑133 (MQ=255)
                                                                                                                            ccGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTg      >  1:64113/1‑139 (MQ=255)
                                                                                                                              ggaggaGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGtt    <  2:51213/139‑1 (MQ=255)
                                                                                                                                aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa  <  2:482660/139‑1 (MQ=255)
                                                                                                                                 ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATa                          <  2:386703/114‑1 (MQ=255)
                                                                                                                                 ggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATa                          >  1:386703/1‑114 (MQ=255)
                                                                                                                                    gTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACaaa                                                                 <  1:62699/72‑1 (MQ=255)
                                                                                                                                    gTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACaaa                                                                 >  2:62699/1‑72 (MQ=255)
                                                                                                                                       |                                                                                                                                   
ATTACAGGCTGGGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTA  >  NC_000913/2470970‑2471236

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: