Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,436 C→A 13.0% Q436K (CAA→AAA)  gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4360CA13.0% 567.6 / 18.2 193Q436K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (71/90);  new base A (28/3);  total (100/93)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGT  >  NC_000913/2470296‑2470585
                                                                                                                                            |                                                                                                                                                     
  atatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                                                                                                            <  2:1746298/150‑1 (MQ=255)
    atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAAAAGCGa                                                                                                                                          >  2:4071395/1‑150 (MQ=255)
     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACGTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATATTTGAAATAAATAAATACAATatagat                                                                                                                                         >  2:1000574/1‑145 (MQ=255)
       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAAACAAAAACGAAGt                                                                                                                                       >  2:1198444/1‑150 (MQ=255)
       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAATCAAGAAAATAGCGATGt                                                                                                                                       >  2:3356531/1‑150 (MQ=255)
         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATAAAATAGCGATGtta                                                                                                                                     >  2:3946658/1‑148 (MQ=255)
         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAAAATCGATGTaa                                                                                                                                     >  2:1909849/1‑150 (MQ=255)
         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGAAAAAAGCGATTTaa                                                                                                                                     >  1:73078/1‑150 (MQ=255)
          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                                    <  1:741723/150‑1 (MQ=255)
           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                                   <  2:2595762/150‑1 (MQ=255)
                  tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGAGATGTAAATCACTAAg                                                                                                                            >  2:3694262/1‑150 (MQ=255)
                   gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                                                                                                           <  1:2793112/150‑1 (MQ=255)
                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAAAAATAAAAAAAATATCGAATTAAAGAAATAAttt                                                                                                                          >  2:2837483/1‑150 (MQ=255)
                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                                                                                                        <  2:500583/150‑1 (MQ=255)
                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAAAAAAAAAAACAAAAGAGATGTAAATGACAAATTTTa                                                                                                                        >  1:1060333/1‑150 (MQ=255)
                       tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAAAAAAAAAAAACAAAAGAGAAGTAAAGACCAAAGATTaa                                                                                                                       >  1:3942719/1‑150 (MQ=255)
                         tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTCAAAGAAAACAAGACAATATCGAATTTAATCAAAAAGTTTAAAt                                                                                                                     >  2:3138509/1‑150 (MQ=255)
                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAACAAAACAAGACAAAAGCGATGTAAAGAACTAAGATTAAAtt                                                                                                                    >  2:3503103/1‑150 (MQ=255)
                            cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                                                                                                                  <  2:1597131/150‑1 (MQ=255)
                             ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAATAAAAGAAAAAAAAACAAAAGCGATTTAAATCAATAAATTTTAATTtcg                                                                                                                 >  1:2162344/1‑150 (MQ=255)
                              gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGGATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                                                                                                                <  1:388155/150‑1 (MQ=255)
                              gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                                                                                                                <  2:1046688/150‑1 (MQ=255)
                               ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAAATACAAGAAAACAACACAAAAGCGAAGTAAAGCACAAACGTTAAATTgcccg                                                                                                               >  2:389809/1‑150 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                              >  2:3098031/1‑150 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                              <  2:1539263/150‑1 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTtcgcgc                                                                                                              >  2:1276990/1‑150 (MQ=255)
                                    gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTTAGCGCCaaa                                                                                                          >  2:917518/1‑149 (MQ=255)
                                      tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAATAAAAAAAGACAAAAGCGATGTAAAGCACAAAGTTTAAAATGCGCGCCACTCa                                                                                                        >  1:3698835/1‑150 (MQ=255)
                                       ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                                                                       <  1:3766376/150‑1 (MQ=255)
                                        taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                      <  1:1911594/150‑1 (MQ=255)
                                        taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                      <  1:1862219/150‑1 (MQ=255)
                                          tataTTTAAGCAAGAATTTATTTGGAATTTATTTTTCTTTCGACCTTTGCTTATACATAACAGTCGATTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                    <  2:420900/147‑1 (MQ=255)
                                          gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                    <  1:2754626/150‑1 (MQ=255)
                                          gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                    <  2:1045010/150‑1 (MQ=255)
                                              tttACGCAAGATTTTATTTGTAATGTATTTTTCGTTCGACGTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATTTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGgcgcgc                                                                                                <  2:1881228/150‑1 (MQ=255)
                                                   gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGAAAATAGAGATTTAAATCACGAAGTTTTAATTTAGAGCCAATAATGGCGCGCCCAAg                                                                                           >  1:2559128/1‑150 (MQ=255)
                                                     aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAAAAGCGATGTTAAGCACTAATTTTAAATTGCGCGCCCATCCTGGGGCGCACAAGCt                                                                                         >  2:3871390/1‑150 (MQ=255)
                                                      agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                        <  2:55296/150‑1 (MQ=255)
                                                      agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                        <  2:2147921/150‑1 (MQ=255)
                                                          tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                                                                    >  2:2380768/1‑150 (MQ=255)
                                                            tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                                                                  <  2:2742188/150‑1 (MQ=255)
                                                             atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                                 <  2:1338718/150‑1 (MQ=255)
                                                                  gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTGAGCGCCAATCATTGCGCGCACAAGCTATAATACCAACCt                                                                            >  2:244924/1‑150 (MQ=255)
                                                                       tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                                       >  2:2256712/1‑150 (MQ=255)
                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                    <  2:1116598/150‑1 (MQ=255)
                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                    <  2:2112396/150‑1 (MQ=255)
                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                    >  2:2666556/1‑150 (MQ=255)
                                                                           tttctttcTACTTTTTATTATAAATAACTTTCTTTTTTTTAATTTATAAAATAATTAAAACAAATCAATAAAATAGCGATTTAAAGCACTATTTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                   <  1:1144259/150‑1 (MQ=255)
                                                                            ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTgtcct                                                                  >  2:3002653/1‑150 (MQ=255)
                                                                              cgttcgACCTTTGCTTATACATAACTGTGTTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCGATTTAAAGCACTAATTTTAATTTGCGCGCAAATCAGGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                                <  1:1124354/150‑1 (MQ=255)
                                                                              cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACATTAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                                <  1:435931/150‑1 (MQ=255)
                                                                                     ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                         <  1:4035298/150‑1 (MQ=255)
                                                                                     ccTTTGCTTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                         <  1:1698216/150‑1 (MQ=255)
                                                                                     ccTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATGGCGCGCCAATCATGGCGCGCACAAGCTATAATGCCAACCTAATTTCTCCTCCTCTTaga                                                         <  1:3908443/150‑1 (MQ=255)
                                                                                     acTTTGCTTATACATAACATTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                         <  1:3965632/149‑1 (MQ=255)
                                                                                            ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                                  <  2:3457426/150‑1 (MQ=255)
                                                                                            ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                                  <  1:4109224/150‑1 (MQ=255)
                                                                                            ttATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                                  <  1:2801714/150‑1 (MQ=255)
                                                                                             tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATCGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                                 <  1:2891970/150‑1 (MQ=255)
                                                                                             tataCATAACACTCGATTTATAAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGATTAAATTGCGCGCCAATCATGGCGCTCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                                 <  2:3966058/150‑1 (MQ=255)
                                                                                                aCATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtc                                              <  1:1287458/150‑1 (MQ=255)
                                                                                                 cATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAACAAGACAAAACCCAAGTAAAGCACAAACTTTAAATTGCGCGCCCAACAAGGCCCGCACAAGCAATAAAACCAACCGAATTTCTCCTCCACGAAGAGTGACGATAtca                                             >  2:2643962/1‑149 (MQ=255)
                                                                                                 cATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAGAAGACAATATAGATGTAAATCACTAAGTTTAAATTTCGCGCCCATAATTGCGCGCACCAGCTATAATACAAACCAAAATTGTCCTCCGGTTAGGGTGACAATAtct                                             >  2:3623804/1‑150 (MQ=255)
                                                                                                    aaCAGTCGTTTTTTTAAATTATAAAAAAATTAAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAACATATAATACCAACCAAAATTCTCCTCCACATAGAGTGACTATATCACCt                                          >  2:306362/1‑150 (MQ=255)
                                                                                                        gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTAAATCTCCTGATa                                      >  2:3024151/1‑150 (MQ=255)
                                                                                                           ttttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                                   <  1:2175866/149‑1 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                  >  1:3638551/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                  >  2:2937071/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                  >  2:684954/1‑150 (MQ=255)
                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAGTGCGCGCCCCTCAGGGCGCGCACAAGCAATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGact                                  >  1:3328834/1‑148 (MQ=255)
                                                                                                             ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                                 <  1:4096477/150‑1 (MQ=255)
                                                                                                               ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                               >  1:2306700/1‑150 (MQ=255)
                                                                                                               ttttAAATTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                               <  1:2600870/150‑1 (MQ=255)
                                                                                                                tttAAATTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                              >  1:3357847/1‑150 (MQ=255)
                                                                                                                    aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                          <  1:2837483/150‑1 (MQ=255)
                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                         <  1:917518/150‑1 (MQ=255)
                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAGTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGTGTGACTATATCTCCTGATAGAATTGCGGTAtt                         <  1:1343883/150‑1 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                       <  2:2558516/150‑1 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                       >  1:1931166/1‑150 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                       <  1:1103984/150‑1 (MQ=255)
                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACAAAGTTTAAATTGCGCGCCAATCATTGCGCGCACAAGCTATAATACCAACCTAAATGCTCCTCCACATAGAGAGACTATATCGCCTGATAGAATTGCGGTATAGa                       >  2:938620/1‑150 (MQ=255)
                                                                                                                         taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                     <  2:1418787/150‑1 (MQ=255)
                                                                                                                         taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                     >  1:1378883/1‑150 (MQ=255)
                                                                                                                         taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCCGATAGAATTGCGGTATTGACt                     <  1:1021342/150‑1 (MQ=255)
                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                  >  2:3160729/1‑150 (MQ=255)
                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                  >  1:2383235/1‑150 (MQ=255)
                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                  >  2:16680/1‑150 (MQ=255)
                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                  >  1:1428651/1‑150 (MQ=255)
                                                                                                                               aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa               >  2:1733141/1‑150 (MQ=255)
                                                                                                                                atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGACTGGCGGGATTGACTATCCAAt              >  1:321532/1‑150 (MQ=255)
                                                                                                                                atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt              >  1:3985765/1‑150 (MQ=255)
                                                                                                                                atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt              >  1:2579471/1‑150 (MQ=255)
                                                                                                                                    aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc          >  1:1158489/1‑150 (MQ=255)
                                                                                                                                    aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc          >  2:1481966/1‑150 (MQ=255)
                                                                                                                                    aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc          >  1:14740/1‑150 (MQ=255)
                                                                                                                                    aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc          >  2:3459995/1‑150 (MQ=255)
                                                                                                                                    aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAGTGACTATCAAATGccc          >  1:1572662/1‑150 (MQ=255)
                                                                                                                                        aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt      <  2:361048/150‑1 (MQ=255)
                                                                                                                                          aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc    <  2:3942671/150‑1 (MQ=255)
                                                                                                                                           tCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg   >  2:2964789/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGT  >  NC_000913/2470296‑2470585

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: