Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,441 C→A 9.7% D437E (GAC→GAA 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,4410CA 9.7% 560.8 / 8.2 175D437E (GAC→GAAgtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (53/98);  new base A (21/3);  total (74/101)
Fisher's exact test for biased strand distribution p-value = 1.79e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

GAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTT  >  NC_000913/2470296‑2470590
                                                                                                                                                 |                                                                                                                                                     
gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAGaaaaaa                                                                                                                                                   >  1:847102/1‑148 (MQ=255)
  atatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                                                                                                                 <  2:3149057/150‑1 (MQ=255)
     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                              <  2:3477838/150‑1 (MQ=255)
       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                                            <  1:1601998/150‑1 (MQ=255)
       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAGCGATGt                                                                                                                                            >  1:281687/1‑150 (MQ=255)
        cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTTTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                                           <  1:1599175/150‑1 (MQ=255)
        cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                                           <  1:3484489/150‑1 (MQ=255)
        cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTTTAAAATAATTAAAATAAATAAATAAAAAAACGATGTa                                                                                                                                           >  1:3997522/1‑150 (MQ=255)
        cAATTTTATGTGCCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAAAAAAGTCACTATCGATGTa                                                                                                                                           >  2:3593708/1‑150 (MQ=255)
         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATAAAAATAAATCAAGAAAAAATGGATGTaa                                                                                                                                          >  1:2322708/1‑150 (MQ=255)
          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                                         <  1:4240604/150‑1 (MQ=255)
           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAATAAATAAAGACAAAATCGAATTaaat                                                                                                                                        >  2:3292210/1‑149 (MQ=255)
             ttATGTGTGATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                                      <  2:3494377/150‑1 (MQ=255)
             ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCa                                                                                                                                      >  2:2438088/1‑150 (MQ=255)
                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                                  <  2:4120594/150‑1 (MQ=255)
                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                                  <  1:140162/150‑1 (MQ=255)
                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                               <  1:367588/150‑1 (MQ=255)
                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                               <  1:37409/150‑1 (MQ=255)
                     cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAATCAAAAAAAAAGCGATGTAAAGAACGAAGttt                                                                                                                              >  1:3501649/1‑150 (MQ=255)
                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGTTTATACATAACAGTCGTTTTTTTAATTTAAAAAATAATTCAAATAAATAAAGAAAAAAGCGATGTAAAGCACTAAGTTTa                                                                                                                             >  1:2073017/1‑150 (MQ=255)
                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTCAAAGAAATCAAGAAAAAAGCGATTTAAAGCACAAAGTTTa                                                                                                                             >  1:381480/1‑150 (MQ=255)
                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAAAAGCGATGTAAAGCACAAAGGTTa                                                                                                                             >  1:969000/1‑150 (MQ=255)
                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                                           <  1:22454/150‑1 (MQ=255)
                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                                                                                                         <  1:2975345/150‑1 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                                   <  1:2328009/150‑1 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAATAAAAAAGCGAAGTAAAGCACAAAGTTTAAAATtcgcgc                                                                                                                   >  2:3374241/1‑150 (MQ=255)
                                    gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                                                               >  2:4179649/1‑150 (MQ=255)
                                     atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAAAAAATACAAAATCGATGTAAAGAACAAAGTTTAAATTTCGCGCCaaac                                                                                                              >  1:326450/1‑148 (MQ=255)
                                      tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAACTGCGCGCCatgat                                                                                                             <  1:3122351/150‑5 (MQ=255)
                                      tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAACTGCGCGCCatgat                                                                                                             <  1:2060045/150‑5 (MQ=255)
                                        taGATATTTAATCAATCATTTATTTGGAATTTATTTTTCGTTAGAACTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAATATAATTCAAAGAAATCAAGACAATTGCGATGTAATGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                           <  1:2494302/150‑1 (MQ=255)
                                        taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                           <  1:2578927/150‑1 (MQ=255)
                                          gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTTCGCGCCAAACATGgc                                                                                                         >  2:671393/1‑150 (MQ=255)
                                           agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAAGTAAAGCACAAATTTTAAATTTCGCGCCAATCATGgcg                                                                                                        >  2:3859931/1‑150 (MQ=255)
                                            gaTTTAAGCACGCATTTATTTTGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                                       <  2:99035/150‑1 (MQ=255)
                                            gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                                       <  2:819997/150‑1 (MQ=255)
                                             aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAAGAAATCAAGACAAAAGCGCTGGAAAGCCCGAAGTTTAAATTGCGCGCCCAACAAGgcgcg                                                                                                      >  2:211957/1‑150 (MQ=255)
                                               ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca                                                                                                    <  2:2282040/150‑1 (MQ=255)
                                                tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                                   <  2:481002/150‑2 (MQ=255)
                                                tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                                   <  2:1969766/150‑2 (MQ=255)
                                                 aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACGAAGTTTAAAGTGCGCGCCATGATTGGCGCGcaat                                                                                                  <  2:4032483/150‑3 (MQ=255)
                                                    caagcaTTTATTTGGAATGAATTTTTCGTTCGACCTATGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAATTCAAACAAATCAATACAATAGCGAATTAAAGCACAAATTTTAAATTGCGCGCCAATCATAGAGCCCACAAGc                                                                                               >  2:4229022/1‑150 (MQ=255)
                                                     aagcaTTTAGTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTAAAAAAAAAATAAAAGAAATCAACACCAAAGCCCTGTAAAGCACAAAGTTTAAAATGCGCCCCAATCAAGGCGCGCCCAAGca                                                                                              >  2:1256067/1‑149 (MQ=255)
                                                        caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGAAAATATCGATTTAAATAACAAAGTTTAAATTTCGCGCCAAACAAGGCGCGCACAAGCtata                                                                                           >  2:64203/1‑150 (MQ=255)
                                                        caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtata                                                                                           <  2:3183885/150‑1 (MQ=255)
                                                         atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACTTAACAGTCGTTTTTTTAATTTTTAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                                          <  2:1270837/150‑1 (MQ=255)
                                                              tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCa                                                                                     <  2:2729210/150‑1 (MQ=255)
                                                              tttGGAAGGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTCTAAAATAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCa                                                                                     <  2:2677275/150‑1 (MQ=255)
                                                                    aTGTATTTTTCGTTCGACCTTTTCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGTCGCGCACATGCTATAATCCCAACCTaa                                                                               <  2:462233/150‑1 (MQ=255)
                                                                         tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAAACAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAAACAAGGCGCGCCCAAGCGATAATAACAACCTAAattat                                                                          >  2:4008466/1‑148 (MQ=255)
                                                                             tcgttcgACCTTTTCTTATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                                      <  1:3616145/150‑1 (MQ=255)
                                                                             tcgttcgACCTTTTATTATACATACCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                                      <  1:3613646/150‑1 (MQ=255)
                                                                             tcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                                      <  2:146526/150‑1 (MQ=255)
                                                                                       tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATAAATAAAAAAGCGATTTAAAGCACTAATTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATTCCAACCAAAATTATACTACACTTagatt                                                            >  2:3570134/1‑148 (MQ=255)
                                                                                       tttGCTTATACATAACAGTCGTTTTTTTAATTTAAAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTCAAGCACTAAGTTTAAAATGCGCGCCAATCATGGCGCGCACAAGCGATAATAACAACCTAATTTCTCCTCCGCGTAGAGt                                                            >  2:136493/1‑150 (MQ=255)
                                                                                         tGCTTATACATCACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATCAAGAAAATAGCGAAGTAAAGCACTAAGTTTAAATTTAGCGCCAATCATGGCGCGCACAAGCTATAATAACAACCAAATTTATACACCTCTTAGAGTgg                                                          >  2:2325844/1‑149 (MQ=255)
                                                                                          tattATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                         <  1:3222612/148‑1 (MQ=255)
                                                                                          gCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                         <  1:2112381/150‑1 (MQ=255)
                                                                                            ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                                       <  2:1988785/150‑1 (MQ=255)
                                                                                                  aTAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtctc                                                 <  1:2864583/150‑1 (MQ=255)
                                                                                                   tAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                                <  1:742030/150‑1 (MQ=255)
                                                                                                    aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAAAAGCGATTTAAAGCACTAAGTTTAAATTTCGCGCCAATCATTGCGCGCACAAGCTATAAAAACAACCAAATTTATACACCACTTAGAGTGACTATATCTCCt                                               >  2:1297119/1‑150 (MQ=255)
                                                                                                      cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                             <  1:603337/150‑1 (MQ=255)
                                                                                                      cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                             <  1:254600/150‑1 (MQ=255)
                                                                                                       aGTCGTTTTTTTAATTTATAAAATAATTGAAAGAAATAAAGACAATAGCGATGTAAAGCAGGTAGTTTAAATTGCGCGGCGATCGTGGCGGGCGCAAGCTGTAATTCCAACCGAATTTGTGGTCCTCTTAGAGTGTCTATATGTCCTGAt                                            >  2:3328715/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAGGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  2:2185182/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  1:1673156/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  2:2230284/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  1:376395/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  2:1055057/1‑150 (MQ=255)
                                                                                                               ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTGATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                    <  1:1282900/150‑1 (MQ=255)
                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                   <  2:3526050/150‑1 (MQ=255)
                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                   >  2:364312/1‑150 (MQ=255)
                                                                                                                 ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGAGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                                  >  2:1336158/1‑150 (MQ=255)
                                                                                                                  tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGAGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                                 <  2:1530959/150‑1 (MQ=255)
                                                                                                                    aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                               >  1:1542506/1‑150 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                            >  1:3876831/1‑150 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                            >  2:197415/1‑150 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                            >  2:1780357/1‑150 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                            >  2:3891863/1‑150 (MQ=255)
                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                           <  1:102206/150‑1 (MQ=255)
                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                       >  2:2027164/1‑150 (MQ=255)
                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                 >  2:4247764/1‑150 (MQ=255)
                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                <  2:2187302/150‑1 (MQ=255)
                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                <  2:1075960/150‑1 (MQ=255)
                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAGTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                <  2:2116640/150‑1 (MQ=255)
                                                                                                                                      aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATTCCCTg             >  2:3029014/1‑150 (MQ=255)
                                                                                                                                       taaaTCAAGACAATAGCTATTTAAAGCACTTAGTTTAAATTGCGCGCCAATCTTGTCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa            <  1:1386311/149‑1 (MQ=255)
                                                                                                                                        aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt           <  1:1292304/150‑1 (MQ=255)
                                                                                                                                         aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcccctgat          >  2:2378955/1‑145 (MQ=255)
                                                                                                                                         aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAACTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATTGATTTGCGGTATTGACTATCAAATGCCCTGAtt          <  1:3888325/150‑1 (MQ=255)
                                                                                                                                         aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCGATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt          <  1:1969889/150‑1 (MQ=255)
                                                                                                                                           tCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg        <  1:1055057/150‑1 (MQ=255)
                                                                                                                                             aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt      <  1:866200/150‑1 (MQ=255)
                                                                                                                                             aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt      >  2:2911354/1‑150 (MQ=255)
                                                                                                                                             aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt      >  2:772083/1‑150 (MQ=255)
                                                                                                                                              aGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttg     >  1:2775379/1‑150 (MQ=255)
                                                                                                                                               gACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgt    <  1:1855797/150‑1 (MQ=255)
                                                                                                                                                 cAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGttt  <  1:1476361/150‑1 (MQ=255)
                                                                                                                                                 cAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGttt  <  1:1126328/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTT  >  NC_000913/2470296‑2470590

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: