Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 18.7% 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,4270CA18.7% 731.5 / 49.2 296Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (71/158);  new base A (59/6);  total (132/164)
Fisher's exact test for biased strand distribution p-value = 9.18e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  NC_000913/2470280‑2470575
                                                                                                                                                   |                                                                                                                                                    
gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaa                                                                                                                                                    >  2:682466/1‑150 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                   <  1:4479700/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTTTTAAATTATTaaaa                                                                                                                                                   >  1:6195124/1‑150 (MQ=255)
   gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCaaata                                                                                                                                                 >  2:1948913/1‑148 (MQ=255)
    taTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaa                                                                                                                                                >  2:4762748/1‑150 (MQ=255)
    taTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaa                                                                                                                                                <  2:1736593/150‑1 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                               <  1:937722/150‑1 (MQ=255)
     aTCTCATTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                               <  1:3083463/150‑1 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                                                                                              >  1:3788556/1‑150 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAt                                                                                                                                              >  2:4390506/1‑150 (MQ=255)
       ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAATaaaac                                                                                                                                             >  1:2256628/1‑148 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                            <  2:3198749/150‑1 (MQ=255)
         ctTTTTGGAATATCTCATTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGTATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                           <  1:4245731/150‑1 (MQ=255)
         ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                           <  1:5971515/150‑1 (MQ=255)
         ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                           <  1:769871/150‑1 (MQ=255)
         ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                           <  2:7007901/150‑1 (MQ=255)
              tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTAAAATAAATCAAGACaa                                                                                                                                      >  1:2003484/1‑150 (MQ=255)
                gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAAGAAaacaaaacaaaa                                                                                                                                    >  2:3448650/1‑148 (MQ=255)
                  atatatCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAATAAAAAAATACAAAAGc                                                                                                                                  >  1:4331245/1‑150 (MQ=255)
                  atatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                                                                                                  <  1:798614/150‑1 (MQ=255)
                     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGACTGTATTTTTCGTTCGACCATTGCTTATACATAACAGTCGTTTTTTTAATTTATAACATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                               <  1:3596924/150‑1 (MQ=255)
                     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                               <  2:4165329/150‑1 (MQ=255)
                      ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATg                                                                                                                              >  2:4583721/1‑150 (MQ=255)
                      ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATAGCGATg                                                                                                                              >  1:983543/1‑150 (MQ=255)
                      ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGAAAAAAGAGAtt                                                                                                                              >  1:6860943/1‑149 (MQ=255)
                        cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                            <  1:5363687/150‑1 (MQ=255)
                           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAATTAAAAGAAATAAAGACAATAGCGATTTAAAg                                                                                                                         >  2:6318216/1‑150 (MQ=255)
                           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                         <  2:1217694/150‑1 (MQ=255)
                           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                         <  2:2049221/150‑1 (MQ=255)
                           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAg                                                                                                                         >  2:479202/1‑150 (MQ=255)
                           ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAAGAAAACAAAACAATAGCGAAGTaaaa                                                                                                                         >  2:6430170/1‑149 (MQ=255)
                                tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                                                                                                    <  2:4241958/150‑1 (MQ=255)
                                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                   <  1:6401335/150‑1 (MQ=255)
                                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTaa                                                                                                                   >  2:2376401/1‑150 (MQ=255)
                                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                   <  1:2339899/150‑1 (MQ=255)
                                   gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                                                                                                 <  2:1553197/150‑1 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAATACAATAGCGATGTAAAGCACTAAGtt                                                                                                                >  2:1338813/1‑150 (MQ=255)
                                       tttttCGGTGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACATTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                                                                                             <  2:7224244/150‑1 (MQ=255)
                                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                            <  2:1008079/150‑1 (MQ=255)
                                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATACAAAGAAAAAAAGACAAAAGCGATGTAAAGAACAAAGTTTaaa                                                                                                            >  1:4058795/1‑150 (MQ=255)
                                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                                                                                          <  1:1714694/150‑1 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                       <  2:2555848/150‑1 (MQ=255)
                                              gggAAGGATATAGAGTTTTAAGCAAGGATTTATTTGGTATGTTTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATATTTCAAAGAAATCATGACTATAGCGATGTAAAGCACTAAGTTTTAATTgcgc                                                                                                      <  2:5771656/150‑1 (MQ=255)
                                              gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAAAAAAAAAAATATCGATGTAAATCACAAAGTTTAAAtttagc                                                                                                      >  1:1104321/1‑146 (MQ=255)
                                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTTAAATAAATAAAAGAAATAAAGAACAAATCGATTTAAATCACTAATGTTAAATTtagcgc                                                                                                    >  1:5636571/1‑150 (MQ=255)
                                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                    <  2:1112874/150‑1 (MQ=255)
                                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                    <  1:6572394/150‑1 (MQ=255)
                                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                    <  1:6692994/150‑1 (MQ=255)
                                                 aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAATAAAAACAAAAGAGATGTAAAGCACAAAGTTTAAATTTAGCGcc                                                                                                   >  1:443050/1‑150 (MQ=255)
                                                 aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAATAAATAAAGAAAAAATCGATGTAAATCAATAAGTTTAAATTGCGCGcc                                                                                                   >  2:5283087/1‑150 (MQ=255)
                                                       ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                                                             <  1:4207828/150‑1 (MQ=255)
                                                          gagaTTTAAGCAAGCTTTTATTTGGAATTTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                          <  2:5730443/150‑1 (MQ=255)
                                                             aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTGAAAGAAATAAATAAAATATGGATGTAAATCAATAAGTTTAAATTTGGCGCCAATCATTgcgcg                                                                                       >  2:5058642/1‑150 (MQ=255)
                                                                    caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAAATAAAAGAAAACAATACAAAAACGATGTAAATCACAAACTATAAAATTAGCCCCAAAAATGTCGCGCCCaacc                                                                                >  2:1860752/1‑148 (MQ=255)
                                                                       gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAAACAAGACAAAATCGGTGTTAAGCACTAAGGTTAAATTGCGCGCCCATCCTGGCGCGCACaaccgat                                                                             >  2:6924136/1‑145 (MQ=255)
                                                                        caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtata                                                                            >  2:5870012/1‑150 (MQ=255)
                                                                         atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                           <  2:4974684/150‑1 (MQ=255)
                                                                          tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGTAATCAAGAGAATATCGATGTAAAGCACTAATTTTAAATTGCGCGCCAAACAATGCGCGCACAAGCTataat                                                                          >  2:4490031/1‑150 (MQ=255)
                                                                             atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTTTTTATAAAATAAATAAAAAAAAAAAAGAAAAAATAGATGTAAAACACAAAGTTTAAATTTAGCGCCCAACAAGGGGCGCACAACCTaaaataac                                                                       >  1:3518158/1‑148 (MQ=255)
                                                                               ttGTAATTTTTTTTTGTTTCGACCTTTGCTTATACATAACAGTGTTTTTTTTTATTTATAAAATAATTCAAATAAATCAAGACATTCGCGATGTAAAGCACTATGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCaa                                                                     <  1:1034123/150‑1 (MQ=255)
                                                                                 ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTAAAAAATAATTCAAAGAAATCAATACAATACCGATGTAAAGCAATAAGTTTTAATTGCGCGCCACAAATGGCGCGCACAAGCTATAATACCAAcc                                                                   >  2:2102563/1‑150 (MQ=255)
                                                                                   attttATTTTTCTTTCGACCTTTGATTATAAATAACTGTCTTTTTTTTAATTTATAAAATAATTAAAATAAATAAATAAAATAGCGTTGTAAAGAAATATTTTTAAATTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCAACCTa                                                                 <  1:584874/146‑1 (MQ=255)
                                                                                       tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                             <  2:2235926/150‑1 (MQ=255)
                                                                                        aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCACCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                            <  1:2719665/150‑1 (MQ=255)
                                                                                        aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTATAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCACCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                            <  1:4872136/150‑1 (MQ=255)
                                                                                          ttttCTTTCGACCTTTGATTATAAATAACTGTTTTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCTATGTAAAGCAATATTTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                          <  1:70888/150‑1 (MQ=255)
                                                                                            ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAACAAACAAATTtctaca                                                        >  2:3572714/1‑147 (MQ=255)
                                                                                              cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                      <  2:3814770/150‑1 (MQ=255)
                                                                                               gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                     <  2:5688849/150‑1 (MQ=255)
                                                                                                     ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  1:6086073/150‑1 (MQ=255)
                                                                                                     ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  2:3890811/150‑1 (MQ=255)
                                                                                                        ttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAAACCAACCTAATTTCTCCTACTCTTAGAGTg                                            >  2:3024686/1‑150 (MQ=255)
                                                                                                          tcTTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                          <  1:4919138/149‑1 (MQ=255)
                                                                                                           cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAACTAATTTATACAACTCTTATAGTGACt                                         >  2:2912622/1‑150 (MQ=255)
                                                                                                           cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                         <  2:5534532/150‑1 (MQ=255)
                                                                                                            ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACAAAGTTTAAAATGCGCGCCCATCAAGGCGCGCACAAGCTATAATAACAACCAAAATTATCCTCCGCTTATCGTGACta                                        >  2:929892/1‑150 (MQ=255)
                                                                                                            ttATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATAAAGAAAAAAACGATTTAAATCACTAAGTTTAAATTTAGCGCCAAACAATGCGCGCACAAACAATAAAAACAAACAAAATTATACACCAATTAGAGTGaaaa                                        >  2:2631000/1‑147 (MQ=255)
                                                                                                              ataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACAAAGTTTAAATTGCGCGCCAAACAAGGCGCGCACAAGCTATAATACCAAACTAAATTCTCCTACACTTAGAGTGACgata                                      >  2:2412531/1‑150 (MQ=255)
                                                                                                               taCATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                     <  1:2870978/150‑1 (MQ=255)
                                                                                                                   tAACAGTCGTTTTTTTAAATTATAAAATAAATCAAAGAAATAAATAAAAAAGCGATGTAAAGCCCAAAGTTTAAATTGCGCGCCAATCATTGCGCGCACAACCAATAAAACCAACCTAATTTATACTCCGCTTAGAGTGACGAAATATcc                                 >  2:2877410/1‑150 (MQ=255)
                                                                                                                       aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAATCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                             >  2:546527/1‑150 (MQ=255)
                                                                                                                        gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                            <  2:7043912/150‑1 (MQ=255)
                                                                                                                           gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                         >  1:3684010/1‑150 (MQ=255)
                                                                                                                           gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                         >  1:5627686/1‑150 (MQ=255)
                                                                                                                            tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                        >  1:3142574/1‑150 (MQ=255)
                                                                                                                            tttttttAATTTATAAAATAAATCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGACAATACTGTCGCGCACAACCTATAATACCAACCTAATATCTCCTACTCTTAGAGTGACTATATCTACTGATAGaaa                        >  1:159246/1‑149 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    <  1:238322/150‑1 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    <  1:1396773/150‑1 (MQ=255)
                                                                                                                                   aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                 >  1:5185049/1‑150 (MQ=255)
                                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  2:3002584/150‑1 (MQ=255)
                                                                                                                                       tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:1813878/150‑1 (MQ=255)
                                                                                                                                       aataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:6084658/149‑1 (MQ=255)
                                                                                                                                         taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt           >  2:1168512/1‑150 (MQ=255)
                                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        >  2:595546/1‑150 (MQ=255)
                                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        >  2:5987196/1‑150 (MQ=255)
                                                                                                                                              taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa      <  1:4992749/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg   >  2:1203673/1‑150 (MQ=255)
                                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  >  1:6607163/1‑150 (MQ=255)
                                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  >  1:338523/1‑150 (MQ=255)
                                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  >  1:2437709/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  NC_000913/2470280‑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: