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.97.8% 611.4 / 4.4 138coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base . (63/72);  total (64/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470306‑2470498
                                                                                                  |                                                                                              
aTTTTATGTGTCATTTTTCGGGGAATTATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtt                                                                                               <  1:371052/100‑1 (MQ=255)
    tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                           <  2:1210327/100‑1 (MQ=255)
       gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        <  1:3501761/100‑4 (MQ=255)
       gtgtCATTTTTCGGGGAAGGAAATAGAGATTTAAGCAAGCATTTATTTGTAATGTATTTTTCGTTCTACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        >  2:807439/1‑97 (MQ=255)
        tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  1:1625231/100‑5 (MQ=255)
        tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  2:379303/100‑5 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     >  1:1728157/1‑100 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     >  1:568547/1‑100 (MQ=255)
          tCATTTTTCGGGGAAGGATATAGAGCTTTAAGCAAGAATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTta                                                                                    <  2:3735990/100‑1 (MQ=255)
          tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTta                                                                                    <  2:2879988/100‑1 (MQ=255)
            aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  <  2:1600531/100‑1 (MQ=255)
            aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  >  2:2964986/1‑100 (MQ=255)
              ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                <  2:890710/100‑1 (MQ=255)
                   ggggCAGGATATAGAGATTTAAGCAAGAATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           >  2:1556973/1‑100 (MQ=255)
                   ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                           >  1:1257632/1‑100 (MQ=255)
                    gggAAGGATATAGAGATTTAAGCTAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                          <  1:3160335/100‑1 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         >  1:1977143/1‑100 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                         <  2:1806246/100‑1 (MQ=255)
                       aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaa                                                                       <  1:3001533/100‑1 (MQ=255)
                       aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaa                                                                       <  2:272081/100‑1 (MQ=255)
                       aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaa                                                                       <  2:1573377/100‑1 (MQ=255)
                        aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                      >  2:204206/1‑100 (MQ=255)
                          gATATAGATATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCTACCTTTGCTTATACATAACATTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     <  2:2802104/99‑1 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     <  2:3255623/99‑1 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     <  1:727889/99‑1 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                    >  1:3910500/1‑100 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                    >  1:2812103/1‑100 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                    >  1:1852656/1‑100 (MQ=255)
                            tataGATATTTAAACAAGCATTTATTTTGAATATATTTTTCGTTCTACCTTTTCTTTTACATAACAGTCTTTTTTTAATTTATAAAATAATTAAAAGa                                                                    <  2:1646854/98‑1 (MQ=255)
                            tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  >  1:912754/1‑100 (MQ=255)
                            tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  <  1:119708/100‑1 (MQ=255)
                               agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  2:2937856/1‑100 (MQ=255)
                               agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               >  2:1244210/1‑100 (MQ=255)
                               agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                               <  1:1461909/100‑1 (MQ=255)
                                 agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                             <  1:3122330/100‑1 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  2:616810/1‑100 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            <  2:819697/100‑1 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  1:1125218/1‑100 (MQ=255)
                                   aTTTAAGCAAGCATTTATTTGGCATGTATTTGTCGTTGGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  2:2646242/100‑1 (MQ=255)
                                   aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  1:2783636/100‑1 (MQ=255)
                                    tttGAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAGCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACa                                                          <  2:917032/100‑1 (MQ=255)
                                      tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                        <  1:970811/100‑1 (MQ=255)
                                         gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     <  2:3172551/100‑1 (MQ=255)
                                         gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     <  2:509108/100‑1 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    <  1:2877075/100‑1 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  1:3705596/1‑100 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAAACAATAGCg                                                    >  1:3594747/1‑100 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAATCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  1:46064/1‑100 (MQ=255)
                                           aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAATAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                   <  1:4911/100‑1 (MQ=255)
                                            agcaTTTATTTGTAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                  <  2:1618530/100‑1 (MQ=255)
                                               atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               >  2:272355/1‑100 (MQ=255)
                                                  tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                            >  1:2425019/1‑100 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCTATCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           <  1:186824/100‑1 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           <  2:433617/100‑1 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           <  2:596784/100‑1 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           <  2:1367758/100‑1 (MQ=255)
                                                     ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                         <  2:701874/100‑1 (MQ=255)
                                                      tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                        >  1:2624526/1‑100 (MQ=255)
                                                         aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     >  1:3834/1‑100 (MQ=255)
                                                         aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     <  1:111035/100‑1 (MQ=255)
                                                          aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAACACTAAGt                                    <  2:1480521/100‑1 (MQ=255)
                                                             tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 <  1:2882889/100‑1 (MQ=255)
                                                              aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                >  2:3320618/1‑100 (MQ=255)
                                                                ttccgatctgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:1561266/91‑1 (MQ=255)
                                                                  ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                            >  2:2188622/1‑100 (MQ=255)
                                                                   tcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                           >  1:669569/1‑100 (MQ=255)
                                                                    cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                             >  2:1434281/1‑97 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtttctc                         >  1:2704760/1‑96 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         >  2:2207917/1‑100 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         <  2:2369526/100‑1 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         >  1:458942/1‑100 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATATAAAACACTAAGTTTAAATTgcgc                         >  1:2676693/1‑100 (MQ=255)
                                                                       tcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                       <  1:1532391/100‑1 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                      >  1:3736763/1‑100 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                      <  1:1185799/100‑1 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                      >  1:1625215/1‑100 (MQ=255)
                                                                               tGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg               >  2:3760394/1‑100 (MQ=255)
                                                                                gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                >  2:3092355/1‑98 (MQ=255)
                                                                                gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              >  2:2446928/1‑100 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  1:1174838/1‑100 (MQ=255)
                                                                                   tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCatgattggcgc           <  1:54511/100‑11 (MQ=255)
                                                                                   tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc           >  1:2491134/1‑100 (MQ=255)
                                                                                   tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc           >  1:2953223/1‑100 (MQ=255)
                                                                                    ataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg          <  2:282159/100‑1 (MQ=255)
                                                                                    ataCATAACAGTCGTTTTTTAATTTAGAAAATTATTTAAAGCAATAACTACGATTGCGATGTAAAGCACTAAGTTTAGTTTGCTCGCCAATCATGgcgcg          <  2:2042283/100‑1 (MQ=255)
                                                                                     taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         <  1:2954847/100‑1 (MQ=255)
                                                                                     taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         >  2:412956/1‑100 (MQ=255)
                                                                                     taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         <  1:494819/100‑1 (MQ=255)
                                                                                      aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca        >  1:3438207/1‑100 (MQ=255)
                                                                                      aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca        <  2:1932917/100‑1 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac       >  1:248488/1‑100 (MQ=255)
                                                                                        aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca      <  1:1177305/100‑1 (MQ=255)
                                                                                         tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  1:3860117/1‑100 (MQ=255)
                                                                                         tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     >  1:3501490/1‑100 (MQ=255)
                                                                                         tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATCTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCCCACaa     >  2:1854268/1‑100 (MQ=255)
                                                                                          aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg    <  1:202127/100‑1 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAATATAATTTAAAGAAATCAAGATAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc   <  1:360316/100‑1 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAAGATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     <  2:1790469/98‑1 (MQ=255)
                                                                                            cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  >  1:1727107/1‑100 (MQ=255)
                                                                                            cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  <  2:2634862/100‑1 (MQ=255)
                                                                                                  |                                                                                              
ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470306‑2470498

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