Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,470,410 (T)7→6 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.94.0% 744.6 / 38.1 185coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (7/4);  new base . (73/101);  total (80/105)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

AATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTA  >  NC_000913/2470305‑2470499
                                                                                                   |                                                                                               
aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGt                                                                                                 >  2:1742552/1‑100 (MQ=255)
    ttATGTGTTATTTTTCGGGGAACGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                             <  2:867600/100‑1 (MQ=255)
    ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                             >  1:1440915/1‑100 (MQ=255)
      aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                           >  2:2502675/1‑99 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         <  1:21920/100‑4 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         <  1:736316/100‑4 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACAAAACAGTCGttttttaat                                                                                         >  2:1702795/1‑97 (MQ=255)
         tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                        <  1:2472042/100‑5 (MQ=255)
          gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                      <  2:2671522/100‑1 (MQ=255)
            cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                    >  1:803438/1‑100 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtata                                                                                   <  1:1005581/100‑1 (MQ=255)
              tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaa                                                                                  <  1:2239966/100‑1 (MQ=255)
               ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaa                                                                                 <  1:833077/100‑1 (MQ=255)
                    ggggAAGGATATATATATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTTCTTATACATAACTGTCG‑TTTTTTAATTTATAAaataat                                                                            <  1:298966/100‑1 (MQ=255)
                     gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAtt                                                                           >  2:848230/1‑100 (MQ=255)
                     gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAtt                                                                           <  2:1551523/100‑1 (MQ=255)
                      ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  2:242184/100‑1 (MQ=255)
                      ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  2:2332059/100‑1 (MQ=255)
                      ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  2:1512647/100‑1 (MQ=255)
                       gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCa                                                                         >  1:1070382/1‑100 (MQ=255)
                        aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaa                                                                        <  2:795460/100‑1 (MQ=255)
                          ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAg                                                                      <  2:2713422/100‑1 (MQ=255)
                           gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAg                                                                      <  2:2172802/100‑1 (MQ=255)
                           gATATAGAGATGTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGa                                                                     <  1:2349669/100‑1 (MQ=255)
                             tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTTAAATAATTCAAAGaa                                                                    >  1:1092356/1‑100 (MQ=255)
                             tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGaaa                                                                   <  1:1768502/100‑1 (MQ=255)
                             tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGaaa                                                                   <  2:929268/100‑1 (MQ=255)
                               taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTTTTTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                  >  1:1329390/1‑100 (MQ=255)
                               taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                 <  2:1526980/100‑1 (MQ=255)
                               taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                 >  1:516960/1‑100 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                >  1:2772591/1‑100 (MQ=255)
                                 gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                               <  2:1212246/100‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                              <  1:949076/99‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             >  2:957691/1‑100 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             <  1:1545074/100‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             >  1:1961948/1‑100 (MQ=255)
                                    aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             <  2:1329962/100‑1 (MQ=255)
                                    aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                            <  2:50215/100‑1 (MQ=255)
                                    aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                            <  2:1075377/100‑1 (MQ=255)
                                     tttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACa                                                           <  1:383494/100‑1 (MQ=255)
                                       tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                         <  1:496229/100‑1 (MQ=255)
                                       tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                         <  2:2406882/100‑1 (MQ=255)
                                        aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                        <  2:2515975/100‑1 (MQ=255)
                                         agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                       <  1:2282382/100‑1 (MQ=255)
                                            aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                    <  1:1373307/100‑1 (MQ=255)
                                             agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                   <  2:2360304/100‑1 (MQ=255)
                                                atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                  <  1:1135167/98‑1 (MQ=255)
                                                  ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                              <  2:1614705/100‑1 (MQ=255)
                                                   tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaatg                                             >  2:1148466/1‑98 (MQ=255)
                                                   tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                             <  2:1167025/100‑1 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCTATCTTTTCTTATACATAACAGTCGTTTTTTT‑ATTTTTTTTATAATTCAAAGAAATCAAGACAATAGCGATGt                                                 >  1:1688605/1‑95 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  1:2651969/100‑1 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  2:2242047/100‑1 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  2:419987/100‑1 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  2:1197093/100‑1 (MQ=255)
                                                     tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                           <  1:660946/100‑1 (MQ=255)
                                                       tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                         >  2:612604/1‑100 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCagat                                        >  2:2275454/1‑97 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                        >  2:2572560/1‑100 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                        >  2:364023/1‑100 (MQ=255)
                                                           aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                     <  2:2597112/100‑1 (MQ=255)
                                                             gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATGATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                   >  1:1854971/1‑100 (MQ=255)
                                                             gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                   >  2:2136956/1‑100 (MQ=255)
                                                              tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGATAATAGCGATGTAAAGCACTAAGTTTa                                  >  1:1906484/1‑100 (MQ=255)
                                                                 ttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                               <  2:504836/100‑1 (MQ=255)
                                                                 ttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                               <  2:2381131/100‑1 (MQ=255)
                                                                  tttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                              >  1:587432/1‑100 (MQ=255)
                                                                  tccgatctgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGATATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                              <  1:1787648/92‑1 (MQ=255)
                                                                   ttcgttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                             <  2:1875209/100‑1 (MQ=255)
                                                                   ttcgttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                             >  1:162673/1‑100 (MQ=255)
                                                                   tgcgttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                             >  2:1966046/3‑100 (MQ=255)
                                                                      gttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAACATAATTCAACGCAATCAAGACAATAGCCATGTAAAGCACTAAGTATAAAttacac                          >  1:2773089/1‑96 (MQ=255)
                                                                       ttcgACCTTTGCTTATTCAAAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                         <  1:1702769/100‑1 (MQ=255)
                                                                       ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                          <  2:1092466/100‑1 (MQ=255)
                                                                       ttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                         <  2:2064908/100‑1 (MQ=255)
                                                                        tcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                        >  1:346779/1‑100 (MQ=255)
                                                                        tcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                        >  2:648419/1‑100 (MQ=255)
                                                                         cgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                       <  2:1052231/100‑1 (MQ=255)
                                                                         cgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                       >  1:2488604/1‑100 (MQ=255)
                                                                          gACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAttagatcgga                      >  2:1787674/1‑92 (MQ=255)
                                                                              tttGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgctattat                  <  1:2082544/100‑7 (MQ=255)
                                                                               ttGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                 <  2:2345527/100‑1 (MQ=255)
                                                                                tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                <  2:2773125/100‑1 (MQ=255)
                                                                                tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                >  2:1384733/1‑100 (MQ=255)
                                                                                  cTTATACATAACAGTCGTTTTTTT‑TTTTTTTATTTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGTACTAAGTTTAAATTGCGCGCCAATCATTgc              >  1:791746/1‑100 (MQ=255)
                                                                                  cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc              >  2:726874/1‑100 (MQ=255)
                                                                                  cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc              >  1:870610/1‑100 (MQ=255)
                                                                                   ttATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCGCTAAGTTTAAATTGCGCGCCAATCATGgcg             >  2:1663727/1‑100 (MQ=255)
                                                                                    tataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc            >  2:466180/1‑100 (MQ=255)
                                                                                    tataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc            <  2:1463638/100‑1 (MQ=255)
                                                                                     ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg           <  2:1400390/100‑1 (MQ=255)
                                                                                     ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg           <  2:461088/100‑1 (MQ=255)
                                                                                      taCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc          >  1:1269399/1‑100 (MQ=255)
                                                                                        cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac        >  1:2687030/1‑100 (MQ=255)
                                                                                        cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac        >  2:1731104/1‑100 (MQ=255)
                                                                                          tAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa      >  2:327045/1‑100 (MQ=255)
                                                                                           aaCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg     <  2:677240/100‑1 (MQ=255)
                                                                                            aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc    <  1:1203495/100‑1 (MQ=255)
                                                                                            aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc    >  1:1159846/1‑100 (MQ=255)
                                                                                              aGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta  <  1:607763/100‑1 (MQ=255)
                                                                                                   |                                                                                               
AATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTA  >  NC_000913/2470305‑2470499

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