Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,472 G→T 7.7% intergenic (+10/‑343) gtrS → / → tfaS serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage/pseudogene, CPS‑53 (KpLE1) prophage; tail fiber assembly protein fragment;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4720GT 7.7% 1478.0 / 30.7 428intergenic (+10/‑343)gtrS/tfaSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage/pseudogene, CPS‑53 (KpLE1) prophage; tail fiber assembly protein fragment;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (159/226);  new base T (41/2);  total (200/228)
Fisher's exact test for biased strand distribution p-value = 1.03e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 437 total aligned reads displayed.

TCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGAT  >  NC_000913/2470323‑2470620
                                                                                                                                                     |                                                                                                                                                    
     gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                                                                 <  1:892754/150‑1 (MQ=255)
       aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                                                               <  2:6203574/150‑1 (MQ=255)
         gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAAGAAATCAAGACAAAAGCGATGGAAAGCACTAAGTTTAAATTTCGCGCCAAt                                                                                                                                             >  2:521913/1‑150 (MQ=255)
             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                                                         <  1:3807783/150‑1 (MQ=255)
             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGATTATACATAACAGTCGTTTTTTTTATTTATAAAAAAATTCAAATAAATAAAGACAAAAGCCGTGTAAAGCACTAAGTTTAAATTTCGGGcccaacaag                                                                                                                                         >  2:4556361/1‑143 (MQ=255)
               gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATAAAGAAAAAAGAGATTTAAATAAATAATTTTAAATTTCGCGCCAATAATTgc                                                                                                                                       >  2:3747313/1‑150 (MQ=255)
                agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAATAAAAACAAAAGCGATGTTAATAACGAAGTTTTAATTTGGGGCCAAACATgggg                                                                                                                                      >  1:7338499/1‑148 (MQ=255)
                agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTTTAAAATAATTAAAAGAAATAAATACAATAGAGATTTAAAGAAATAAGTTTTAATTGAGGGCCCATCATGgcg                                                                                                                                      >  1:6895230/1‑150 (MQ=255)
                agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg                                                                                                                                      >  2:7438704/1‑150 (MQ=255)
                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                                                                     <  1:1634329/150‑1 (MQ=255)
                       ataatgcaTTTATTTGGATTTTATTTTTCGTTCGACCTTTGCTTATAAATAACAGTCTTTTTTTTATTTTATAAAATAATTAAAAGAAATCAAGACAATAGCGTTGTAAAGCACTAATTTTATATTGCGCGCCAATCATGGCGCGCACaa                                                                                                                               <  1:6229246/145‑1 (MQ=255)
                       agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                                                               <  2:1819666/150‑1 (MQ=255)
                          aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATTAAATAAATAAAAGAAAAAAAGAAAAAAGAGATTTAAATGAAAAATTTTAAATTTAGCCCCAAACAAGGCGCGCCCAAGca                                                                                                                            >  1:3304514/1‑149 (MQ=255)
                            gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtat                                                                                                                          <  2:3828499/150‑1 (MQ=255)
                             caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAATAAAAAAAGACAATAGCGATGTAAATCAAAAATTTTAAATTTCGCCCCAAAAATTGCGCGCACaaacaaaa                                                                                                                         >  2:1017672/1‑144 (MQ=255)
                                   tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGAAAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATAAATGCGCGCACAAGCTataataaaa                                                                                                                   >  2:8792843/1‑147 (MQ=255)
                                              tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATAGCGATTTAAAGCACTAAGTTTAAATTGAGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCGAAATtct                                                                                                        >  2:2291029/1‑150 (MQ=255)
                                                tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCATGACAATAGCGATGTAAAGCACTATGTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                                                      <  1:2534211/150‑1 (MQ=255)
                                                     ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAATCAAGACAAAAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctactcctc                                                                                                 >  2:5501422/1‑150 (MQ=255)
                                                        gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                                                              <  1:10026911/150‑1 (MQ=255)
                                                            tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAATAAAATAGCGATGTAAAGCAAAAAGTTTAAATTGCGCGCCCATAATTGCGCGCACAAGCTATAATACCAACCTAATTTATCCAACTCTTATAGt                                                                                          >  2:4659034/1‑150 (MQ=255)
                                                               gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCTCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                                                       <  2:592040/150‑1 (MQ=255)
                                                                cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                                                      <  2:3011908/150‑1 (MQ=255)
                                                                 ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                                                                     <  1:18365/150‑1 (MQ=255)
                                                                 ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                                                                     <  1:112981/150‑1 (MQ=255)
                                                                    taCATAACCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                                                                  <  1:7751754/150‑1 (MQ=255)
                                                                    taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                                                                  <  2:2215204/150‑1 (MQ=255)
                                                                    taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                                                                  <  2:6708459/150‑1 (MQ=255)
                                                                        tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                                                              <  2:4993876/150‑1 (MQ=255)
                                                                        tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                                                              <  1:841498/150‑1 (MQ=255)
                                                                           cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                                                           <  1:8216175/150‑1 (MQ=255)
                                                                            aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                                                                          <  2:8121422/150‑1 (MQ=255)
                                                                             gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                                                         <  2:2508420/150‑1 (MQ=255)
                                                                             gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                                                         <  1:860652/150‑1 (MQ=255)
                                                                                 tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                                                     >  1:2307268/1‑150 (MQ=255)
                                                                                 tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                                                     >  2:3157/1‑150 (MQ=255)
                                                                                   tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTATAGTGACTATATCTCCTGATAGAATTg                                                                   >  1:9453466/1‑150 (MQ=255)
                                                                                   tttttAATTAATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCGCTTAGAGTGTCTATATGTCCTGATAGAATTg                                                                   >  1:5505011/1‑150 (MQ=255)
                                                                                      ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                                                                <  2:4111294/150‑1 (MQ=255)
                                                                                      ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                                                                <  1:8829762/150‑1 (MQ=255)
                                                                                       tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                                                               <  1:1911146/150‑1 (MQ=255)
                                                                                        aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                                                              >  2:4000654/1‑150 (MQ=255)
                                                                                         aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAGGGAAAGCACTAAGTTTCAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                                                             <  2:5057221/150‑1 (MQ=255)
                                                                                            tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                          >  2:5769707/1‑150 (MQ=255)
                                                                                            tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                          >  1:9233934/1‑150 (MQ=255)
                                                                                            tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                          >  1:2088376/1‑150 (MQ=255)
                                                                                             ataAAATAATTCAAAGAAATCAAGACAATCGCGATGTCAAGCACTAAGTTGAAATTGCGCGCCAATCATGGCGCGCACGCGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATCTCTCCTGATAGACGGGCGGTATTgcg                                                         >  2:3420194/1‑148 (MQ=255)
                                                                                             ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                                                         >  2:4770231/1‑150 (MQ=255)
                                                                                                aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                      >  1:5161380/1‑150 (MQ=255)
                                                                                                 aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                                                     >  2:5398610/1‑150 (MQ=255)
                                                                                                  ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa                                                    <  2:629281/150‑1 (MQ=255)
                                                                                                   taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa                                                   <  1:4420317/150‑1 (MQ=255)
                                                                                                    atttaaaaGAAATCAAGACTATAGCGATGTAAAGCTCTTTTTTTAAATTGCGCGACAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                                                  <  1:3747313/145‑1 (MQ=255)
                                                                                                    aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                                                  <  2:6077747/150‑1 (MQ=255)
                                                                                                     atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt                                                 <  2:8090450/150‑1 (MQ=255)
                                                                                                      ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                                                >  1:8051777/1‑150 (MQ=255)
                                                                                                       tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                                               <  2:7654607/150‑1 (MQ=255)
                                                                                                        cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                                              >  2:2969700/1‑150 (MQ=255)
                                                                                                        cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                                              <  2:8463124/150‑1 (MQ=255)
                                                                                                        cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCCCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                                              <  2:9045772/150‑1 (MQ=255)
                                                                                                         aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGtcc                                             >  2:8323325/1‑150 (MQ=255)
                                                                                                         aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc                                             >  1:3906467/1‑150 (MQ=255)
                                                                                                            gAAATCAAGACAATAGCGATTTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGAGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTATAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa                                          <  1:4799633/150‑1 (MQ=255)
                                                                                                             aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt                                         <  1:6179695/150‑1 (MQ=255)
                                                                                                              aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt                                        >  2:4504525/1‑150 (MQ=255)
                                                                                                               aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc                                       <  2:3304514/150‑1 (MQ=255)
                                                                                                                tCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg                                      <  1:4895093/150‑1 (MQ=255)
                                                                                                                 caatacaataTCTATGTAAAGCACTTAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTTTTGACTATCAAATGCCCTGATTCgt                                     <  2:197897/150‑1 (MQ=255)
                                                                                                                 cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt                                     <  1:65463/150‑1 (MQ=255)
                                                                                                                  aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt                                    >  1:1392912/1‑150 (MQ=255)
                                                                                                                  aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt                                    <  1:2807658/150‑1 (MQ=255)
                                                                                                                  aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt                                    <  2:1525402/150‑1 (MQ=255)
                                                                                                                   aGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttg                                   <  1:992202/150‑1 (MQ=255)
                                                                                                                    gACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgt                                  >  2:285494/1‑150 (MQ=255)
                                                                                                                         tAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTAtt                             <  1:5772695/150‑1 (MQ=255)
                                                                                                                         tAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTAtt                             <  1:7923075/150‑1 (MQ=255)
                                                                                                                          aGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTg                            <  2:9871478/150‑1 (MQ=255)
                                                                                                                          aGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTg                            >  1:7664488/1‑150 (MQ=255)
                                                                                                                            cGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTa                          >  2:925256/1‑150 (MQ=255)
                                                                                                                             gATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTaa                         <  2:7113217/150‑1 (MQ=255)
                                                                                                                                   aaGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATAtctc                   <  1:9047114/150‑1 (MQ=255)
                                                                                                                                   aaGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATAtctc                   >  2:7290918/1‑150 (MQ=255)
                                                                                                                                   aaGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATAtctc                   >  1:3217873/1‑150 (MQ=255)
                                                                                                                                      cACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcctc                <  1:3456401/150‑1 (MQ=255)
                                                                                                                                       aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGATTGACTATATCTACTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct               <  2:2358143/150‑1 (MQ=255)
                                                                                                                                        cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCGCTGATTCGTTGTTTATTGTCATATCTCctctc              >  2:9600251/1‑149 (MQ=255)
                                                                                                                                         tAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTAt             >  2:1672212/1‑150 (MQ=255)
                                                                                                                                           aGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCt           <  2:1385406/150‑1 (MQ=255)
                                                                                                                                           aGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCt           <  1:1243266/150‑1 (MQ=255)
                                                                                                                                            gTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTg          >  1:6340306/1‑150 (MQ=255)
                                                                                                                                             tttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGc         >  2:3656666/1‑150 (MQ=255)
                                                                                                                                             tttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGc         >  2:2606091/1‑150 (MQ=255)
                                                                                                                                             tttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGc         >  2:2319188/1‑150 (MQ=255)
                                                                                                                                              ttAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCa        >  1:8330614/1‑150 (MQ=255)
                                                                                                                                              ttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCa        <  2:6504195/150‑1 (MQ=255)
                                                                                                                                                aaaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCGTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGa      <  1:1672225/150‑1 (MQ=255)
                                                                                                                                                aaaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGa      <  1:4673981/150‑1 (MQ=255)
                                                                                                                                                 aaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGAc     >  2:8628130/1‑150 (MQ=255)
                                                                                                                                                   ttGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGa   <  2:7519718/150‑1 (MQ=255)
                                                                                                                                                    tGCGCGCCAATCATGGCGCGCACAAGCTATAATACCTACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCTGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGAt  <  2:6929485/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGAT  >  NC_000913/2470323‑2470620

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: