Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,431 G→T 11.0% R434I (AGA→ATA)  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,4310GT11.0% 718.3 / 5.7 252R434I (AGA→ATA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base G (106/107);  new base T (22/3);  total (142/110)
Fisher's exact test for biased strand distribution p-value = 2.13e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

GTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGA  >  NC_000913/2470283‑2470580
                                                                                                                                                    |                                                                                                                                                     
gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGa                                                                                                                                                      >  1:2211012/1‑150 (MQ=255)
     tctTGTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTCAAAGAAATCa                                                                                                                                                 <  2:1155081/150‑1 (MQ=255)
      ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAATCAAAAAAAACaa                                                                                                                                                >  2:111485/1‑150 (MQ=255)
       tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAAAAAATAAAg                                                                                                                                               >  2:354063/1‑150 (MQ=255)
       tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATaaat                                                                                                                                               >  2:551364/1‑149 (MQ=255)
         tttGGAATATCTCAATGTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                                                                                                             <  1:3396036/150‑1 (MQ=255)
           tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGACGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACaa                                                                                                                                           >  1:2471410/1‑150 (MQ=255)
           tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACaa                                                                                                                                           <  2:2559379/150‑1 (MQ=255)
            ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACaaa                                                                                                                                          >  1:1531131/1‑149 (MQ=255)
            ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAATTAAAAAAAAAAAACACaaa                                                                                                                                          >  1:3552322/1‑149 (MQ=255)
              aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                                                                                                        <  1:1974631/150‑1 (MQ=255)
                  tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                    <  2:1058762/150‑1 (MQ=255)
                    tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACACAATCGATGt                                                                                                                                  >  2:2375162/1‑150 (MQ=255)
                    tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAAAAGCGATGt                                                                                                                                  >  1:498995/1‑150 (MQ=255)
                         tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                                                                                                             <  2:3332919/150‑1 (MQ=255)
                              gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                        <  1:3517991/150‑1 (MQ=255)
                                gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGAGATGTAAAGCACTAAGt                                                                                                                      >  2:961314/1‑150 (MQ=255)
                                gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAAAATCGAAGTAAAGCACTaatt                                                                                                                      >  1:2762221/1‑148 (MQ=255)
                                 tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAATACAAAAGCGATGTAAAGCACTAAGtt                                                                                                                     >  1:566103/1‑150 (MQ=255)
                                   aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                                                                                                   <  2:3470911/150‑1 (MQ=255)
                                     ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTATAAAAAAAATAAAAGAAATAAAGAGAATAGGGATGTGAAGGACTTAGtgttaa                                                                                                                 >  1:737735/1‑145 (MQ=255)
                                           gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCGTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATCATTAAAAAAACGAGGTACAGCACTAAGTTTAAAATgcgc                                                                                                           >  2:2705369/1‑150 (MQ=255)
                                             gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                         <  2:2445039/150‑1 (MQ=255)
                                                ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAAAAAAAATCGATTTAAAACACAAAGTTTTATTTTGGCGCCaa                                                                                                      >  1:2656125/1‑150 (MQ=255)
                                                ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAAAAAAAAATAAAAAAGAGATGGAAAGCACAAAGTTTTAATTTAGCGccca                                                                                                      >  1:1084248/1‑148 (MQ=255)
                                                 gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAAACAACAAAAAACCCATGTTAATCAAGAAGTTTTAATTTCGCGCCaaa                                                                                                     >  1:110816/1‑149 (MQ=255)
                                                   tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                                                                                                   >  2:2284216/1‑150 (MQ=255)
                                                    ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                                                                  <  2:1535889/150‑1 (MQ=255)
                                                     taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                 <  1:419422/150‑1 (MQ=255)
                                                          aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACAAAGTTTAAATTTCGCGCCAAACATGgcgcg                                                                                            >  2:1601271/1‑150 (MQ=255)
                                                             tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                         <  2:3715150/150‑2 (MQ=255)
                                                               agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAAAAAAAAAAACAAAAACGATGTAAATAACAAAGTTTAAATTTGGCGCCCATCATGGCGCGCCCaa                                                                                       >  1:402958/1‑150 (MQ=255)
                                                               agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAATAAATAAAGACAAAATAGATGTTAAACACAAATTTTTAATTTCGAGCCAAACATGGCGCGCACaa                                                                                       >  1:2833669/1‑150 (MQ=255)
                                                                gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTAAAAAATAATTAAAAGAAATAAATAAAATAGAGATGTAAAGCACAAAGTTTAAATTTCGCGCCAATAATGGCGCGCACAAg                                                                                      >  2:1286052/1‑150 (MQ=255)
                                                                gcaagcaTTTATTTGGAATGTATCTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg                                                                                      <  2:3240213/150‑1 (MQ=255)
                                                                 caagcaTTTATTTGTAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                     <  2:338463/150‑1 (MQ=255)
                                                                 caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                     <  2:1318407/150‑1 (MQ=255)
                                                                   agcaTTTATTTGGAATGTATTTTTCGTTCGACATTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                   <  2:3389244/150‑1 (MQ=255)
                                                                      atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                                <  2:491222/150‑1 (MQ=255)
                                                                          ttctGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATAGCGATGTAAATCACAAAGTTTAAATTGCGCGCCAATCATGGCGCGCACCAGCAATAAAAcc                                                                            >  2:3885975/4‑150 (MQ=255)
                                                                           tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGAGCGCCAAACATGGCGCGCACAAGCGATAATACCa                                                                           >  2:2957054/1‑150 (MQ=255)
                                                                           tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGAAAAAATCCATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCAAGGCGCGCACAAGCAAAAATAcaa                                                                           >  2:1921258/1‑148 (MQ=255)
                                                                             tGGAATGTATTTTTCGTTCTCCCTTGGCTAATACATAACAGCCTTTTTTTTAATTAAAAAAATAATTCAAAGAAATCAAGAAAATAGCTTTGTAAAGCACTAATTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAAc                                                                         <  1:3704326/150‑1 (MQ=255)
                                                                                  tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAt                                                                    <  2:1852157/150‑1 (MQ=255)
                                                                                       ttttCGTTCGACCTTTGCTTATACATAACCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGTCGCGCACAAGCTATAATACCAACCTAATTtctc                                                               <  2:2923847/150‑1 (MQ=255)
                                                                                       ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                               <  2:3704055/150‑1 (MQ=255)
                                                                                        tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                              <  2:3767688/150‑1 (MQ=255)
                                                                                        tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAAGAAAAAAAAAAAAAAGCGATGTAAACCAATAACTTTAAAATGCGCGCCAATCAAGGCGCGCACAAGCTATAAAAACAAACAAATTtctac                                                              >  2:464648/1‑148 (MQ=255)
                                                                                            gttcgACCTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                          <  1:2164489/150‑1 (MQ=255)
                                                                                             ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAATGCGCGCCAATCATGGCGCGCACAACCAATAAAACCAACCAAATTTatcctcctc                                                         >  2:2748073/1‑150 (MQ=255)
                                                                                             ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATATCGATGTAAATCACTAAGTTTAAATTTCGCGCCAATCAATGCGCGCACAATCTATAAAAACAACCGAATTTatcctacac                                                         >  2:3146259/1‑146 (MQ=255)
                                                                                              tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                        <  1:2950875/150‑1 (MQ=255)
                                                                                              tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                        <  1:1654750/150‑1 (MQ=255)
                                                                                                gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                      <  2:2118133/150‑1 (MQ=255)
                                                                                                  ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                    <  2:789610/150‑1 (MQ=255)
                                                                                                      ttcTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                                <  1:2125412/148‑1 (MQ=255)
                                                                                                      tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATCTCTCCTCCTCTTAGAGTGa                                                <  1:3283608/150‑1 (MQ=255)
                                                                                                      tGCTTATACATAAAAGTCTTTTTTTTAATTTATATAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCATTCTTGGCGCGCACTAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                                <  1:2529792/150‑1 (MQ=255)
                                                                                                       gCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                               <  1:2382810/150‑1 (MQ=255)
                                                                                                       gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                               <  2:1434763/150‑1 (MQ=255)
                                                                                                       gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                               <  2:1255255/150‑1 (MQ=255)
                                                                                                        cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAAAAGCGATGTAAAGCAATAAGATTAAATTTAGCGCCAATCATGGCGCGCACAAGCTATAAAAACAAACAAAATTCACCTCCTCAAAGAGTGAca                                              >  2:2070571/1‑149 (MQ=255)
                                                                                                        cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                              <  1:3093319/150‑1 (MQ=255)
                                                                                                            taCATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGAAAATATAGATGTAAATCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATAGCAACCAAATTTCTCCTCCTCTTAGTGTGTCGatat                                          >  2:2834570/1‑150 (MQ=255)
                                                                                                              cATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct                                        <  2:1820766/150‑1 (MQ=255)
                                                                                                                 aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                     <  1:1501715/150‑1 (MQ=255)
                                                                                                                     gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACGAAGTTTAAATTGCGCGCCAATCATGGCGCGCCCAAGCTAAAATACCAACCTAATTTCTCCTACTCTAAGAGTGACTATATCTACTGATa                                 >  2:1952454/1‑150 (MQ=255)
                                                                                                                     gTCGTATTTTTAATTTAAAAAAAAAATAAAATAAAACAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGAGCCCCAATCATGGCGCGCACAACCTAAAAAACCAACCAAAAATCTCCTCAACTAAGAGTGACTATATCTCCAGATa                                 >  2:3640999/1‑150 (MQ=255)
                                                                                                                         tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                             >  1:1780931/1‑150 (MQ=255)
                                                                                                                         tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                             >  1:134016/1‑150 (MQ=255)
                                                                                                                           tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAGTCATGGCGCGCACAAGCTGTAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATGTCTCCTGATAGAATTg                           >  1:681500/1‑150 (MQ=255)
                                                                                                                           tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                           >  2:2145175/1‑150 (MQ=255)
                                                                                                                            ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                          >  1:2220538/1‑150 (MQ=255)
                                                                                                                            ttttAATTGATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                          <  2:3399804/150‑1 (MQ=255)
                                                                                                                             tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTACTTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                         <  2:1710693/150‑1 (MQ=255)
                                                                                                                             tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                         >  2:578368/1‑150 (MQ=255)
                                                                                                                             tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                         >  2:3827286/1‑150 (MQ=255)
                                                                                                                              ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                        <  1:3474460/150‑1 (MQ=255)
                                                                                                                                 aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                     >  2:198210/1‑150 (MQ=255)
                                                                                                                                 aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCCCCTGATAGAATTGCGGTAt                     >  2:389061/1‑150 (MQ=255)
                                                                                                                                  tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                    <  2:3901851/150‑1 (MQ=255)
                                                                                                                                   ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                   >  1:3253359/1‑150 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  <  1:372332/150‑1 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  <  2:828172/150‑1 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  >  1:285654/1‑150 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  <  1:654759/150‑1 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  >  2:1054851/1‑150 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  >  1:2285738/1‑150 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  <  2:288529/150‑1 (MQ=255)
                                                                                                                                    tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                  >  2:2828081/1‑150 (MQ=255)
                                                                                                                                      taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                >  1:1237605/1‑150 (MQ=255)
                                                                                                                                       aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa               >  1:1459762/1‑150 (MQ=255)
                                                                                                                                        aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt              <  1:3140249/150‑1 (MQ=255)
                                                                                                                                         aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc             <  1:1187538/150‑1 (MQ=255)
                                                                                                                                         aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc             <  2:3377071/150‑1 (MQ=255)
                                                                                                                                          ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa            >  2:1367964/1‑150 (MQ=255)
                                                                                                                                                 aaaGAAATCAAGACAATAGCGATGTAATGCACTAAGTTTAAATTGCGCGCCAATGATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAATTGccg     >  1:3871932/1‑149 (MQ=255)
                                                                                                                                                 aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc     >  2:1553224/1‑150 (MQ=255)
                                                                                                                                                 aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc     >  1:324509/1‑150 (MQ=255)
                                                                                                                                                 aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc     >  1:3169334/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGA  >  NC_000913/2470283‑2470580

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