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.96.7% 772.0 / 18.5 182coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base . (82/93);  total (89/93)
Fisher's exact test for biased strand distribution p-value = 1.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 202 total aligned reads displayed.

AATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470305‑2470498
                                                                                                   |                                                                                              
 aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtt                                                                                               >  1:2518055/1‑100 (MQ=255)
   tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                             <  2:454035/100‑1 (MQ=255)
      aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          <  1:804545/100‑2 (MQ=255)
      aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          <  1:2757558/100‑2 (MQ=255)
      aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          <  2:442502/100‑2 (MQ=255)
       tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                         <  2:1972638/100‑3 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        <  2:1177123/100‑4 (MQ=255)
        gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        <  1:2036873/100‑4 (MQ=255)
          gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                     >  2:1746115/1‑100 (MQ=255)
          gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                     >  2:2867435/1‑100 (MQ=255)
            cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                   >  2:2667287/1‑100 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtata                                                                                  <  1:1662172/100‑1 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtata                                                                                  <  1:187369/100‑1 (MQ=255)
             aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtata                                                                                  >  2:2685059/1‑100 (MQ=255)
              tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaa                                                                                 <  1:1576261/100‑1 (MQ=255)
               ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaa                                                                                >  2:1254503/1‑100 (MQ=255)
                tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                               >  2:2205848/1‑100 (MQ=255)
                tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                               >  1:2907348/1‑100 (MQ=255)
                  tCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaata                                                                             <  1:1438172/100‑1 (MQ=255)
                    ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaataat                                                                           <  2:797272/100‑1 (MQ=255)
                      ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                         <  1:1392148/100‑1 (MQ=255)
                      ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                         >  2:2325747/1‑100 (MQ=255)
                       gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCa                                                                        <  2:825170/100‑1 (MQ=255)
                        aaGGATATAGTTATTTAAGCAATCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGATTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaa                                                                       <  2:1999915/100‑1 (MQ=255)
                         aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaaa                                                                      <  1:540830/100‑1 (MQ=255)
                          ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTGTTCGACCTTTACTTATACATAACAGTCG‑TTTTTTAATTTATAAaataat                                                                           <  2:521532/94‑1 (MQ=255)
                          ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAg                                                                     <  2:1887650/100‑1 (MQ=255)
                           gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGa                                                                    <  1:254887/100‑1 (MQ=255)
                             tataGAGTTTTAAGCAAGTTTTTATTTGGAATGTATTTTTCTTTTGATCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGaa                                                                   <  2:2289344/99‑1 (MQ=255)
                             tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  >  2:2638662/1‑100 (MQ=255)
                              atgtatgttttAGCAATCTTTTATTTGGTATTTATTTTTTGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaaa                                                                      <  2:1653974/88‑1 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               <  1:2155125/100‑1 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  1:1655134/1‑100 (MQ=255)
                                agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  2:613814/1‑100 (MQ=255)
                                 gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCTTTCGACCTTTACTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                <  2:1761003/98‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGATCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                             <  2:1808174/99‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                              <  2:1104231/98‑1 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  2:675550/1‑100 (MQ=255)
                                   gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  1:886463/1‑100 (MQ=255)
                                     tttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACa                                                          <  1:2308126/100‑1 (MQ=255)
                                      ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACaa                                                         >  2:641814/1‑100 (MQ=255)
                                       tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                        >  2:947224/1‑100 (MQ=255)
                                       tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                        <  1:369348/100‑1 (MQ=255)
                                        aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                       >  1:425901/1‑100 (MQ=255)
                                           caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  2:888950/1‑100 (MQ=255)
                                           caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  1:994589/1‑100 (MQ=255)
                                           caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  1:603023/1‑100 (MQ=255)
                                             agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                  <  2:1011328/100‑1 (MQ=255)
                                              gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  1:2561912/100‑1 (MQ=255)
                                              gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  1:289722/100‑1 (MQ=255)
                                                 tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                              <  1:2645385/100‑1 (MQ=255)
                                                   tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                            >  2:568546/1‑100 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATGGCGATGTAAAGc                                           <  1:2008073/100‑1 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           >  2:1331569/1‑100 (MQ=255)
                                                    atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           <  1:1541577/100‑1 (MQ=255)
                                                      ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                         <  1:1077829/100‑1 (MQ=255)
                                                       tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCACt                                        >  2:1283775/1‑100 (MQ=255)
                                                       tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                        >  2:1359833/1‑100 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                       >  2:2797466/1‑100 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                       >  1:2386386/1‑100 (MQ=255)
                                                        ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                       <  1:1125495/100‑1 (MQ=255)
                                                          aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     <  1:1321232/100‑1 (MQ=255)
                                                          aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     <  2:616725/100‑1 (MQ=255)
                                                             gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                  <  1:3010345/100‑1 (MQ=255)
                                                              tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 >  1:1182276/1‑100 (MQ=255)
                                                               aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                <  1:768122/100‑1 (MQ=255)
                                                               aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                <  1:664414/100‑1 (MQ=255)
                                                                 ttttTGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  2:2554805/100‑1 (MQ=255)
                                                                 ttttCTTTTGACCTTTACTTATACATAACAATCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  2:2937486/100‑1 (MQ=255)
                                                                 ttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:938286/100‑1 (MQ=255)
                                                                 ttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:443765/100‑1 (MQ=255)
                                                                       ttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                        <  1:635087/100‑1 (MQ=255)
                                                                        tcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                       >  1:3049080/1‑100 (MQ=255)
                                                                         cgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                       >  1:1879415/1‑99 (MQ=255)
                                                                          gACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAGTCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTgcgcgc                       >  2:250654/1‑98 (MQ=255)
                                                                            ccTTTGCTTATACATAACAGTCGATTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTGCGCGCCaa                    >  2:343319/1‑100 (MQ=255)
                                                                            ccTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                    >  2:2094048/1‑99 (MQ=255)
                                                                            ccTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                   >  2:766706/1‑100 (MQ=255)
                                                                             cTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                  <  2:2694290/100‑1 (MQ=255)
                                                                             cTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAACACTAAGTTTAAATTGCGCACCAAt                   >  2:272936/1‑99 (MQ=255)
                                                                              tttGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                 <  2:2229412/100‑1 (MQ=255)
                                                                               ttGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAACTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                <  1:720085/100‑1 (MQ=255)
                                                                               gatcTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                <  1:2812182/97‑1 (MQ=255)
                                                                                tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg               <  1:1507308/100‑1 (MQ=255)
                                                                                tGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg               >  1:1275287/1‑100 (MQ=255)
                                                                                 gCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              <  1:1687573/100‑1 (MQ=255)
                                                                                 gCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              >  1:1747939/1‑100 (MQ=255)
                                                                                  catataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  1:2827862/3‑100 (MQ=255)
                                                                                  cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:1208466/1‑100 (MQ=255)
                                                                                  cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:855351/1‑100 (MQ=255)
                                                                                   ttATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg            >  1:1318671/1‑100 (MQ=255)
                                                                                    tataCATAACAGTCGTTTTTTT‑TTTTTTTAAATAATTCAAAGAAATCAAGACAATAGTTTTGTAAAGTATTAAGTTTAAATTGCGCGCCAATCAt                >  2:1387989/1‑95 (MQ=255)
                                                                                    tataCATAACAGTCGTTTTTTT‑ATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATTgcgc           >  2:2156838/1‑100 (MQ=255)
                                                                                        cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac       <  2:41241/100‑1 (MQ=255)
                                                                                         aTAACAGTCGTTTTTTT‑TTTTTTTTTTTAATTCAAAGAAATCAAGACAATAGCGTTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca      >  2:94126/1‑100 (MQ=255)
                                                                                          tAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  2:1529020/1‑100 (MQ=255)
                                                                                           aaCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg    >  2:2819735/1‑100 (MQ=255)
                                                                                            aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaattta   <  1:2937517/100‑6 (MQ=37)
                                                                                            aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc   <  2:2561328/100‑1 (MQ=255)
                                                                                             cAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  >  1:242305/1‑100 (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: ‑