Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,441 C→A 9.9% D437E (GAC→GAA gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4410CA 9.9% 661.9 / 9.9 215D437E (GAC→GAAgtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (87/96);  new base A (29/2);  total (117/98)
Fisher's exact test for biased strand distribution p-value = 7.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTT  >  NC_000913/2470294‑2470589
                                                                                                                                                   |                                                                                                                                                     
tGGAATATCTCAATTTTTTGTGTCTTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTTTTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACaa                                                                                                                                                     <  1:1499210/150‑1 (MQ=255)
   aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAATCAAAACAAAAg                                                                                                                                                  >  1:2335999/1‑150 (MQ=255)
   aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTCAAAGAAAAAAATACAAAAg                                                                                                                                                  >  1:3639820/1‑150 (MQ=255)
   aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATAAAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAGACAATat                                                                                                                                                  >  1:2325510/1‑149 (MQ=255)
        ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAAAAAAAAAAAAAAAAAGCGAtt                                                                                                                                            >  1:1293661/1‑149 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAAAGATGt                                                                                                                                           >  1:2915238/1‑150 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATAAATACAAAAACGATGt                                                                                                                                           >  1:3403465/1‑150 (MQ=255)
          cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                                          <  1:1577202/150‑1 (MQ=255)
             ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                                       <  2:1124448/150‑1 (MQ=255)
                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAAATCAAAGAAATAAAGACAAAAGCGATGTAAAGAACTAAGtt                                                                                                                              >  2:871635/1‑150 (MQ=255)
                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATATTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                              <  1:1725296/150‑1 (MQ=255)
                        aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAATAAATAAAGACAATAGCGATGTAAAGCACGAAGTTTa                                                                                                                            >  2:800358/1‑150 (MQ=255)
                         tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAATAAAATATCGATTTAAATCAATAAGTTTaa                                                                                                                           >  2:2603909/1‑150 (MQ=255)
                            ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                                                                                                        <  1:3848553/150‑1 (MQ=255)
                             tCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATTAAATAAATAAAAAAAATAAAAAAAAAAGAGATGTAAAGAACGAAGTTTaaaatt                                                                                                                       >  1:1803043/1‑147 (MQ=255)
                              cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                                                                                                                      <  1:2682620/150‑1 (MQ=255)
                              cGGGGAAGGATATAGAGATCTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                                                                                                                      <  1:2683322/150‑1 (MQ=255)
                               ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                                     <  2:320977/150‑1 (MQ=255)
                               ggggAAGGATATAGAGAGTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                                     <  1:1854493/150‑1 (MQ=255)
                                gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTGAAAGAAATCAAGAAAAAAAAGATGGAAAGAAAAAATTTTAAAtttggc                                                                                                                    >  1:2565451/1‑146 (MQ=255)
                                  gAATGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                                  <  2:3986384/150‑1 (MQ=255)
                                  gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTTAAATAATTAAAAGAAATCAAGAGAATAGCGGGGTGATGCACAAAGTTTAAATTgcgcgc                                                                                                                  >  1:2656867/1‑150 (MQ=255)
                                  gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                                  <  2:152005/150‑1 (MQ=255)
                                    aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                                <  1:3822337/150‑1 (MQ=255)
                                     ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                               <  1:2037907/150‑1 (MQ=255)
                                           agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAAGAAAAAAGCGATTTAAATCAAAAAGTTTAAATTTGGCGCCCATCAAgg                                                                                                         >  1:482910/1‑150 (MQ=255)
                                            gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                        <  1:3225786/150‑1 (MQ=255)
                                             agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAAACAATACAATAGCGCTGTTAAGAACTAAGGTTAAATTTAGCGCCCATCCTGgcg                                                                                                       >  1:1691536/1‑150 (MQ=255)
                                              gaTTTAAGCACGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                                      <  1:3276111/150‑1 (MQ=255)
                                               aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAACTCAAGAAAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATGgcgcg                                                                                                     >  2:334832/1‑150 (MQ=255)
                                                  tAAGCAATCATTTATTTTGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                                  <  2:2133585/150‑2 (MQ=255)
                                                    agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAAATAAAAAAAAAAAAAAAAAAATCGATTTAAAACAAAAAAGTTAAATTTAGCGCCAAACAATGCGCGCACaa                                                                                                >  1:3265314/1‑150 (MQ=255)
                                                    agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                                <  2:2764574/150‑1 (MQ=255)
                                                      caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                              <  2:579379/150‑1 (MQ=255)
                                                      caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATTAAATAAATCAAAAAAAAAAAGAAAATAGCGTTGTAAAGCAATAAGTGTAAATTTAGCGCCCCGCAGGGGGCGCCCAAGc                                                                                              >  1:527198/1‑150 (MQ=255)
                                                        agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAAAAAAAAAAAAAGCGATTGTAAACAAAAATATTTAAATTAGCGCCAAACAATGCGCGCACaaacaa                                                                                            >  1:3715668/1‑146 (MQ=255)
                                                           atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAAAAAATAAAATAAAAAAAAAAAAAAGAGCTGTAAAAAACGAAGATTAAATTTAGCGCCACACAGGGCGCGCACCAGcaaaaa                                                                                         >  1:2892655/1‑145 (MQ=255)
                                                            tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAAACAACACAAAACCGAAGTAAAGCACTAAGTTTAAATTTCGCGCCCAACATGGCGCCCCCAACCAataaa                                                                                        >  2:1963353/1‑149 (MQ=255)
                                                                    gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGGGCCAATCATGGCGCGCACAAGCTATAATAACAACCt                                                                                >  2:2410897/1‑150 (MQ=255)
                                                                            ctcccgctcgACCTCTGCTCACCCATAACCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGACGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                        <  1:2195978/143‑1 (MQ=255)
                                                                             tttCGTTCGACCTTTTCTTATAAATAACAGTCTTTTTTTTAATTTATAAATTAATTCAAATAAATCAAGACAATAGCTTTGTAAAGCACTAATTTTAAATTGCGCGCCAATCAGGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                       <  1:1534448/150‑1 (MQ=255)
                                                                             tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                       <  2:66524/150‑1 (MQ=255)
                                                                             tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                       <  2:3446026/150‑1 (MQ=255)
                                                                                 gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATGGCGCGCCCAAGCTATAATACCAACCTAATTTctcctcct                                                                   >  2:626418/1‑150 (MQ=255)
                                                                                  ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAATAGCGAAGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATTGCGCGCACAAGCTAAAAAACCAACCTAATTTctccaccac                                                                  >  2:2728018/1‑148 (MQ=255)
                                                                                     gACGTTTGATTATACATAAGATTTTTTTTTTTAATTTATAAAATAATTAAATTAAATCTAGACAATAGCGATGTATAGCACTATTTTTAAATTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                               <  1:3380641/150‑1 (MQ=255)
                                                                                     gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                               <  2:1293661/150‑1 (MQ=255)
                                                                                      aCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTag                                                              <  2:1197745/150‑1 (MQ=255)
                                                                                       ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                             <  2:3418220/150‑1 (MQ=255)
                                                                                       ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                             >  2:338699/1‑150 (MQ=255)
                                                                                       ccTTTGCTTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                             <  1:2293194/150‑1 (MQ=255)
                                                                                             cTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                       <  2:2974101/150‑1 (MQ=255)
                                                                                             cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACAAAGTTTAAATTTCGCGCCAATCATTGCGCGCACAAGCAATAATACCAAACTAATTTATACTACTCATAGATTGACt                                                       >  2:3566682/1‑150 (MQ=255)
                                                                                              ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAGACAAAATCGAATTAAAGCACAAAGTTTAAATTGCGCGCCAATAATTGCGCGCACAATCTATAATAACAACCAAATTTATCCTCCTCGTAGAGTGAcaa                                                      >  2:417056/1‑148 (MQ=255)
                                                                                               tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                                     <  2:3625531/150‑1 (MQ=255)
                                                                                                        cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                            <  1:4133542/150‑1 (MQ=255)
                                                                                                        cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                            <  2:2338967/150‑1 (MQ=255)
                                                                                                              tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                      >  1:1738428/1‑150 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAACTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAGACCAACCTACTTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                  <  2:3210419/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                  <  2:3995906/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                  >  1:377910/1‑150 (MQ=255)
                                                                                                                     aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                               >  1:1547428/1‑150 (MQ=255)
                                                                                                                      aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                              >  2:2733518/1‑150 (MQ=255)
                                                                                                                       tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACCAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATCGAATTGCGGTAtt                             >  1:2240838/1‑150 (MQ=255)
                                                                                                                        ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                            >  2:3001511/1‑150 (MQ=255)
                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                           <  1:957786/150‑1 (MQ=255)
                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                          <  2:2325510/150‑1 (MQ=255)
                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                          <  2:2325106/150‑1 (MQ=255)
                                                                                                                           taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                         >  1:3921231/1‑150 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                        <  1:3989092/150‑1 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                        >  2:1874173/1‑150 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATGGACTa                        >  2:2063687/1‑150 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAGAGAATTGCGGTATTGACTa                        >  2:3681285/1‑150 (MQ=255)
                                                                                                                              aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                      >  1:884203/1‑150 (MQ=255)
                                                                                                                               ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTCAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa                     <  2:3912861/150‑1 (MQ=255)
                                                                                                                                taatTCAAAGAAATAAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGAGCCAATCATGGCGCCCACAAGCTATAATACCAACCAAATTAAACATCCTAATAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa                    >  2:137742/1‑150 (MQ=255)
                                                                                                                                 aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                   <  1:3653811/150‑1 (MQ=255)
                                                                                                                                 aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                   <  2:3456126/150‑1 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                 >  1:1274184/1‑150 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                 >  1:1413740/1‑150 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACCAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                 <  2:149811/150‑1 (MQ=255)
                                                                                                                                    tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAGACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                <  2:2886151/150‑1 (MQ=255)
                                                                                                                                     cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc               <  1:2410897/150‑1 (MQ=255)
                                                                                                                                     cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc               >  1:1421199/1‑150 (MQ=255)
                                                                                                                                      aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc              <  2:4058258/150‑1 (MQ=255)
                                                                                                                                      aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc              >  2:719345/1‑150 (MQ=255)
                                                                                                                                      aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc              >  2:3790924/1‑150 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt             >  2:3255218/1‑150 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt             >  2:2168924/1‑150 (MQ=255)
                                                                                                                                         gAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa           <  2:2489591/150‑1 (MQ=255)
                                                                                                                                             tCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg       >  1:160405/1‑150 (MQ=255)
                                                                                                                                             tCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg       >  2:847652/1‑150 (MQ=255)
                                                                                                                                              cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt      <  1:1963353/150‑1 (MQ=255)
                                                                                                                                               aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     <  2:1747859/150‑1 (MQ=255)
                                                                                                                                               aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     >  1:4063897/1‑150 (MQ=255)
                                                                                                                                               aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     >  1:1094572/1‑150 (MQ=255)
                                                                                                                                               aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     >  2:2382862/1‑150 (MQ=255)
                                                                                                                                               aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCACTGATTCgtt     >  2:3542307/1‑150 (MQ=255)
                                                                                                                                                 gACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgt   >  2:2720531/1‑150 (MQ=255)
                                                                                                                                                  aCAATAGAGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCAGAAGCTATAATACCAACCTTATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgtt  <  2:522658/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTT  >  NC_000913/2470294‑2470589

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: