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.7% 393.9 / 6.3 93coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (39/50);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 104 total aligned reads displayed.

TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470307‑2470498
                                                                                                 |                                                                                              
ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttt                                                                                              <  1:75166/100‑1 (MQ=255)
 tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAATCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtt                                                                                               <  2:2155697/98‑1 (MQ=255)
 tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                             <  1:1985724/100‑1 (MQ=255)
 tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                             <  2:1792967/100‑1 (MQ=255)
 ccgatctgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                             <  1:1056520/94‑1 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          <  1:331965/100‑2 (MQ=255)
     tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                         >  2:1233484/1‑98 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  2:1249508/100‑5 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  2:1494271/100‑5 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttagatcgg                                                                                      >  2:1056516/1‑93 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     >  1:2382487/1‑100 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTCCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     <  1:1638471/100‑1 (MQ=255)
         ttaTTTTTCGAGTAAGGATATAAATCTTTAAGCAATCATTTATTTGAAATGTATTTTTCGTACGACCTTTGCTTAGACATAACAGTCGTTTTTTAATTta                                                                                    <  2:1408987/98‑1 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATCACAGTCGTTTTTTAATTtat                                                                                   >  1:1991667/1‑100 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                   >  1:136162/1‑100 (MQ=255)
           aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATCTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  <  2:107956/100‑1 (MQ=255)
            tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                 <  1:145656/100‑1 (MQ=255)
               ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaat                                                                              <  1:932350/100‑1 (MQ=255)
                  ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           >  2:370198/1‑100 (MQ=255)
                  ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           >  2:1052515/1‑100 (MQ=255)
                   gggAAGGATATAGAGCTTTAAGCAAGCGTTTATCTGGAATGTATTTTAAGATTCACTTTTGCAAATATATAACAGGCGGAGTTTAAGTTAGAAAATAAtt                                                                          >  2:438119/1‑100 (MQ=255)
                   gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                          <  1:381801/100‑1 (MQ=255)
                   gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                          >  1:1646627/1‑100 (MQ=255)
                    ggAAGGATATATAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         <  2:1117598/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         >  2:471297/1‑100 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         <  1:2108281/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         <  1:1315054/100‑1 (MQ=255)
                     gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         <  2:1694061/99‑1 (MQ=255)
                        ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     <  1:1595631/100‑1 (MQ=255)
                         gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCTTTCTACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                    >  2:78994/1‑100 (MQ=255)
                         gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                    >  1:1072853/1‑100 (MQ=255)
                          atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaa                                                                   <  2:359925/100‑1 (MQ=255)
                          atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaa                                                                   >  1:598302/1‑100 (MQ=255)
                           tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGATCTTTGCGTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                    <  2:1426900/98‑1 (MQ=255)
                           tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  <  2:1613993/100‑1 (MQ=255)
                            ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                 <  1:1158093/100‑1 (MQ=255)
                             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                >  1:1379652/1‑100 (MQ=255)
                              agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  1:1721363/1‑100 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  1:1311725/1‑100 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTAGCTTATACATAACAGTCCTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                              <  2:1113564/98‑1 (MQ=255)
                                  aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  2:2100083/100‑1 (MQ=255)
                                       agtaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGATTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                      <  2:2388078/97‑1 (MQ=255)
                                       agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                      >  2:237376/1‑100 (MQ=255)
                                        gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCGAGACAATAGc                                                     <  1:78995/100‑1 (MQ=255)
                                        gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     <  2:11460/100‑1 (MQ=255)
                                        gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     <  2:1098758/100‑1 (MQ=255)
                                         ctagctTTTATTTGGAATGAATTTTTTGTTCAACCTTTGCTTATACATGACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  2:365573/91‑1 (MQ=255)
                                         caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  1:281724/1‑100 (MQ=255)
                                          aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                   <  1:153296/100‑1 (MQ=255)
                                          aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                   <  1:2332014/100‑1 (MQ=255)
                                          aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                   <  1:664544/100‑1 (MQ=255)
                                           agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                  >  2:1832158/1‑100 (MQ=255)
                                            gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  2:467038/100‑1 (MQ=255)
                                             caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  2:92969/99‑1 (MQ=255)
                                             caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                >  1:1076221/1‑100 (MQ=255)
                                              atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  2:2171084/98‑1 (MQ=255)
                                              atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               <  2:1072849/100‑1 (MQ=255)
                                              atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               <  2:121567/100‑1 (MQ=255)
                                              atttatttGGAATGTATTTTTCGTTCGACCATTGATTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  2:314009/98‑1 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                              >  1:2027074/1‑100 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                              >  1:1183547/1‑100 (MQ=255)
                                                  atttGGAATGTATTTTTTTTTTTTTTTTTTCTTATACATAACAGTCGTTTTTTAATTTATAAAGTAATTCAAAGAAATTAAGACAATAGCGATGTAAAg                                            >  2:1314796/1‑99 (MQ=255)
                                                   tttGGAATGTATTTTTTTTTCGTTCTTTCCATATACATAGCAGTCATTTTTTATTTTATTTACTAATTCAAAGAAATTAAAACAATAACGTTGTAAAga                                           >  2:1094513/1‑98 (MQ=255)
                                                     tGGAATGTATTTTTTTTTTTTTTTTTACTTATACATAACAGTCGTTTTTTAATTTTTAAAATAATTCAAAGAAATTAATACAATAGCGATGTAAAGCACt                                        >  2:410281/1‑100 (MQ=255)
                                                     tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                        >  2:2197177/1‑100 (MQ=255)
                                                     tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                        >  2:1965584/1‑100 (MQ=255)
                                                        aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     <  2:644985/100‑1 (MQ=255)
                                                        aaTGTATTTTTCGTTAGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     <  2:1670919/100‑1 (MQ=255)
                                                           gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                  >  1:99280/1‑100 (MQ=255)
                                                             aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTTTTTTTTAATAATTCAAAGAAATCAAGACAATAGCTATGTAAATCACTAAGTTTaa                                >  2:725144/1‑100 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:165645/100‑1 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:2357255/100‑1 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:2155289/100‑1 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         <  1:1368992/100‑1 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         <  1:2097370/100‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                        <  1:2279602/100‑1 (MQ=255)
                                                                        gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATCAAATAATTCAAAGATATCAAGTCAATAAGGATGCTATGAACTAAGTTTAAAt                              >  2:905056/1‑91 (MQ=255)
                                                                        gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATTGCGATGTAAACCACTAAGTTTAAATTGCGCGCCa                     >  2:1418678/1‑100 (MQ=255)
                                                                        gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                      <  2:2313740/99‑1 (MQ=255)
                                                                        gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                     >  1:1469164/1‑100 (MQ=255)
                                                                            tttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                 >  1:303521/1‑100 (MQ=255)
                                                                              tGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg               <  1:2084333/100‑1 (MQ=255)
                                                                               gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              <  1:2284488/100‑1 (MQ=255)
                                                                                  tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc           >  2:1534014/1‑100 (MQ=255)
                                                                                   ataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg          <  1:1458254/100‑1 (MQ=255)
                                                                                   ataCATAACAGTCGTCTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg          <  1:2351711/100‑1 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         >  2:249243/1‑100 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         >  1:2005467/1‑100 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         >  1:192095/1‑100 (MQ=255)
                                                                                     aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca        <  2:267053/100‑1 (MQ=255)
                                                                                     aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca        <  1:1302581/100‑1 (MQ=255)
                                                                                     aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca        >  2:522802/1‑100 (MQ=255)
                                                                                      cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac       >  2:978570/1‑100 (MQ=255)
                                                                                      cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac       >  1:1880902/1‑100 (MQ=255)
                                                                                        tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca      >  1:2036106/1‑100 (MQ=255)
                                                                                        tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  1:864024/1‑100 (MQ=255)
                                                                                        tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  1:2216680/1‑100 (MQ=255)
                                                                                         aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg    <  1:677494/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTAATTTGTAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  <  1:445844/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  <  2:966478/100‑1 (MQ=255)
                                                                                                 |                                                                                              
TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470307‑2470498

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