Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,480 A→T 6.0% intergenic (+18/‑335) 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,4800AT 6.0% 405.8 / 4.5 138intergenic (+18/‑335)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 A (51/77);  new base T (5/5);  total (56/82)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

ATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAA  >  NC_000913/2470333‑2470627
                                                                                                                                                   |                                                                                                                                                   
atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATACTTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                                                                                                                                                   <  2:1138492/150‑1 (MQ=255)
 tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATTTAAAGCACTAAGTTTAAATTGCTCGCCAATCa                                                                                                                                                  >  2:901153/1‑150 (MQ=255)
 tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                                                                                                                                                  <  1:1136798/150‑1 (MQ=255)
   taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAAAAATCAAGACAATAGCGATGTAAATCACTAATTTTAAATTGCGCGCCaaaaatg                                                                                                                                                >  2:357025/1‑145 (MQ=255)
   taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATTCATg                                                                                                                                                >  2:940080/1‑150 (MQ=255)
    agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                                                                               <  1:412437/150‑1 (MQ=255)
     gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCCATCATGgc                                                                                                                                              >  2:764425/1‑150 (MQ=255)
      agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg                                                                                                                                             <  1:346746/150‑1 (MQ=255)
      agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg                                                                                                                                             <  1:984150/150‑1 (MQ=255)
       gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCGATTTAAAGCAATAAGTTTAAATTTCGCGCCAATAATTgcgc                                                                                                                                            >  1:1642123/1‑150 (MQ=255)
         tttAAGCAATCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc                                                                                                                                          <  1:662573/150‑1 (MQ=255)
          ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGca                                                                                                                                         <  1:1700249/150‑1 (MQ=255)
           tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                                                                        <  1:476766/150‑2 (MQ=255)
            aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca                                                                                                                                       >  2:1246108/1‑150 (MQ=255)
               caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                                                                    <  2:350993/150‑1 (MQ=255)
                aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt                                                                                                                                   <  1:513658/150‑1 (MQ=255)
                 agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCACAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                                                                  <  2:152866/150‑1 (MQ=255)
                 agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                                                                  <  1:498144/150‑1 (MQ=255)
                 agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                                                                  <  1:491293/150‑1 (MQ=255)
                 agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                                                                  <  2:1179800/150‑1 (MQ=255)
                  ccaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtat                                                                                                                                 <  2:517212/149‑1 (MQ=255)
                       tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAGTTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCCCCAATCATGGCGCGCACAAGCTATAATAc                                                                                                                            <  1:597170/150‑1 (MQ=255)
                       tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGGCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                                                                                                            <  1:885726/150‑1 (MQ=255)
                        atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                                                                           >  2:1796671/1‑150 (MQ=255)
                           tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                                                                        <  2:586324/150‑1 (MQ=255)
                            ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                                                                       >  2:691936/1‑150 (MQ=255)
                            ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAca                                                                                                                       >  1:1089749/1‑149 (MQ=255)
                             gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCt                                                                                                                      <  2:256305/150‑1 (MQ=255)
                              aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTa                                                                                                                     <  1:15633/150‑1 (MQ=255)
                               aTGTATTTTTCGTTCTACCTTTGCTAAAACATAACAGTCGTTTTTTTAATTTATAAAATATTTCAAACAAATCAAGACAATAGCGATGTAAAGCACTAAGATTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTaa                                                                                                                    <  1:924474/150‑1 (MQ=255)
                                tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACATGCAATAATACCAACCTAAt                                                                                                                   >  2:53037/1‑150 (MQ=255)
                                  tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                                                                                 <  1:1005366/150‑1 (MQ=255)
                                         cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                                                                          <  2:27695/150‑1 (MQ=255)
                                             cgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                                                                      <  1:1310898/150‑1 (MQ=255)
                                              gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAAAAAATAACGATGTAAAGCACTAATTTTAAATTGCGCGCCAAAAATTGCGCGCACAAGCAATAATACCAACCTAATTTATACTACTCATa                                                                                                     >  1:1350637/1‑150 (MQ=255)
                                                ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAACTTTAAATTGCGAGCAAATAATGGCGCCCACAACCTATAATACCACCCTAATTTCTCCTCCTCTTaca                                                                                                   >  1:779194/1‑148 (MQ=255)
                                                  tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGt                                                                                                 >  2:251608/1‑150 (MQ=255)
                                                   ttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg                                                                                                <  2:404987/150‑1 (MQ=255)
                                                        tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAATCTATAATACCAACCTAATTTCTCATCCTCTTAGAGTGACtat                                                                                           >  2:1295658/1‑150 (MQ=255)
                                                         ataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtata                                                                                          <  1:438563/150‑1 (MQ=255)
                                                           aCATAACAGTAGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtc                                                                                        <  1:1269/150‑1 (MQ=255)
                                                              tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                                                                     >  2:340267/1‑150 (MQ=255)
                                                               aaCAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAACAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAAAACCAACCTAATTTATCCTCCTCTTAGAGTGACTATATGTCCt                                                                                    >  2:1434697/1‑150 (MQ=255)
                                                                  aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                                                                                 <  1:455260/150‑1 (MQ=255)
                                                                     cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                                                                              <  2:1064475/150‑1 (MQ=255)
                                                                     cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCACTAAGAGTGACTATATCTCCTCATAGa                                                                              >  1:1199817/1‑150 (MQ=255)
                                                                         tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                                                                          <  2:273239/150‑1 (MQ=255)
                                                                          ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                                                         <  1:417989/150‑1 (MQ=255)
                                                                          ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCATGATAGAATTGc                                                                         >  2:657849/1‑150 (MQ=255)
                                                                          ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATAATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                                                         <  2:894709/150‑1 (MQ=255)
                                                                          ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGAGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                                                         <  1:263943/150‑1 (MQ=255)
                                                                            ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                                                                       <  2:1078557/150‑1 (MQ=255)
                                                                            ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                                                                       <  1:938728/150‑1 (MQ=255)
                                                                             tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                                                                      <  1:1129511/150‑1 (MQ=255)
                                                                             tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                                                                      >  2:594697/1‑150 (MQ=255)
                                                                                tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCATCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                                                                   >  2:406348/1‑150 (MQ=255)
                                                                                  tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                 <  1:1791058/150‑1 (MQ=255)
                                                                                  tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                 >  1:308033/1‑150 (MQ=255)
                                                                                  tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                 >  1:1789774/1‑150 (MQ=255)
                                                                                          aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                                                         <  2:1767899/150‑1 (MQ=255)
                                                                                                aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                   <  1:172318/150‑1 (MQ=255)
                                                                                                aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                   <  1:858267/150‑1 (MQ=255)
                                                                                                aaGAAATCAAGACAATAGCGATGAAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                   >  1:439666/1‑150 (MQ=255)
                                                                                                 aGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGAGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg                                                  <  1:220750/150‑1 (MQ=255)
                                                                                                     aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc                                              <  1:1308026/150‑1 (MQ=255)
                                                                                                       cAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGAGCGCCAATCATGTCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt                                            <  1:357025/150‑1 (MQ=255)
                                                                                                       cAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCATTCATGGCGCGCACAAGCTATAATACCACCCTAATTTCTCATCCTCTAAGAGTGACTATATCTGCTGATAGAATTGCGGTATAGAGTATCACATGCACTGATTctt                                            >  2:738183/1‑148 (MQ=255)
                                                                                                        aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCACACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt                                           <  2:881244/150‑1 (MQ=255)
                                                                                                         aGACATTAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAATCATGGAGCGAAGAAGCTATAATAACAACATAATTTCTCCTCCTCTAAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCtttg                                          <  1:267292/150‑1 (MQ=255)
                                                                                                           aCAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgtt                                        >  1:1333272/1‑150 (MQ=255)
                                                                                                           aCAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgtt                                        >  1:687532/1‑150 (MQ=255)
                                                                                                                cgctaTGTAAAGCATTAAGTTTAAATTGGGCGCCATTCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGTATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTg                                   <  1:851897/146‑1 (MQ=255)
                                                                                                                aGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTg                                   <  1:1806735/150‑1 (MQ=255)
                                                                                                                aGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTg                                   >  1:618629/1‑150 (MQ=255)
                                                                                                                 gCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGt                                  <  2:1636493/150‑1 (MQ=255)
                                                                                                                 gCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGt                                  <  2:643933/150‑1 (MQ=255)
                                                                                                                 gCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGt                                  >  1:1702216/1‑150 (MQ=255)
                                                                                                                   taTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTaa                                <  1:833405/149‑1 (MQ=255)
                                                                                                                   gATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTaa                                >  1:1777932/1‑150 (MQ=255)
                                                                                                                   gATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGATATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTaa                                <  2:939119/150‑1 (MQ=255)
                                                                                                                     tGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAata                              <  1:1192482/150‑1 (MQ=255)
                                                                                                                      gTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAatat                             >  2:1484739/1‑150 (MQ=255)
                                                                                                                        aaaGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATAtct                           >  1:1020566/1‑150 (MQ=255)
                                                                                                                          aGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcc                         <  1:1821664/150‑1 (MQ=255)
                                                                                                                             aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct                      >  1:1714122/1‑150 (MQ=255)
                                                                                                                             aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct                      >  2:1170017/1‑150 (MQ=255)
                                                                                                                              cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTa                     >  2:300767/1‑150 (MQ=255)
                                                                                                                              cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTa                     >  1:1082622/1‑150 (MQ=255)
                                                                                                                                 aGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCt                  <  2:1734947/150‑1 (MQ=255)
                                                                                                                                  gTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTg                 >  1:1004454/1‑150 (MQ=255)
                                                                                                                                    ttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCa               <  2:1221659/150‑1 (MQ=255)
                                                                                                                                    ttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCa               <  1:1051136/150‑1 (MQ=255)
                                                                                                                                    ttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCa               <  2:1664133/150‑1 (MQ=255)
                                                                                                                                       ttgtGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTTCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATGTCCTCTATCTGCAGAc            >  1:531098/4‑150 (MQ=255)
                                                                                                                                       aaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGAGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGAc            <  1:901153/150‑1 (MQ=255)
                                                                                                                                        tgtGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTTATATCTCCTCTATCTGCAGACg           >  1:1520842/3‑150 (MQ=255)
                                                                                                                                          tGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGAt         >  1:1396759/1‑150 (MQ=255)
                                                                                                                                            cgcgcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGAGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCATTATTCTTTGTTTATTGTAATATCTCCTATATCTGCAGACGATaa       >  2:460657/1‑150 (MQ=255)
                                                                                                                                              ngcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACt     >  2:1304949/2‑150 (MQ=255)
                                                                                                                                                 cAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAGAGATTTGGGGTATTGACTATCATATGCCCTGATTCGGTGTTTATTGTAATATATCCTCTATCTGCTGACGATATCTGaa  >  2:1723140/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAA  >  NC_000913/2470333‑2470627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: