Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 19.6% 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,4270CA19.6% 792.0 / 66.7 321Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (80/169);  new base A (67/5);  total (147/174)
Fisher's exact test for biased strand distribution p-value = 1.01e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470280‑2470576
                                                                                                                                                   |                                                                                                                                                     
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:5243709/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:1826591/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:2143973/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATaaaa                                                                                                                                                    >  2:3047125/1‑150 (MQ=255)
  agtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaaat                                                                                                                                                   >  2:8052165/1‑149 (MQ=255)
   gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGa                                                                                                                                                  >  2:5348600/1‑150 (MQ=255)
   gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGa                                                                                                                                                  >  1:6118539/1‑150 (MQ=255)
   gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGa                                                                                                                                                  >  2:2491625/1‑150 (MQ=255)
    taTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGaa                                                                                                                                                 >  1:7961748/1‑150 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  2:1759032/150‑1 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAaaataaat                                                                                                                                               >  2:5818766/1‑150 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGaaaa                                                                                                                                               >  2:5655027/1‑149 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTCAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                                                                                               <  2:2296268/150‑1 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                             <  1:1555415/150‑1 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAtaa                                                                                                                                             >  1:2938222/1‑148 (MQ=255)
           ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                          <  1:4258342/150‑1 (MQ=255)
           ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                          <  2:6046691/150‑1 (MQ=255)
               ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACaaa                                                                                                                                      >  1:5446737/1‑149 (MQ=255)
                taaTATCTCATTTTTATTTGTAATTTTTCGGGGAAGTATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                                                                                                     <  1:6365169/149‑1 (MQ=255)
                gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATa                                                                                                                                     >  2:162360/1‑150 (MQ=255)
                gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAATAAATCAAGaaaaaa                                                                                                                                     >  1:334425/1‑148 (MQ=255)
                 aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAg                                                                                                                                    >  2:5380128/1‑150 (MQ=255)
                   tatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                                                                                                  <  1:1808739/150‑1 (MQ=255)
                    atCTCAATTTTATGTGTCATTTTTCGGGGAAGTATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                 <  1:2277135/150‑1 (MQ=255)
                    atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                 <  1:4635681/150‑1 (MQ=255)
                     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  1:1454566/150‑1 (MQ=255)
                           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                          <  1:4481032/150‑1 (MQ=255)
                            tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAAAAATCAAGTCAATATCGATGTAAAgg                                                                                                                         >  2:3492596/1‑149 (MQ=255)
                               aTGTGTCATTGTTCGGGGAAGGATATAGAGAATTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                                                                                                      <  2:6005037/150‑1 (MQ=255)
                                  tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAATAAAAAATCGATGTGAATCaaaaag                                                                                                                   >  2:8192058/1‑145 (MQ=255)
                                  tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAAAAAACACCAAAGCCATGTAAAACACaaat                                                                                                                   >  2:1460103/1‑149 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAAAAGCGATGTAAAGCACTAAGtt                                                                                                                 >  2:6421684/1‑150 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGtt                                                                                                                 >  1:484576/1‑150 (MQ=255)
                                     cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATATCGATGTAAAGCACTAAGttt                                                                                                                >  2:6820393/1‑150 (MQ=255)
                                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                                                                                               <  2:6545274/150‑1 (MQ=255)
                                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                                                                                               <  1:3732212/150‑1 (MQ=255)
                                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  1:7722326/150‑1 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                        <  2:2766820/150‑1 (MQ=255)
                                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                     <  1:1079428/150‑1 (MQ=255)
                                                   ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                  >  2:7089745/1‑150 (MQ=255)
                                                   ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                  >  2:7380855/1‑150 (MQ=255)
                                                   ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAATCAAGAAAAAAGCGATGTAAAGCAATAAGTTTAAATTTCGCGCCaa                                                                                                  >  2:3349805/1‑150 (MQ=255)
                                                    gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                                                 >  2:3223863/1‑150 (MQ=255)
                                                       ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAAATAAAAGAAAAAAACAAAAAAGCGAATTAAAGCCCCAAGATTAAATTGCCCGccccacca                                                                                              >  1:1261641/1‑144 (MQ=255)
                                                          gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAATAGCGATGTAAATCACAAAGTTTAAATTTCGCGCCAATCATGgc                                                                                           >  2:6664364/1‑150 (MQ=255)
                                                             aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTTCGCGCCAAACATGgcgcg                                                                                        >  2:7459072/1‑150 (MQ=255)
                                                             aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATCAAGACAAAAGAGATGTAAATCACTAAGTTTAAATTTCGCGCCAATCAAGgcgcg                                                                                        >  2:3969646/1‑150 (MQ=255)
                                                                 aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTATAAAAAAAATAAAAAAAATAAAAAAAAAATAGAAGTAAAGAACAAAGGTTAAATTTGGGGCCCAGAATGGCGCGgacc                                                                                    >  1:7104890/1‑149 (MQ=255)
                                                                    caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAGAAAAAAGAGATTTAAAGAACTAAGTTTTAATTTGGCGCCCATCAAGgcgcgccccacc                                                                                 >  2:777729/1‑144 (MQ=255)
                                                                     aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATAAAATATAGATGTAAAACACTAAGTTTAAATTTCGCGACAATCATGGCGCGCAAAACCt                                                                                >  2:8249089/1‑150 (MQ=255)
                                                                     aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAAATTATAAAATAAAAAAAAGAAAAAAAGACAAAAAAAATGTAAAGAACAAATTTTAAATTTAGCGCCCCACAAGGCGCGcaccaact                                                                                >  1:7775130/1‑145 (MQ=255)
                                                                      agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                               >  2:4023319/1‑150 (MQ=255)
                                                                      agcaTTTACTTGGAATGTATTTGTCGTTCGACCTCTCCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatttaaa                                                                               <  2:5089257/150‑8 (MQ=255)
                                                                        caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGACATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtata                                                                             <  2:3188015/150‑1 (MQ=255)
                                                                         atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTATAAAAAAAATAAAAAAAAAAAAAAAAAAAGCGAAGTAAATCACTAAGATTAAATATAGCCCCCCTCAAGGCGCGCCCaaccaaaaa                                                                            >  1:2745349/1‑143 (MQ=255)
                                                                                tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                     <  2:5510343/150‑1 (MQ=255)
                                                                                     tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAt                                                                <  1:4309295/150‑1 (MQ=255)
                                                                                         tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAAAAAGACAAAAGCGATTTAAAGCACAAAGTTTAAATTTGGCGCCAATCATGGCGCGCCCAAGCTATAATAAAAACCTAAattaa                                                            >  2:6306390/1‑148 (MQ=255)
                                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  2:2361607/150‑1 (MQ=255)
                                                                                             tcgttcgACCTTTGCTTATACATAACTGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGAAAAGCACTACGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                        <  1:979914/150‑1 (MQ=255)
                                                                                             tcgttcgACCTTTGATTATACATAACTGTCTTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                        <  1:8052165/150‑1 (MQ=255)
                                                                                             acgttcgACCTTTGATTATACATAAAAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTACTTTTAACTTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                        <  1:3648660/149‑1 (MQ=255)
                                                                                              cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  2:2039666/150‑1 (MQ=255)
                                                                                              cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  2:4879726/150‑1 (MQ=255)
                                                                                                ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAAACAAACAAAACAACACAAAAGCCATTTAAAGCACTAAGTTTAAATTTCGCGCCAAACAATGCGCCCCCAAGCAATAATACCAACCAAAATTcaccaccac                                                     >  2:5683890/1‑148 (MQ=255)
                                                                                                   gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                  <  2:4568732/150‑1 (MQ=255)
                                                                                                   gACCTGTTCTTATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATCGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                  <  1:4049582/150‑1 (MQ=255)
                                                                                                     ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                <  2:695601/150‑1 (MQ=255)
                                                                                                          gCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAATAAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTTAGCGCCAAAAATTGAGCGCACAAGCAATAAAACAAACCAAATTTAACCTACACTTAGAGTGAc                                           >  2:3127967/1‑150 (MQ=255)
                                                                                                           cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                          <  1:1192233/150‑1 (MQ=255)
                                                                                                              ataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtata                                       >  2:4711216/1‑150 (MQ=255)
                                                                                                                     aCAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAAACATTGCGCGCACAAGCTATAAAACCAACCTAAGTTATCCTCCTCTTAGAGTGACTATATATccag                                >  2:5415902/1‑148 (MQ=255)
                                                                                                                        gTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                             <  2:4289311/150‑1 (MQ=255)
                                                                                                                        gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGTTGTAAAGCACTAAGTTTAAAGTGCGCGCCAATCATGGCGCGCACAAGCTATAAGACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAtc                             <  1:1375246/150‑2 (MQ=255)
                                                                                                                         tCGTTTTTTTAATTTATAAAATAATACAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATAATGGCGCGCACAAGCTATAAAACCAACCTAATATCACCTACTCTTAGAGTGACTATATCTCATGATAg                            >  2:5183089/1‑150 (MQ=255)
                                                                                                                          cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                           <  1:6476068/150‑1 (MQ=255)
                                                                                                                           gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCCCAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          >  2:6228101/1‑150 (MQ=255)
                                                                                                                           gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          <  2:7287043/150‑1 (MQ=255)
                                                                                                                              tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       >  1:5880867/1‑150 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:6820393/150‑1 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:2585832/150‑1 (MQ=255)
                                                                                                                                   aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                  >  1:5358675/1‑150 (MQ=255)
                                                                                                                                    aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                 <  1:6501152/150‑1 (MQ=255)
                                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  1:1917338/150‑1 (MQ=255)
                                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  1:2894421/150‑1 (MQ=255)
                                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              >  2:7313405/1‑150 (MQ=255)
                                                                                                                                       tagaaaaTAATTCAAAGAACTCAAGACAATAGCGAGGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  2:7903634/147‑1 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  1:7291398/150‑1 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  1:1646051/1‑150 (MQ=255)
                                                                                                                                          aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa           >  1:6726035/1‑150 (MQ=255)
                                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         <  1:5655027/150‑1 (MQ=255)
                                                                                                                                             ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  2:114134/1‑150 (MQ=255)
                                                                                                                                              taatTCAAAGAAATCAAGACAATAGAGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa       <  2:7736112/150‑1 (MQ=255)
                                                                                                                                               aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  2:4889094/150‑1 (MQ=255)
                                                                                                                                               aatTCAAAGAAATAAAGACAATAGCGATTTAAAGCACTAAGTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  1:5348693/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    >  2:1957760/1‑150 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  2:1295959/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  1:6871716/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTATAAATAGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCAAATTTCTCATCCTCTAAGAGTGACTATATATCCTGATAGAATTGCGGTATTGACTATCAAAtc    >  2:4443678/1‑149 (MQ=255)
                                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc   >  2:3755737/1‑150 (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: