Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 16.4% Q433K (CAA→AAA)  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,4270CA16.4% 1085.1 / 68.6 405Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (108/217);  new base A (76/3);  total (185/220)
Fisher's exact test for biased strand distribution p-value = 6.31e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 416 total aligned reads displayed.

GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470280‑2470576
                                                                                                                                                   |                                                                                                                                                     
gtagtaTCTCTTTTTGGAATATCTCAATTTTATGGGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaa                                                                                                                                                     <  1:1417633/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATACCAGTCGTTTTTTTAATTTATAAAATAATTaaaa                                                                                                                                                    >  1:7621017/1‑150 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTTATTTTAAAATAATTaaaa                                                                                                                                                    >  1:3683460/1‑150 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAATTCaaa                                                                                                                                                    >  1:6276235/1‑150 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  2:2939626/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:996259/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:2792507/150‑1 (MQ=255)
    taTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCACGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaa                                                                                                                                                 <  2:3949323/150‑1 (MQ=255)
    taTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaa                                                                                                                                                 <  1:6768287/150‑1 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  1:4525408/150‑1 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  1:8478972/150‑1 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                                                                                               >  2:6642127/1‑150 (MQ=255)
       ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGaaaaa                                                                                                                                              >  2:6968576/1‑148 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACACAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                             <  1:6686578/150‑1 (MQ=255)
         ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                            <  1:4997479/150‑1 (MQ=255)
         ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                            <  2:4820440/150‑1 (MQ=255)
            tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                                                                                                         <  1:7148627/150‑1 (MQ=255)
              tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAAACaa                                                                                                                                       >  2:4600983/1‑150 (MQ=255)
               ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAt                                                                                                                                      >  1:5731339/1‑150 (MQ=255)
                 aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAg                                                                                                                                    >  1:9543051/1‑150 (MQ=255)
                     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  2:7898075/150‑1 (MQ=255)
                     actcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGACATTTAAGCAAGAATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  1:2623405/149‑1 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                              <  1:6519770/150‑1 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                              >  2:5825143/1‑150 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGt                                                                                                                              >  2:8243362/1‑150 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAAAAAAGACAAAATAGATGt                                                                                                                              >  1:3393797/1‑150 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTATTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                              >  2:4240691/1‑150 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTGTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAAATCAAAGAAAAAAAGACAAAAGCAATGt                                                                                                                              >  2:9352787/1‑150 (MQ=255)
                         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                                                                                            <  2:4673656/150‑1 (MQ=255)
                         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAGCGATTTaa                                                                                                                            >  2:3288789/1‑150 (MQ=255)
                          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATTTaaa                                                                                                                           >  1:1666269/1‑150 (MQ=255)
                          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           >  1:4818057/1‑150 (MQ=255)
                             ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  1:6479751/150‑1 (MQ=255)
                             ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  1:7932702/150‑1 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                 <  1:3800800/150‑1 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAAAAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                 >  1:5124877/1‑150 (MQ=255)
                                     cATTTTTCGGTGATGGATATAGAGATTTAAGCAAGCATTTATTTGGATTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCATAGAAATCAAGACAATAGCGATGTAAAGCTCTAAGttt                                                                                                                <  1:963355/150‑1 (MQ=255)
                                     cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                <  1:2338884/150‑1 (MQ=255)
                                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  2:1157926/150‑1 (MQ=255)
                                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  1:5483270/150‑1 (MQ=255)
                                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                                                                                           <  2:1285243/150‑1 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                        <  2:6058023/150‑1 (MQ=255)
                                              gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                                                                                                       <  1:8190373/150‑1 (MQ=255)
                                                     atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAAAAGCGATGTAAAGCAATAATTTTAAATTTCGCGCCAAta                                                                                                >  2:5469895/1‑149 (MQ=255)
                                                         agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAAATTAGCGCCAATCATgg                                                                                            >  2:1708948/1‑150 (MQ=255)
                                                          gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTTCGCGCCAATCATGgc                                                                                           >  2:5293308/1‑150 (MQ=255)
                                                          gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTTTAAAATAATTAAAAGAAATAAAGAAAAAAACGAAGTAAAGAACAAAGTTTTAATTTCGGGCCCAGCATGgc                                                                                           >  1:6077747/1‑150 (MQ=255)
                                                           agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAATAAAAACAAAAGCGATGTTAATAACGAAGTTTTAATTTGGGGCCAAACATgggg                                                                                          >  1:7338499/1‑148 (MQ=255)
                                                           agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTTTAAAATAATTAAAAGAAATAAATACAATAGAGATTTAAAGAAATAAGTTTTAATTGAGGGCCCATCATGgcg                                                                                          >  1:6895230/1‑150 (MQ=255)
                                                           agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg                                                                                          <  1:7866184/150‑1 (MQ=255)
                                                                 aatcaagctTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaat                                                                                    <  2:6303002/147‑3 (MQ=255)
                                                                 aagcaatcTTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaat                                                                                    <  2:2291653/150‑3 (MQ=255)
                                                                  agcaagaaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                   <  2:1820788/150‑1 (MQ=255)
                                                                     aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAATAAAGACAAAAGAGATGTAAAGCACAAAGTTTAAATTTCGCGCCAATCAAGGCGCGCACAAGcg                                                                                >  2:267731/1‑149 (MQ=255)
                                                                         atttatttGGAATGTATTTTTCGTTCGACCTTTTCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                            <  2:3818710/150‑1 (MQ=255)
                                                                         atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGATataa                                                                            >  2:10033056/1‑150 (MQ=255)
                                                                             atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                        <  2:9642588/150‑1 (MQ=255)
                                                                             atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAAAAATAAAGACAAAAGGGTTGTTAATAAATAATTTTTAAATTAGCGCCAAAAAAGGGGCGCCCAAGCAaaaataaa                                                                        >  1:2243359/1‑148 (MQ=255)
                                                                              tttGGAATTTATTTTTCGTTCGACCTTTTATTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACATTTGCTATGTAAAGCAATAAGTTTAATTTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCa                                                                       <  1:1192924/150‑1 (MQ=255)
                                                                              tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGAAAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATAAATGCGCGCACAAGCTataataaaa                                                                       >  2:8792843/1‑147 (MQ=255)
                                                                                        aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                             <  2:5186718/150‑1 (MQ=255)
                                                                                               gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAAACAAGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcca                                                      >  2:7491073/1‑149 (MQ=255)
                                                                                                ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAATCAAGAAAAAAGCGATGGTAATCACAAAATTTAAATATCGCGCCAATCATCGCCCGCACCAGCTAGAAAACCAACCAAAATTatactcctc                                                     >  2:213895/1‑150 (MQ=255)
                                                                                                ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTTCGCGCCAATAATGGCGAGCAAAAGCTATAATAACAAACAAATTTcacctcctc                                                     >  2:4063556/1‑150 (MQ=255)
                                                                                                 tctacttttGCTTATACATAACTGTTGTTTTTTTAATTTATAAAATAATACAAATAAATAAATACAATAGCTATGTAAAGAACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                    <  1:3288789/144‑1 (MQ=255)
                                                                                                 tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTAGCGCCAATCATGGCGCGCACAAGCTATAATACCAACATAATTTATACTActct                                                    >  2:6233003/1‑150 (MQ=255)
                                                                                                   gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                  <  2:1489529/150‑1 (MQ=255)
                                                                                                   gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                  <  1:10026911/150‑1 (MQ=255)
                                                                                                   gACCTTTGCTTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                  <  1:3856627/150‑1 (MQ=255)
                                                                                                      cTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTagag                                               <  2:3683460/150‑1 (MQ=255)
                                                                                                           cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                          <  2:4416613/150‑1 (MQ=255)
                                                                                                           cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                          <  2:3011908/150‑1 (MQ=255)
                                                                                                            ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACCAGCTCTAATACCAACCTAATTTCTCCTCCTCGTAGAGTGACta                                         >  2:9090889/1‑150 (MQ=255)
                                                                                                            ttATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAATAAATAAAAAAAAAAGCGATGTAAATGACTAAGTTTAAATTGAGCGCCAAAAATTGCGCGCACAAGCTATAAAAACAAACAAAATTGAACTCCTCTTAGAGTGACta                                         >  2:8540692/1‑150 (MQ=255)
                                                                                                                   tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                  <  1:841498/150‑1 (MQ=255)
                                                                                                                   tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                  <  2:4993876/150‑1 (MQ=255)
                                                                                                                       tttcttttttttAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCTATGTAAAGCACTTTGTTTAAATTTCGCGCCACTCATGTCGCGCACAAGCTATAATACCATCCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                              <  1:8959889/145‑1 (MQ=255)
                                                                                                                        gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                             <  2:2508420/150‑1 (MQ=255)
                                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTCATTTCTCCTCCTCTTAGAGTGACTATATATCCTGATAGAAt                         >  1:9263079/1‑150 (MQ=255)
                                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                         >  2:8131514/1‑150 (MQ=255)
                                                                                                                              tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       >  1:9444146/1‑150 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  2:3118204/150‑1 (MQ=255)
                                                                                                                                 ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                    >  1:8328811/1‑150 (MQ=255)
                                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  1:1762317/150‑1 (MQ=255)
                                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  2:7327874/150‑1 (MQ=255)
                                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              >  1:2088376/1‑150 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATCGCGATGTCAAGCACTAAGTTGAAATTGCGCGCCAATCATGGCGCGCACGCGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATCTCTCCTGATAGACGGGCGGTATTgcg             >  2:3420194/1‑148 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  2:7548712/1‑150 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  2:6152501/150‑1 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  1:5264637/1‑150 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  1:8780341/1‑150 (MQ=255)
                                                                                                                                           aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  1:8846633/1‑150 (MQ=255)
                                                                                                                                           aaaTAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          <  2:6863268/150‑1 (MQ=255)
                                                                                                                                             ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATGATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  1:5942899/1‑150 (MQ=255)
                                                                                                                                             ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        <  2:5517691/150‑1 (MQ=255)
                                                                                                                                               aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  2:7044030/150‑1 (MQ=255)
                                                                                                                                                atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt     <  2:8090450/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    >  1:4381155/1‑150 (MQ=255)
                                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc   >  1:1891443/1‑150 (MQ=255)
                                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  <  2:351468/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470280‑2470576

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: