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→GAA) gtrSserotype‑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.

TGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTT  >  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)
        ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAAAAAAAAAAAAAAAAAGCGA‑tt                                                                                                                                            >  1:1293661/1‑149 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAAAGA‑TGt                                                                                                                                           >  1:2915238/1‑150 (MQ=255)
         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATAAATACAAAAACGA‑TGt                                                                                                                                           >  1:3403465/1‑150 (MQ=255)
          cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTa                                                                                                                                          <  1:1577202/150‑1 (MQ=255)
             ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAg                                                                                                                                       <  2:1124448/150‑1 (MQ=255)
                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAAATCAAAGAAATAAAGACAAAAGCGA‑TGTAAAGAACTAAGtt                                                                                                                              >  2:871635/1‑150 (MQ=255)
                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATATTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGtt                                                                                                                              <  1:1725296/150‑1 (MQ=255)
                        aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAATAAATAAAGACAATAGCGA‑TGTAAAGCACGAAGTTTa                                                                                                                            >  2:800358/1‑150 (MQ=255)
                         tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAATAAAATATCGA‑TTTAAATCAATAAGTTTaa                                                                                                                           >  2:2603909/1‑150 (MQ=255)
                            ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAAtt                                                                                                                        <  1:3848553/150‑1 (MQ=255)
                             tCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATTAAATAAATAAAAAAAATAAAAAAAAAAGAGA‑TGTAAAGAACGAAGTTTaaaatt                                                                                                                       >  1:1803043/1‑147 (MQ=255)
                              cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTgc                                                                                                                      <  1:2682620/150‑1 (MQ=255)
                              cGGGGAAGGATATAGAGATCTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTgc                                                                                                                      <  1:2683322/150‑1 (MQ=255)
                               ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTgcg                                                                                                                     <  2:320977/150‑1 (MQ=255)
                               ggggAAGGATATAGAGAGTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTgcg                                                                                                                     <  1:1854493/150‑1 (MQ=255)
                                gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTGAAAGAAATCAAGAAAAAAAAGA‑TGGAAAGAAAAAATTTTAAAtttggc                                                                                                                    >  1:2565451/1‑146 (MQ=255)
                                  gAATGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                                  <  2:3986384/150‑1 (MQ=255)
                                  gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTTAAATAATTAAAAGAAATCAAGAGAATAGCGG‑GGTGATGCACAAAGTTTAAATTgcgcgc                                                                                                                  >  1:2656867/1‑150 (MQ=255)
                                  gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                                  <  2:152005/150‑1 (MQ=255)
                                    aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                                <  1:3822337/150‑1 (MQ=255)
                                     ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                               <  1:2037907/150‑1 (MQ=255)
                                           agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAAGAAAAAAGCGA‑TTTAAATCAAAAAGTTTAAATTTGGCGCCCATCAAgg                                                                                                         >  1:482910/1‑150 (MQ=255)
                                            gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                        <  1:3225786/150‑1 (MQ=255)
                                             agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAAACAATACAATAGCGC‑TGTTAAGAACTAAGGTTAAATTTAGCGCCCATCCTGgcg                                                                                                       >  1:1691536/1‑150 (MQ=255)
                                              gaTTTAAGCACGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                                      <  1:3276111/150‑1 (MQ=255)
                                               aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAACTCAAGAAAATATCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATGgcgcg                                                                                                     >  2:334832/1‑150 (MQ=255)
                                                  tAAGCAATCATTTATTTTGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                                  <  2:2133585/150‑2 (MQ=255)
                                                    agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAAATAAAAAAAAAAAAAAAAAAATCGA‑TTTAAAACAAAAAAGTTAAATTTAGCGCCAAACAATGCGCGCACaa                                                                                                >  1:3265314/1‑150 (MQ=255)
                                                    agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                                <  2:2764574/150‑1 (MQ=255)
                                                      caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                              <  2:579379/150‑1 (MQ=255)
                                                      caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATTAAATAAATCAAAAAAAAAAAGAAAATAGCGT‑TGTAAAGCAATAAGTGTAAATTTAGCGCCCCGCAGGGGGCGCCCAAGc                                                                                              >  1:527198/1‑150 (MQ=255)
                                                        agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAAAAAAAAAAAAAGCGATTGTAAACAAAAATATTTAAATT‑AGCGCCAAACAATGCGCGCACaaacaa                                                                                            >  1:3715668/1‑146 (MQ=255)
                                                           atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAAAAAATAAAATAAAAAAAAAAAAAAGAGC‑TGTAAAAAACGAAGATTAAATTTAGCGCCACACAGGGCGCGCACCAGcaaaaa                                                                                         >  1:2892655/1‑145 (MQ=255)
                                                            tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAAACAACACAAAACCGA‑AGTAAAGCACTAAGTTTAAATTTCGCGCCCAACATGGCGCCCCCAACCAataaa                                                                                        >  2:1963353/1‑149 (MQ=255)
                                                                    gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTTCGGGCCAATCATGGCGCGCACAAGCTATAATAACAACCt                                                                                >  2:2410897/1‑150 (MQ=255)
                                                                            ctcccgctcgACCTCTGCTCACCCATAACCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑CGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                        <  1:2195978/143‑1 (MQ=255)
                                                                             tttCGTTCGACCTTTTCTTATAAATAACAGTCTTTTTTTTAATTTATAAATTAATTCAAATAAATCAAGACAATAGCTT‑TGTAAAGCACTAATTTTAAATTGCGCGCCAATCAGGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                       <  1:1534448/150‑1 (MQ=255)
                                                                             tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                       <  2:66524/150‑1 (MQ=255)
                                                                             tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                       <  2:3446026/150‑1 (MQ=255)
                                                                                 gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATATCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATGGCGCGCCCAAGCTATAATACCAACCTAATTTctcctcct                                                                   >  2:626418/1‑150 (MQ=255)
                                                                                  ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAATAGCGA‑AGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATTGCGCGCACAAGCTAAAAAACCAACCTAATTTctccaccac                                                                  >  2:2728018/1‑148 (MQ=255)
                                                                                     gACGTTTGATTATACATAAGATTTTTTTTTTTAATTTATAAAATAATTAAATTAAATCTAGACAATAGCGA‑TGTATAGCACTATTTTTAAATTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                               <  1:3380641/150‑1 (MQ=255)
                                                                                     gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                               <  2:1293661/150‑1 (MQ=255)
                                                                                      aCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTag                                                              <  2:1197745/150‑1 (MQ=255)
                                                                                       ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                             <  2:3418220/150‑1 (MQ=255)
                                                                                       ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                             >  2:338699/1‑150 (MQ=255)
                                                                                       ccTTTGCTTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                             <  1:2293194/150‑1 (MQ=255)
                                                                                             cTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                       <  2:2974101/150‑1 (MQ=255)
                                                                                             cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACAAAGTTTAAATTTCGCGCCAATCATTGCGCGCACAAGCAATAATACCAAACTAATTTATACTACTCATAGATTGACt                                                       >  2:3566682/1‑150 (MQ=255)
                                                                                              ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAGACAAAATCGA‑ATTAAAGCACAAAGTTTAAATTGCGCGCCAATAATTGCGCGCACAATCTATAATAACAACCAAATTTATCCTCCTCGTAGAGTGAcaa                                                      >  2:417056/1‑148 (MQ=255)
                                                                                               tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                                     <  2:3625531/150‑1 (MQ=255)
                                                                                                        cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                            <  1:4133542/150‑1 (MQ=255)
                                                                                                        cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                            <  2:2338967/150‑1 (MQ=255)
                                                                                                              tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                      >  1:1738428/1‑150 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAACTCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAGACCAACCTACTTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                  <  2:3210419/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                  <  2:3995906/150‑1 (MQ=255)
                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                  >  1:377910/1‑150 (MQ=255)
                                                                                                                     aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                               >  1:1547428/1‑150 (MQ=255)
                                                                                                                      aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                              >  2:2733518/1‑150 (MQ=255)
                                                                                                                       tttATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACCAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATCGAATTGCGGTAtt                             >  1:2240838/1‑150 (MQ=255)
                                                                                                                        ttATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                            >  2:3001511/1‑150 (MQ=255)
                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                           <  1:957786/150‑1 (MQ=255)
                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                          <  2:2325510/150‑1 (MQ=255)
                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                          <  2:2325106/150‑1 (MQ=255)
                                                                                                                           taAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                         >  1:3921231/1‑150 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                        <  1:3989092/150‑1 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                        >  2:1874173/1‑150 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATGGACTa                        >  2:2063687/1‑150 (MQ=255)
                                                                                                                            aaaaTAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAGAGAATTGCGGTATTGACTa                        >  2:3681285/1‑150 (MQ=255)
                                                                                                                              aataatTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                      >  1:884203/1‑150 (MQ=255)
                                                                                                                               ataatTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTCAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa                     <  2:3912861/150‑1 (MQ=255)
                                                                                                                                taatTCAAAGAAATAAAGACAAAAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGAGCCAATCATGGCGCCCACAAGCTATAATACCAACCAAATTAAACATCCTAATAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa                    >  2:137742/1‑150 (MQ=255)
                                                                                                                                 aatTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                   <  1:3653811/150‑1 (MQ=255)
                                                                                                                                 aatTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa                   <  2:3456126/150‑1 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                 >  1:1274184/1‑150 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                 >  1:1413740/1‑150 (MQ=255)
                                                                                                                                   ttCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACCAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                 <  2:149811/150‑1 (MQ=255)
                                                                                                                                    tCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAGACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                <  2:2886151/150‑1 (MQ=255)
                                                                                                                                     cAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc               <  1:2410897/150‑1 (MQ=255)
                                                                                                                                     cAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc               >  1:1421199/1‑150 (MQ=255)
                                                                                                                                      aaaGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc              <  2:4058258/150‑1 (MQ=255)
                                                                                                                                      aaaGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc              >  2:719345/1‑150 (MQ=255)
                                                                                                                                      aaaGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc              >  2:3790924/1‑150 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt             >  2:3255218/1‑150 (MQ=255)
                                                                                                                                       aaGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt             >  2:2168924/1‑150 (MQ=255)
                                                                                                                                         gAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa           <  2:2489591/150‑1 (MQ=255)
                                                                                                                                             tCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg       >  1:160405/1‑150 (MQ=255)
                                                                                                                                             tCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg       >  2:847652/1‑150 (MQ=255)
                                                                                                                                              cAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt      <  1:1963353/150‑1 (MQ=255)
                                                                                                                                               aaGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     <  2:1747859/150‑1 (MQ=255)
                                                                                                                                               aaGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     >  1:4063897/1‑150 (MQ=255)
                                                                                                                                               aaGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     >  1:1094572/1‑150 (MQ=255)
                                                                                                                                               aaGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt     >  2:2382862/1‑150 (MQ=255)
                                                                                                                                               aaGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCACTGATTCgtt     >  2:3542307/1‑150 (MQ=255)
                                                                                                                                                 gACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgt   >  2:2720531/1‑150 (MQ=255)
                                                                                                                                                  aCAATAGAGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCAGAAGCTATAATACCAACCTTATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttgtt  <  2:522658/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGA‑TGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTT  >  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: ‑