Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,441 C→A 5.9% 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 5.9% 1167.9 / 12.2 373D437E (GAC→GAAgtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (143/191);  new base A (36/3);  total (179/194)
Fisher's exact test for biased strand distribution p-value = 8.50e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

TGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTT  >  NC_000913/2470294‑2470590
                                                                                                                                                   |                                                                                                                                                     
  gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTCAAATAAATAAAGTCAATa                                                                                                                                                   >  2:6386049/1‑150 (MQ=255)
  gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATCAAAGAAAACAAGaaaaaa                                                                                                                                                   >  2:1819907/1‑148 (MQ=255)
   aaTATCTCAATTTTATTTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAAGTAATCAAGAAAATaa                                                                                                                                                  >  1:6509038/1‑149 (MQ=255)
   aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAAATTATAAAAAAATTCAAAGAAAAAAAGACaaaac                                                                                                                                                  >  1:3004233/1‑147 (MQ=255)
      atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                               <  1:3163208/150‑1 (MQ=255)
       tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                              <  1:1745790/150‑1 (MQ=255)
        ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAATAAAAAATCGAtt                                                                                                                                             >  2:7893010/1‑149 (MQ=255)
        ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATAAACAAAAAAGCCATg                                                                                                                                             >  2:6507356/1‑150 (MQ=255)
        ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTAAAAAATAAATAAAATAAAACAATAAAAAAGCGATg                                                                                                                                             >  1:8391008/1‑150 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTAAAAGAAATAAAGACAAAATCGATGt                                                                                                                                            >  2:80240/1‑150 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGt                                                                                                                                            >  1:6947772/1‑150 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTAAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAAAAGCGATGt                                                                                                                                            >  2:2817711/1‑150 (MQ=255)
           aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                                                                                                          <  1:7791498/150‑1 (MQ=255)
           aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAAAATCGATGTaa                                                                                                                                          >  2:316251/1‑150 (MQ=255)
             ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                                        <  1:7581636/150‑1 (MQ=255)
               ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                                      <  1:4693329/150‑1 (MQ=255)
               ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCa                                                                                                                                      >  2:2664027/1‑150 (MQ=255)
                tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                                                                                                                     <  2:2374952/150‑1 (MQ=255)
                 aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                                                                                                                    <  2:5127131/150‑1 (MQ=255)
                 aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGCCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                                                                                                                    <  2:880049/150‑1 (MQ=255)
                         tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                                                                                                            <  2:7168074/150‑1 (MQ=255)
                         tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAAAAAATAAAAAATCGATGTAAAGAACTAAGTTTaa                                                                                                                            >  1:2007810/1‑150 (MQ=255)
                            ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAATAAAATAAAAAAAAAAAAAATCGATGTAAAGCACAAAAGTTaaaat                                                                                                                         >  2:2253312/1‑148 (MQ=255)
                            ttCGGGGAAGGAAATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGATTTTTTAATTTATAAAATAATTAAAAAAAATAAAAAACAAAACGATGTAAATCACAAACTTTAAAtt                                                                                                                         >  1:5047637/1‑150 (MQ=255)
                                    aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                                 <  2:3584475/150‑1 (MQ=255)
                                    aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                                 <  2:4260656/150‑1 (MQ=255)
                                          taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                           <  1:7070731/150‑1 (MQ=255)
                                            gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAAAATCGATGTAAAGCACAAAGTTTAAATTTCGCGCCAAACAAGgc                                                                                                         >  2:4020940/1‑150 (MQ=255)
                                               aTTTAAGCAAGCATTTATTTGGAATGTATTTTTTGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg                                                                                                      <  2:6615125/150‑1 (MQ=255)
                                                  tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac                                                                                                   <  2:5633674/150‑1 (MQ=255)
                                                    agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATTAAAAAAATAAAAAAAAAAAAAAAAAAAACGATTTAAAGCGAAAATTTTAAAATTAGCGCCCAGCCAGGCGCGCACaa                                                                                                 >  1:785914/1‑150 (MQ=255)
                                                      caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATCAAGAAAAAAGCGATGTAAAGCACTAATTTTAAATTTAGCGCCAATCATTGCGCGCACaacc                                                                                               >  2:1858346/1‑148 (MQ=255)
                                                       aagcaTGTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACAAAACAGTCGTTTTTTTTATTTAATAAAAAATTCAAATAAAAAGAGTCAAAAAAGAATTAAAGAACCAAATTTAGATTTAGCGCCCAAAATGCCGCGCCCaaaag                                                                                              >  1:6712139/1‑147 (MQ=255)
                                                            tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                                                                         >  2:4795338/1‑150 (MQ=255)
                                                            tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAAATAAAATAAATAAAAACAAAATCGATGTAAAACACAAAGTTTAAATTTAGCGCCAAGCAAGGCGCGCCCAACCGataaa                                                                                         >  1:5244309/1‑149 (MQ=255)
                                                                 ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCaa                                                                                    <  2:4125485/150‑1 (MQ=255)
                                                                      aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTaa                                                                               <  2:7447972/150‑1 (MQ=255)
                                                                       tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAATAAAGAAAAAATAGATTTAAATCACAAAATTTAAATTTAGCGCCAAACAATGCCCGCACAAGCAATAAAACCaaacaaaa                                                                              >  2:3681701/1‑144 (MQ=255)
                                                                         tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCAAAttt                                                                            >  2:711578/1‑150 (MQ=255)
                                                                         tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAATAAAGAAAAAAACGATGTAAATCACGAAGTTTAAATTTCGCGCCAATCAATGCGCGCACAAGCGATAATACCAACCAAAttt                                                                            >  2:2505673/1‑150 (MQ=255)
                                                                          aTTTTTAGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                                           <  2:4632560/150‑1 (MQ=255)
                                                                           tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtct                                                                          <  2:7719576/150‑1 (MQ=255)
                                                                            ttttCGTTCGACCTTTGCTTATACATAACCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                         <  2:3462/150‑1 (MQ=255)
                                                                                 gttctACCTTTTATTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCTATGTAAAGCAATATTTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                                    <  1:2235842/150‑1 (MQ=255)
                                                                                 gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                                    <  2:2421713/150‑1 (MQ=255)
                                                                                  ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTATTTTATAAAATAATTCAAATAAATCAAGACAATAGCTATGTAAAGCAATAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcctc                                                                   <  1:4108069/150‑1 (MQ=255)
                                                                                   tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATAATGGCGCGCACAAGCTATAATACCAACCCAAATTCTCCTCcact                                                                  >  2:8078014/1‑150 (MQ=255)
                                                                                    cgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAATAGCGATGGAAAGCACTAAGGTTAAATTGCGCGCCAATCAAGGCGCGCACAAGCTATAAAACCAACCAAAATTctcctcctcat                                                                 >  2:1202851/1‑148 (MQ=255)
                                                                                     gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                                <  1:3051389/150‑1 (MQ=255)
                                                                                           tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                                          <  2:2879661/150‑1 (MQ=255)
                                                                                            gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                         <  2:3026566/150‑1 (MQ=255)
                                                                                             cTTATACATAACAGTCTTTTTTTTATTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                        <  1:6559040/150‑1 (MQ=255)
                                                                                                 taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                                    <  1:7111708/150‑1 (MQ=255)
                                                                                                 taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATAGTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCGTAGAGTGACTatat                                                    >  2:4730485/1‑150 (MQ=255)
                                                                                                  aCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtc                                                   <  1:1437247/150‑1 (MQ=255)
                                                                                                       aCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTg                                              <  2:8307227/150‑1 (MQ=255)
                                                                                                          ttcttttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                           <  1:5708282/146‑1 (MQ=255)
                                                                                                            cGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACAAAGTTTAAATTGCGCGCCCCTCATGGCGCGCCCAAGCTATAATACCAACCGAATTTCTCCTCCTCTTAGAGTGACTATAACTCCGGATata                                         >  1:6280827/1‑148 (MQ=255)
                                                                                                             gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                                        >  2:6124473/1‑150 (MQ=255)
                                                                                                              tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  1:1130887/1‑150 (MQ=255)
                                                                                                              tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  2:2963891/1‑150 (MQ=255)
                                                                                                              tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                       >  1:6463117/1‑150 (MQ=255)
                                                                                                                 ttttattttaTAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAATCAATATGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTACTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                    <  1:2867084/150‑1 (MQ=255)
                                                                                                                 ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                    <  1:3433952/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                   <  2:5416188/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                   <  2:5578604/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                   <  2:4905783/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                   <  2:1888897/150‑1 (MQ=255)
                                                                                                                       tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                              >  2:4504594/1‑150 (MQ=255)
                                                                                                                       tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                              >  1:1584375/1‑150 (MQ=255)
                                                                                                                       tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                              >  2:7807256/1‑150 (MQ=255)
                                                                                                                        ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                             >  2:7273053/1‑150 (MQ=255)
                                                                                                                        ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                             >  2:1707820/1‑150 (MQ=255)
                                                                                                                        ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                             >  2:2719448/1‑150 (MQ=255)
                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                            >  1:8050842/1‑150 (MQ=255)
                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                           >  2:1105850/1‑150 (MQ=255)
                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                           <  2:3789617/150‑1 (MQ=255)
                                                                                                                            aaaaTAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                         <  1:2654560/150‑1 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                         <  1:958583/150‑1 (MQ=255)
                                                                                                                              aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                       <  1:5933064/150‑1 (MQ=255)
                                                                                                                              aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCCCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                       <  2:1458670/150‑1 (MQ=255)
                                                                                                                               ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa                      >  1:5460555/1‑150 (MQ=255)
                                                                                                                                taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa                     <  1:3402631/150‑1 (MQ=255)
                                                                                                                                 aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                    >  1:2481429/1‑150 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                  >  1:4299051/1‑150 (MQ=255)
                                                                                                                                    tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                 >  1:1895410/1‑150 (MQ=255)
                                                                                                                                    tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                 >  2:842483/1‑150 (MQ=255)
                                                                                                                                    tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                 <  1:7899717/150‑1 (MQ=255)
                                                                                                                                     cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCCTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                >  2:7598171/1‑150 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt              <  1:4377016/150‑1 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt              >  1:7300917/1‑150 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGCATTGACTATCAAATGCCCt              >  2:6980009/1‑150 (MQ=255)
                                                                                                                                        aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg             >  1:5791351/1‑150 (MQ=255)
                                                                                                                                         gAAATCAAGACTATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATTCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa            <  2:3671270/150‑1 (MQ=255)
                                                                                                                                            aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc         <  1:5337058/150‑1 (MQ=255)
                                                                                                                                            aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc         >  2:6364555/1‑150 (MQ=255)
                                                                                                                                               aaGACAATCGCGATGTATAGCTCTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt      <  2:3065735/150‑1 (MQ=255)
                                                                                                                                                aGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttg     >  2:5381258/1‑150 (MQ=255)
                                                                                                                                                   cAATCGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCCATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGttt  <  1:6049543/150‑1 (MQ=255)
                                                                                                                                                   cAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGttt  <  1:1053972/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTT  >  NC_000913/2470294‑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: