Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 20.5% Q433K (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,4270CA20.5% 589.3 / 50.1 236Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (72/108);  new base A (51/4);  total (124/112)
Fisher's exact test for biased strand distribution p-value = 7.56e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  NC_000913/2470279‑2470575
                                                                                                                                                    |                                                                                                                                                    
tGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCa                                                                                                                                                     <  2:2210769/150‑1 (MQ=255)
tGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCa                                                                                                                                                     <  2:1885258/150‑1 (MQ=255)
  tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaaa                                                                                                                                                   >  1:957910/1‑150 (MQ=255)
  tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaaa                                                                                                                                                   >  2:3093319/1‑150 (MQ=255)
  tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                   <  2:707622/150‑1 (MQ=255)
    gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGa                                                                                                                                                 >  1:2211012/1‑150 (MQ=255)
         tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                            <  2:3679607/150‑1 (MQ=255)
         tctTGTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTCAAAGAAATCa                                                                                                                                            <  2:1155081/150‑1 (MQ=255)
          ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAATCAAAAAAAACaa                                                                                                                                           >  2:111485/1‑150 (MQ=255)
          ctTTGTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGTAATTTATTTTTAGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                           <  2:1372084/150‑1 (MQ=255)
           tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATaaat                                                                                                                                          >  2:551364/1‑149 (MQ=255)
             tttGGAATATCTCAATGTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                                                                                                        <  1:3396036/150‑1 (MQ=255)
                ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACaaa                                                                                                                                     >  1:1531131/1‑149 (MQ=255)
                      tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                               <  2:1058763/150‑1 (MQ=255)
                        tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAATAAAAAAAGACAAAAGCGATGt                                                                                                                             >  1:3497829/1‑150 (MQ=255)
                        tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACACAATCGATGt                                                                                                                             >  2:2375162/1‑150 (MQ=255)
                        tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAACAAGACAAAAGCGATGt                                                                                                                             >  1:2322608/1‑150 (MQ=255)
                        tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAAAAGCGATGt                                                                                                                             >  1:498995/1‑150 (MQ=255)
                        tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAAAAAAAAAAAAAAAAGCGAAGt                                                                                                                             >  1:1255255/1‑150 (MQ=255)
                          aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                                                                                           <  1:3903370/150‑1 (MQ=255)
                            ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                         <  1:3146486/150‑1 (MQ=255)
                                aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                                                                                                     <  2:1935042/150‑1 (MQ=255)
                                 tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                                                                                                    <  1:725783/150‑1 (MQ=255)
                                  gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                   <  1:3517991/150‑1 (MQ=255)
                                   tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAACAAAAAAACACAAAAGCGAAGTAAAGAACAAAg                                                                                                                  >  2:1823914/1‑150 (MQ=255)
                                    gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGAGATGTAAAGCACTAAGt                                                                                                                 >  2:961314/1‑150 (MQ=255)
                                    gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCACGAAGt                                                                                                                 >  2:126537/1‑150 (MQ=255)
                                     tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAATACAAAAGCGATGTAAAGCACTAAGtt                                                                                                                >  1:566103/1‑150 (MQ=255)
                                      cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGAAAAAAGCGATGTAAAGCACTAAGttt                                                                                                               >  1:626788/1‑150 (MQ=255)
                                       aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                                                                                              <  2:3470911/150‑1 (MQ=255)
                                           ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                                                                                          <  2:1362276/150‑1 (MQ=255)
                                             cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                                                                                                        <  2:3740101/150‑1 (MQ=255)
                                               gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAATAAATACAAAAGCGATGTAAAGCACTAAGTTTAAAtttagc                                                                                                      >  1:418805/1‑146 (MQ=255)
                                               gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCGTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATCATTAAAAAAACGAGGTACAGCACTAAGTTTAAAATgcgc                                                                                                      >  2:2705369/1‑150 (MQ=255)
                                                 gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                    <  2:2445039/150‑1 (MQ=255)
                                                    ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAAAAAAAATCGATTTAAAACACAAAGTTTTATTTTGGCGCCaa                                                                                                 >  1:2656125/1‑150 (MQ=255)
                                                    ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAAAAAAAAATAAAAAAGAGATGGAAAGCACAAAGTTTTAATTTAGCGccca                                                                                                 >  1:1084248/1‑148 (MQ=255)
                                                     gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAAACAACAAAAAACCCATGTTAATCAAGAAGTTTTAATTTCGCGCCaaa                                                                                                >  1:110816/1‑149 (MQ=255)
                                                      atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTCATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                                                                                               >  2:2561049/1‑150 (MQ=255)
                                                       tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                                                                                              >  2:2284216/1‑150 (MQ=255)
                                                           gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAATTACAATAAAAAAAGAAAAAAACGATGTAAAGCACTAAGTTTTAATTTGGCGCCCATCCTggg                                                                                          >  1:3245494/1‑149 (MQ=255)
                                                                 aaagcaagcAATTATTTGGAATGTAATTTTCGTTCGACCTTTGCTTATACATAACAGTCCTTTTTTTAAATTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                    <  2:2035178/149‑2 (MQ=255)
                                                                   agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAATAAATAAAGACAAAATAGATGTTAAACACAAATTTTTAATTTCGAGCCAAACATGGCGCGCACaa                                                                                  >  1:2833669/1‑150 (MQ=255)
                                                                       agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                              <  2:3389067/150‑1 (MQ=255)
                                                                       agcaTTTATTTGGAATGTATTTTTCGTTCGACATTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                              <  2:3389244/150‑1 (MQ=255)
                                                                          atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                           <  2:491222/150‑1 (MQ=255)
                                                                           tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTGAAATAAATAAATTCAATAGGGATGTAAAGGAATTATTTTTATTTTCGGGCCAATCATTGCGCGCACAATCTttaat                                                                          >  2:1584164/1‑150 (MQ=255)
                                                                           tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAATTAAAAAAAAAAAAGAAAAAATAGAAGTAAACAACAAAGTTTAAATTTAGCGCCCAACATTGCGCGCACAAGcaaaaaa                                                                          >  1:153307/1‑144 (MQ=255)
                                                                              ttctGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATAGCGATGTAAATCACAAAGTTTAAATTGCGCGCCAATCATGGCGCGCACCAGCAATAAAAcc                                                                       >  2:3885975/4‑150 (MQ=255)
                                                                              ttctGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAAAAGCGATGTAAAGCACTAAGTTTAAAATGCGCGCCAATCATGGCGCGCACAAGCTataaaaac                                                                       >  2:3268397/4‑146 (MQ=255)
                                                                                 tGGAATGTATTTTTCGTTCTCCCTTGGCTAATACATAACAGCCTTTTTTTTAATTAAAAAAATAATTCAAAGAAATCAAGAAAATAGCTTTGTAAAGCACTAATTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAAc                                                                    <  1:3704326/150‑1 (MQ=255)
                                                                                       ttATTTTTCGTTCGACCTTTGATTATAAATAACAGACGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                              <  1:2416395/149‑1 (MQ=255)
                                                                                       gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                              <  2:895382/150‑1 (MQ=255)
                                                                                           ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                          <  2:3704055/150‑1 (MQ=255)
                                                                                            tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                         <  2:3767688/150‑1 (MQ=255)
                                                                                                gttcgACCTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                     <  1:2164489/150‑1 (MQ=255)
                                                                                                 ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATATCGATGTAAATCACTAAGTTTAAATTTCGCGCCAATCAATGCGCGCACAATCTATAAAAACAACCGAATTTatcctacac                                                    >  2:3146259/1‑146 (MQ=255)
                                                                                                   cgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                  <  2:283287/150‑1 (MQ=255)
                                                                                                    gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                 <  2:2118133/150‑1 (MQ=255)
                                                                                                    gACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTCCAAGAAAACAAGACAATAGCGAGGTCAAGCCCGAAGATTAAAATGCGCCCCACGCATGGCGCGCCCCAGCGGAAATACCAACCGAATTTCTCCACCCCTTa                                                 >  2:1975916/1‑150 (MQ=255)
                                                                                                      ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  2:789610/150‑1 (MQ=255)
                                                                                                        tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGt                                             <  1:3564384/150‑1 (MQ=255)
                                                                                                          ttcTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                           <  1:2125412/148‑1 (MQ=255)
                                                                                                          tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATCTCTCCTCCTCTTAGAGTGa                                           <  1:3283608/150‑1 (MQ=255)
                                                                                                          tGCTTATACATAAAAGTCTTTTTTTTAATTTATATAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCATTCTTGGCGCGCACTAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                           <  1:2529792/150‑1 (MQ=255)
                                                                                                           gctaaaaCATAACAGTCTTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGAAAATTGCTATTTAAAGCACTAAGTTTAAATTTCGCGCAAATCATGCCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                          <  1:1711822/144‑1 (MQ=255)
                                                                                                           gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGATTGAc                                          <  2:292346/150‑1 (MQ=255)
                                                                                                             ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCCCAAGCTATAATACCAACCAAATTTCTCCTCCTATTAGAGGGACta                                        >  2:2928672/1‑150 (MQ=255)
                                                                                                             ttATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAAGAAAAAAAGAAAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATAATTGCGCGCACAAGCTATAATAACAACCAAAATTATCCTCCGCTTGGAGTGACta                                        >  2:399143/1‑150 (MQ=255)
                                                                                                                  cATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct                                   <  2:1820766/150‑1 (MQ=255)
                                                                                                                   aTAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAAAAAAGAAAAAAGCGATGTAAATCACAAAGGTTAAAATTCGCGCCAAAAATGGCGCGCCCAAGCAAAAATACCAAACAAATATATACACCTCTTAGAGTGACtaaaaatc                                  >  2:1503497/1‑144 (MQ=255)
                                                                                                                   aTAACAGTCGTTTTTTTAAATTAAAAAATAAATCAAATAAATAAAGAAAAAAGCGATGTAAATCAATAATTTTAAATTTCGCCCCAATCAAGGCGCCCACCAGCAAAAAAACCAAACAAAATTCTCCACCTCATAGAGTGACAATAtctc                                  >  2:3427341/1‑150 (MQ=255)
                                                                                                                     aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                <  1:1501715/150‑1 (MQ=255)
                                                                                                                        aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                             <  2:3848055/150‑1 (MQ=255)
                                                                                                                             tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                        >  1:1780931/1‑150 (MQ=255)
                                                                                                                             tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                        >  1:2123506/1‑150 (MQ=255)
                                                                                                                             tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                        >  1:134016/1‑150 (MQ=255)
                                                                                                                               tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAGTCATGGCGCGCACAAGCTGTAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATGTCTCCTGATAGAATTg                      >  1:681500/1‑150 (MQ=255)
                                                                                                                               tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                      <  2:2133495/150‑1 (MQ=255)
                                                                                                                                ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                     >  1:2220538/1‑150 (MQ=255)
                                                                                                                                 tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTACTTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    <  2:1710693/150‑1 (MQ=255)
                                                                                                                                 tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    >  2:578368/1‑150 (MQ=255)
                                                                                                                                     aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCCCCTGATAGAATTGCGGTAt                >  2:389061/1‑150 (MQ=255)
                                                                                                                                     aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAAACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                >  2:1294651/1‑150 (MQ=255)
                                                                                                                                      tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt               >  2:1237055/1‑150 (MQ=255)
                                                                                                                                      tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt               <  2:3901851/150‑1 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  1:3086720/1‑150 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  1:2898583/1‑150 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  1:285654/1‑150 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  2:288529/150‑1 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  2:1054851/1‑150 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  1:2285738/1‑150 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  2:2828081/1‑150 (MQ=255)
                                                                                                                                           aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa          >  1:3138944/1‑150 (MQ=255)
                                                                                                                                           aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa          >  1:3817910/1‑150 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        <  2:3377071/150‑1 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        <  1:1187538/150‑1 (MQ=255)
                                                                                                                                                   tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  >  2:3610875/1‑150 (MQ=255)
                                                                                                                                                   tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  >  1:2709706/1‑150 (MQ=255)
                                                                                                                                                   tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  >  2:3358844/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  NC_000913/2470279‑2470575

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: