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.98.7% 721.5 / 1.3 157coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (72/83);  total (74/83)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 179 total aligned reads displayed.

AATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470305‑2470498
                                                                                                   |                                                                                              
   tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                             <  2:143735/100‑1 (MQ=255)
   tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                             <  2:2099893/100‑1 (MQ=255)
    ttATGTGTCATTTTTCGGGGAAGGATATCGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                            <  1:3051056/100‑1 (MQ=255)
    ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                            >  2:790332/1‑100 (MQ=255)
      aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          >  2:2172392/1‑99 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        >  1:1528853/1‑97 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        >  2:2924831/1‑97 (MQ=255)
         tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       >  2:840127/1‑96 (MQ=255)
         tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       >  2:2856428/1‑96 (MQ=255)
         tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       >  1:2741960/1‑96 (MQ=255)
         tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  1:343462/100‑5 (MQ=255)
         tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  1:2311866/100‑5 (MQ=255)
          gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     >  2:207919/1‑100 (MQ=255)
          gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     <  2:1489969/100‑1 (MQ=255)
          gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     >  2:2623990/1‑100 (MQ=255)
            cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                   >  2:297312/1‑100 (MQ=255)
            cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                   >  1:636784/1‑100 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  >  2:1036641/1‑100 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  <  1:914394/100‑1 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  <  2:1271681/100‑1 (MQ=255)
              tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  <  1:187805/99‑1 (MQ=255)
              tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                 <  2:1753970/100‑1 (MQ=255)
              tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                 <  2:1180464/100‑1 (MQ=255)
               ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                <  1:2701549/100‑1 (MQ=255)
                   cgtggAAGGATATAGAGACTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTAGACCTTAGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                            <  1:2073228/97‑1 (MQ=255)
                    ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           >  2:1549421/1‑100 (MQ=255)
                    ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           >  1:2462906/1‑100 (MQ=255)
                    ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           <  2:2823721/100‑1 (MQ=255)
                      ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         >  2:1794402/1‑100 (MQ=255)
                         aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                      >  1:2421709/1‑100 (MQ=255)
                          ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     <  2:2281192/100‑1 (MQ=255)
                          ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     >  1:644812/1‑100 (MQ=255)
                             tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  >  2:95701/1‑100 (MQ=255)
                              ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                 <  2:385383/100‑1 (MQ=255)
                               taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                >  2:2350046/1‑100 (MQ=255)
                               taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                >  1:920531/1‑100 (MQ=255)
                               taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                <  2:16638/100‑1 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATTCATAACAGTCGTTTTTTCATTTATAAAATAATTCAAAGAAATc                                                                >  1:849165/1‑99 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  2:1003087/1‑100 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  2:1176694/1‑100 (MQ=255)
                                 gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  2:2748880/1‑100 (MQ=255)
                                 gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  1:2772481/1‑99 (MQ=255)
                                  agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                             <  2:2806464/100‑1 (MQ=255)
                                  agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                             <  2:1475832/100‑1 (MQ=255)
                                   gaTTTAGGCAAGCATTTATTTTGAATGTATTTTTCTTTCGACCTTTACTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                 <  1:1002253/95‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                              <  1:1598689/98‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               <  1:3028746/97‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  2:1694016/1‑100 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            <  2:2517839/100‑1 (MQ=255)
                                    aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  2:636789/100‑1 (MQ=255)
                                     tttAAGCAAGCATTTATTTGTAATGTATTTTTTGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  1:657919/99‑1 (MQ=255)
                                       tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                        <  1:1036635/100‑1 (MQ=255)
                                         agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                      >  1:2398189/1‑100 (MQ=255)
                                          gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     <  2:2741976/100‑1 (MQ=255)
                                          gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     >  1:1257024/1‑100 (MQ=255)
                                           caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  1:1080688/1‑100 (MQ=255)
                                            aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                   <  1:790327/100‑1 (MQ=255)
                                             agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGAt                                                  >  1:457671/1‑100 (MQ=255)
                                                atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  1:2494413/98‑1 (MQ=255)
                                                atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               >  1:2942709/1‑100 (MQ=255)
                                                atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               <  1:1505668/100‑1 (MQ=255)
                                                  ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                             <  2:3033559/100‑1 (MQ=255)
                                                   tatttGGAATGTATTTTTCGTTCGATTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAATTAATTCAAAGAAATCAAGACAATAGCGATGt                                                >  1:1391942/1‑96 (MQ=255)
                                                     tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                          >  2:1855502/1‑100 (MQ=255)
                                                      ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCan                                         <  2:1611975/100‑2 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                       >  1:1312665/1‑100 (MQ=255)
                                                           aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                    >  2:798914/1‑100 (MQ=255)
                                                            tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                   >  1:2919449/1‑100 (MQ=255)
                                                             gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                  >  2:897677/1‑100 (MQ=255)
                                                              tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 >  2:447284/1‑100 (MQ=255)
                                                              tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 >  1:497633/1‑100 (MQ=255)
                                                               aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                >  2:2143363/1‑100 (MQ=255)
                                                                tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                               <  2:2996368/100‑1 (MQ=255)
                                                                tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                               <  2:497637/100‑1 (MQ=255)
                                                                tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                               <  2:1342644/100‑1 (MQ=255)
                                                                      gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGTTGTAAAGCACTAAGTTTAAAttac                           >  1:1237695/1‑96 (MQ=255)
                                                                      gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         >  1:2705408/1‑100 (MQ=255)
                                                                      gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         <  1:1549413/100‑1 (MQ=255)
                                                                       ttcgATCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                        <  1:2914975/100‑1 (MQ=255)
                                                                        tcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                       <  2:2661495/100‑1 (MQ=255)
                                                                        tcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                       <  2:2306422/100‑1 (MQ=255)
                                                                          gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                     <  2:1404099/100‑1 (MQ=255)
                                                                           aCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                    <  2:2546918/100‑1 (MQ=255)
                                                                              tttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                 <  2:2664058/100‑1 (MQ=255)
                                                                                 gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              <  1:1442767/100‑1 (MQ=255)
                                                                                 gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              <  1:2350033/100‑1 (MQ=255)
                                                                                 gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              <  2:1720967/100‑1 (MQ=255)
                                                                                  cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:321311/1‑100 (MQ=255)
                                                                                  cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:1786371/1‑100 (MQ=255)
                                                                                   ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCTAATCATGgcg            >  2:1955438/1‑100 (MQ=255)
                                                                                   ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg            <  1:1719323/100‑1 (MQ=255)
                                                                                      taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         >  2:1638292/1‑100 (MQ=255)
                                                                                        cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac       >  1:2204813/1‑100 (MQ=255)
                                                                                        cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac       >  1:282083/1‑100 (MQ=255)
                                                                                         aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca      >  1:180952/1‑100 (MQ=255)
                                                                                         aTAACACTCGTTGTTTAATTTATAAAATAATTCAAAGATATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGca        <  1:423729/98‑1 (MQ=255)
                                                                                          tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  2:2430359/1‑100 (MQ=255)
                                                                                          tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  2:408331/1‑100 (MQ=255)
                                                                                           aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg    <  2:2108088/100‑1 (MQ=255)
                                                                                             cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  <  2:1077456/100‑1 (MQ=255)
                                                                                                   |                                                                                              
AATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470305‑2470498

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