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.81.7% 288.7 / 82.0 104coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (8/11);  new base . (40/45);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 115 total aligned reads displayed.

TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTAT  >  NC_000913/2470307‑2470500
                                                                                                 |                                                                                                
ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttt                                                                                                <  1:1398766/100‑1 (MQ=255)
 tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                               <  2:2177359/100‑1 (MQ=255)
  ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                              <  2:926861/100‑1 (MQ=255)
  ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                              <  2:2078910/100‑1 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                            <  1:1552360/100‑2 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTAAACATAACAGTCGtttttta                                                                                            <  2:2113713/100‑2 (MQ=255)
      gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                          >  1:1195545/1‑97 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                         <  1:170028/100‑5 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                         >  2:923610/1‑96 (MQ=255)
        ntCATTTTTCGGGGAATGATATAGAGATTTAAGCAACCCTTTATTGCGAATGTATTTTTCGTTCGCCCTTTGCTTATACATAACAGTCGttttttacttt                                                                                        >  2:1442659/2‑95 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                       <  2:2349713/100‑1 (MQ=255)
         tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTTGTTTTTTAATTta                                                                                      >  2:2378411/1‑100 (MQ=255)
         tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTta                                                                                      >  2:1339394/1‑100 (MQ=255)
            tttttCGGTGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                   <  2:1311409/100‑1 (MQ=255)
             ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                  <  2:2054231/100‑1 (MQ=255)
              tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaaa                                                                                 <  1:687694/100‑1 (MQ=255)
               ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaat                                                                                <  1:835965/100‑1 (MQ=255)
                  ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                             >  1:2023827/1‑100 (MQ=255)
                  ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                             >  1:1943350/1‑100 (MQ=255)
                  ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                             >  2:2159385/1‑100 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAATATAATTc                                                                           <  1:1371177/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                           <  1:352476/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                           <  1:2312894/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                           <  1:2119853/100‑1 (MQ=255)
                     gAAGGTTATATAGATTTAAGCAAGCTTTTATTTGGAATGTATTTTTCGTTAGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                            <  2:1256627/98‑1 (MQ=255)
                       aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                        >  1:699690/1‑100 (MQ=255)
                        ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                       >  1:2416526/1‑100 (MQ=255)
                        ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                       <  1:2352179/100‑1 (MQ=255)
                           tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                    >  1:357416/1‑100 (MQ=255)
                           tataAAAATTTCAGCAGTCCTTTATTTTTCATGTATTTTTCGTTCGACCTATGTTTATACAGGACAGTCGTTTTTTAATTTATAAAATAAtt                                                                            <  2:2277572/92‑1 (MQ=255)
                              agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                 <  1:806612/100‑1 (MQ=255)
                              agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                 >  1:1560788/1‑100 (MQ=255)
                               gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTTTTTTTTCTTATACATAACAGTCGTTTTTTAATTTTTATTATAGTTCAAAGAAAt                                                                   >  2:2177752/1‑97 (MQ=255)
                               gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                >  1:2242264/1‑100 (MQ=255)
                               gagaTATAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                 <  2:751523/99‑1 (MQ=255)
                                agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                               <  1:1771797/100‑1 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                               <  2:1806853/99‑1 (MQ=255)
                                  aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                             <  1:2474270/100‑1 (MQ=255)
                                       agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                        <  1:309919/100‑1 (MQ=255)
                                       agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                        <  1:1455263/100‑1 (MQ=255)
                                            gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATg                                                   >  2:1339051/1‑100 (MQ=255)
                                             caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCCAAGAAATCAAGACAATCGCGATGt                                                  >  2:584385/1‑100 (MQ=255)
                                              atttatttGGAATGTATTTTTTTTTTTTTCTTTTCTTATACATAACAGTCGTTTTTTAATTTATAATATAATTCAAAGAAATCAAGACATTATCGATGt                                                  >  2:894401/1‑99 (MQ=255)
                                              atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                 <  1:1884631/100‑1 (MQ=255)
                                               tttatttGGAATGTATTTTTTTTTCTTCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                  >  2:837984/1‑98 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                <  1:1446983/100‑1 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                <  1:1593320/100‑1 (MQ=255)
                                               tttatttGAAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGATTATTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                <  2:1457659/100‑1 (MQ=255)
                                                 tatttGGTATATATATTTCGTTCGACATTTGTTTAAACATAACAGTTGTTTTTTAATTTATAAAATAATTCAAAGAATTCAAGACAATAGCGATGTAAAg                                              <  2:916375/100‑1 (MQ=255)
                                                 tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                              >  1:1761182/1‑100 (MQ=255)
                                                 tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                              >  1:343325/1‑100 (MQ=255)
                                                 tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                              <  2:1761173/100‑1 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                             <  1:923615/100‑1 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                             <  2:2183545/100‑1 (MQ=255)
                                                   tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                            >  1:856119/1‑100 (MQ=255)
                                                    ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                           >  1:586382/1‑100 (MQ=255)
                                                       gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                        <  1:2514371/100‑1 (MQ=255)
                                                        aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                       <  2:1195540/100‑1 (MQ=255)
                                                           gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                    >  1:1811337/1‑100 (MQ=255)
                                                            tATTTTTCGTTCTATTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                   >  2:43657/1‑100 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                   >  1:1154309/1‑100 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCACTAAGtt                                     >  2:1823962/1‑98 (MQ=255)
                                                              tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                 <  2:199025/100‑1 (MQ=255)
                                                               ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                                <  1:2385916/100‑1 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCG                            >  1:2350227‑M1/1‑100 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCG                            >  2:539422‑M1/1‑100 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCG                            <  1:1060307‑M1/100‑1 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCGTTTTTTATTTTTTTTAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGC                           >  2:2318249‑M1/1‑100 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGC                           >  1:1832577‑M1/1‑100 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGC                           >  1:55078‑M1/1‑100 (MQ=255)
                                                                     ttagacgTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGC                           <  2:1369543‑M1/96‑1 (MQ=255)
                                                                         aCCTTTGCTTATAGATAACAGTCGTTTTTTAATTTATATAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                      >  2:1607784/1‑100 (MQ=255)
                                                                           cTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                    >  2:201353/1‑100 (MQ=255)
                                                                            tttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAtgat                   >  2:1636639‑M1/1‑96 (MQ=255)
                                                                            tttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                   >  1:173863/1‑100 (MQ=255)
                                                                                cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGATAATATCGATGTAAAGGACTAACTTTAAATTGAGCGCCAATCAt                  <  2:702139‑M2/97‑2 (MQ=255)
                                                                                cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggc               <  2:1145675‑M2/100‑5 (MQ=255)
                                                                                cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtgg                <  2:368539‑M2/99‑4 (MQ=255)
                                                                                 ttATACAGAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGCTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtgg                >  2:1385191‑M2/1‑95 (MQ=255)
                                                                                  tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggcgc             <  1:1685164‑M2/100‑7 (MQ=255)
                                                                                  tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggcgc             >  1:1119403‑M2/1‑94 (MQ=255)
                                                                                  tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggcg              <  2:1627327‑M2/99‑6 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggcgcgc           <  1:1110184‑M2/100‑9 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggcgcgc           <  1:1141869‑M2/100‑9 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTAATTTATAAAATAAGTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTACGTTTAAATTGCGCGCCAATCAtggcgcgc           <  2:1565769‑M2/100‑9 (MQ=255)
                                                                                     aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAtggcgcgca          <  2:36146‑M2/100‑10 (MQ=255)
                                                                                      cATAACATTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac         >  1:1368474/1‑100 (MQ=255)
                                                                                      cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTCAGTTTAAATTGCGCGCCAATCAtagcgcgc           <  2:790777‑M2/98‑9 (MQ=255)
                                                                                      cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac         >  1:639479/1‑100 (MQ=255)
                                                                                      cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTACGCGCCAATCATGGCGCGcac         >  2:1233977/1‑100 (MQ=255)
                                                                                       aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGAttggcgcgcaat       >  1:1262876‑M2/1‑88 (MQ=255)
                                                                                       aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGAttggcgcgcaat       >  1:321235‑M2/1‑88 (MQ=255)
                                                                                       aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca        <  1:160234/100‑1 (MQ=255)
                                                                                       aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca        <  1:1546723/100‑1 (MQ=255)
                                                                                         aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGAttggcgcgcaattt     >  1:1366350‑M2/1‑86 (MQ=255)
                                                                                          aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     <  1:1627336/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    <  1:1073666/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    >  1:1675317/1‑100 (MQ=255)
                                                                                            aGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGAttggcgcgcaatttaaa  >  2:1743187‑M2/1‑83 (MQ=255)
                                                                                            aGTCGTTTTTTATGTTATAGAATAAGTCAAAGACAGCAAGACAATAGCGATGTAAAGCGCGAAGCTTAAATTGAACGCCAATCATGGCGCGAACAAGCt    <  2:1282769/99‑1 (MQ=255)
                                                                                                 |                                                                                                
TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTAT  >  NC_000913/2470307‑2470500

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: