Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,410 (T)7→6 100% coding (1280/1332 nt) gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4040T.95.5% 562.6 / 12.0 132coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (54/72);  total (56/76)
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

Only 100 of 156 total aligned reads displayed.

TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTA  >  NC_000913/2470307‑2470499
                                                                                                 |                                                                                               
 ttttTGTGTCATTTTTCGGGGAAGGATATAGCGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttn                                                                                               <  1:2959343/99‑2 (MQ=255)
  ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                             <  2:2783846/100‑1 (MQ=255)
   tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                            <  2:2733691/100‑1 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                           <  2:1294339/100‑2 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                           <  1:2002405/100‑2 (MQ=255)
     tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                          <  2:3075866/100‑3 (MQ=255)
     tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                          <  2:755969/100‑3 (MQ=255)
     tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                          <  2:705532/100‑3 (MQ=255)
     tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                          <  2:3320416/100‑3 (MQ=255)
      gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGGATGCATTTATTTGTAATGTATTTTTCGTTCGTCCTTTGCTTATACATAACAGTCGttttttaat                                                                                         <  2:827937/100‑4 (MQ=255)
      gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         <  1:1164626/100‑4 (MQ=255)
      gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         <  2:2195397/100‑4 (MQ=255)
      gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         >  2:2970339/1‑97 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                        <  2:2875051/100‑5 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                        <  2:3353547/100‑5 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTTTTTTTACTTATACATAACAGTCGTTTTTTTAttt                                                                                       >  1:39598/1‑100 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttttttt                                                                                       >  1:3028447/1‑96 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                      <  2:2222722/100‑1 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                      >  2:1034691/1‑100 (MQ=255)
         tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTta                                                                                     <  1:2535146/100‑1 (MQ=255)
         tCATTTTTCAGGTAAGGATATAGGGATGTAAGCAAGCAGTTATTTGCAATGTATTTTTTGTTCGAACTTTGCTTATACATACCAGTCGttttttgatt                                                                                        <  1:210163/98‑5 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                    >  1:1461906/1‑100 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                    <  2:3432355/100‑1 (MQ=255)
             ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                 <  2:2368085/100‑1 (MQ=255)
             ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                 <  1:2884246/100‑1 (MQ=255)
              tttCGGGGGAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaaa                                                                                >  2:369614/1‑100 (MQ=255)
               ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTACAAaat                                                                               <  1:1811010/100‑1 (MQ=255)
                 cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                             >  1:149517/1‑100 (MQ=255)
                   gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                           >  2:1868343/1‑100 (MQ=255)
                   gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                           <  2:2635099/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                          <  1:734139/100‑1 (MQ=255)
                       aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCTTTTTATTTTTACTTATACATAACATTCGTTTTTTTATTTATAAAATAATTCaa                                                                        >  1:68620/1‑99 (MQ=255)
                       aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                       <  2:2803501/100‑1 (MQ=255)
                       aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                       <  1:3605471/100‑1 (MQ=255)
                         gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                     >  1:3240932/1‑100 (MQ=255)
                           tataGATATTTAAGCAATTATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTAATTTATAAAATAATTCAAAg                                                                      <  1:3540810/97‑1 (MQ=255)
                             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                  <  1:1419384/99‑1 (MQ=255)
                             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                 >  1:3046957/1‑100 (MQ=255)
                             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                 >  2:1851416/1‑100 (MQ=255)
                              agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                >  1:1424453/1‑100 (MQ=255)
                              agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                <  2:3258302/100‑1 (MQ=255)
                               gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATACTTCAAAGAAATc                                                                 >  2:3171182/1‑98 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             >  1:521019/1‑100 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             >  2:3462108/1‑100 (MQ=255)
                                  aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                            >  2:1791198/1‑100 (MQ=255)
                                   tttAAGCAAGCATTTATTTGGAATGTATTTTTCTTTTTATTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAATATAATTCAAAGAAATCAAGAc                                                            >  1:3393742/1‑99 (MQ=255)
                                    ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATTAAATAATTCAAAGAAATCAAGACaa                                                          >  1:3569517/1‑100 (MQ=255)
                                     tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                         >  2:2857580/1‑100 (MQ=255)
                                      aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                        <  2:360395/100‑1 (MQ=255)
                                       agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                       >  2:322365/1‑100 (MQ=255)
                                         caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTCTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                     <  2:2665562/100‑1 (MQ=255)
                                          aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                    <  2:444131/100‑1 (MQ=255)
                                            gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                  <  2:3181902/100‑1 (MQ=255)
                                            gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                  <  2:661490/100‑1 (MQ=255)
                                             caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATTGCGATGt                                                 <  2:1995285/100‑1 (MQ=255)
                                              atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATCACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                <  1:2420264/100‑1 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGATTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGAGATGTa                                                <  1:1006771/99‑1 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  1:1481009/100‑1 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  2:2334154/100‑1 (MQ=255)
                                                     tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                         <  2:552456/100‑1 (MQ=255)
                                                      ggAATGTATTTTTCGTTTGACCTTTGCTTATACATAACAGTAGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                        <  1:909041/100‑1 (MQ=255)
                                                       gAATGTATTTTTCGTTCGACCTTTACTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                        <  1:948891/99‑1 (MQ=255)
                                                        aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                      <  1:3339936/100‑1 (MQ=255)
                                                          tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                    <  1:2008518/100‑1 (MQ=255)
                                                           gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                   >  2:3301811/1‑100 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                  <  2:1792151/100‑1 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                  <  1:1566519/100‑1 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                  <  2:2666283/100‑1 (MQ=255)
                                                             ntttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                 >  1:1017508/2‑100 (MQ=255)
                                                             aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                 >  2:2186733/1‑100 (MQ=255)
                                                             aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                 <  2:887289/100‑1 (MQ=255)
                                                             aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                 <  1:485468/100‑1 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                               <  2:2230439/100‑1 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                               >  1:3227498/1‑100 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                               >  2:3689845/1‑100 (MQ=255)
                                                                 ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                             <  1:1535729/100‑1 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                           >  2:1471674/1‑100 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                           >  1:3086194/1‑100 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATCTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                           >  1:266076/1‑100 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                          <  1:1770893/100‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                         <  2:983285/100‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                         <  2:345317/100‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                         >  2:1284591/1‑100 (MQ=255)
                                                                      tcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTtcgcgc                        >  1:606554/1‑100 (MQ=255)
                                                                         aCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                       <  1:984910/98‑1 (MQ=255)
                                                                         aCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                     >  1:3466716/1‑100 (MQ=255)
                                                                          ccTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCTCTCCAAt                    >  1:1239032/1‑100 (MQ=255)
                                                                           cTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                   <  2:472070/100‑1 (MQ=255)
                                                                               gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg               <  2:1131064/100‑1 (MQ=255)
                                                                               gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg               <  1:3305212/100‑1 (MQ=255)
                                                                                cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc              >  2:3254592/1‑100 (MQ=255)
                                                                                  tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc            >  2:533564/1‑100 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc          >  2:2520170/1‑100 (MQ=255)
                                                                                     aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca         >  1:2915920/1‑100 (MQ=255)
                                                                                        tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa      >  2:3698296/1‑100 (MQ=255)
                                                                                        tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa      >  2:1675281/1‑100 (MQ=255)
                                                                                         aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg     >  2:1390516/1‑100 (MQ=255)
                                                                                         aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg     <  2:19606/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt   <  2:3381746/100‑1 (MQ=255)
                                                                                            aGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta  <  2:2124117/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTA  >  NC_000913/2470307‑2470499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: