Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,482 C→A 9.4% intergenic (+20/+595) yfdI → / ← yfdK CPS‑53 (KpLE1) prophage; serotype specific glucosyl transferase/CPS‑53 (KpLE1) prophage; putative tail fiber assembly protein YfdK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4820CA 9.4% 462.1 / 24.1 142intergenic (+20/+595)yfdI/yfdKCPS‑53 (KpLE1) prophage; serotype specific glucosyl transferase/CPS‑53 (KpLE1) prophage; putative tail fiber assembly protein YfdK
Reads supporting (aligned to +/- strand):  ref base C (53/75);  new base A (6/8);  total (59/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 142 total aligned reads displayed.

TAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATG  >  NC_000913/2470336‑2470630
                                                                                                                                                  |                                                                                                                                                    
taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                                                                   >  1:1030086/1‑150 (MQ=255)
 agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                                                                                  <  2:1068148/150‑1 (MQ=255)
  gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                                                                 >  2:451390/1‑150 (MQ=255)
     aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg                                                                                                                                              <  1:435518/150‑1 (MQ=255)
         aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaat                                                                                                                                          <  2:478566/150‑3 (MQ=255)
          agcaagcaTTTATTTGGAATGTATTTTTCGTTCGTCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                                                                         <  2:1170684/150‑1 (MQ=255)
          agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca                                                                                                                                            <  2:1073920/147‑1 (MQ=255)
          agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca                                                                                                                                            >  1:1073920/1‑147 (MQ=255)
           gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaattt                                                                                                                                        <  1:1014000/150‑5 (MQ=255)
           gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg                                                                                                                                        >  1:1230561/1‑150 (MQ=255)
           gcaagcaTTTATTTGGAATGTATTTTTCGTTCCACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaattt                                                                                                                                        <  1:154339/150‑5 (MQ=255)
            caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                                                                       <  2:1194966/150‑1 (MQ=255)
            caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAATGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaattta                                                                                                                                       <  2:649540/150‑6 (MQ=255)
            caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                                                                       <  1:1159903/150‑1 (MQ=255)
             aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTAGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatttaa                                                                                                                                      <  1:890295/150‑7 (MQ=255)
                 atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                                                                                  >  1:610633/1‑150 (MQ=255)
                  tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                                                                                                                 <  1:1285863/150‑1 (MQ=255)
                  tttatttGGAATGTATTTTTCGTTCCACCTTTGCTTTAACATATCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                                                                                                                 <  1:1079997/150‑1 (MQ=255)
                    tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                                                                                                               >  2:830133/1‑150 (MQ=255)
                    tatttGGAATGTATTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                                                                              >  2:391829/1‑150 (MQ=255)
                     atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                                                                              <  1:39136/150‑1 (MQ=255)
                     atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                                                                              <  2:78650/150‑1 (MQ=255)
                       ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCaa                                                                                                                            >  2:118190/1‑150 (MQ=255)
                         ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                                                                          >  1:899555/1‑150 (MQ=255)
                         ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                                                                          <  2:70902/150‑1 (MQ=255)
                           aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTa                                                                                                                        >  2:90379/1‑150 (MQ=255)
                           aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTa                                                                                                                        >  2:325310/1‑150 (MQ=255)
                               tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                                                                                    >  2:894029/1‑150 (MQ=255)
                                aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                                                                                   <  2:415877/150‑1 (MQ=255)
                                 tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtct                                                                                                                  >  2:1277514/1‑150 (MQ=255)
                                 tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtct                                                                                                                  <  2:1162591/150‑1 (MQ=255)
                                   tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGGGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                                                                >  1:177083/1‑150 (MQ=255)
                                    ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcct                                                                                                               <  2:123034/150‑1 (MQ=255)
                                       gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                                                                            <  1:608868/150‑1 (MQ=255)
                                        ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcctc                                                                                                           >  2:185982/1‑150 (MQ=255)
                                           gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGAGCGCCAATCATGGCGCGCACAAGCAATAATACCAACCTAATTTCTCATCATCTTa                                                                                                        >  2:19282/1‑150 (MQ=255)
                                           gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                                                                        <  1:1102112/150‑1 (MQ=255)
                                           caCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                                                                        <  1:321479/149‑1 (MQ=255)
                                            aCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCTCGCCAATCATGGCGAGCACAAGCTATAATACCAACCTAATTTCTCATCCTCTTag                                                                                                       <  2:1048240/150‑1 (MQ=255)
                                            aCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTag                                                                                                       <  2:152246/150‑1 (MQ=255)
                                             ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                                                                      >  2:1121415/1‑150 (MQ=255)
                                             ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                                                                      <  2:462966/150‑1 (MQ=255)
                                                ttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg                                                                                                   <  2:24984/150‑1 (MQ=255)
                                                ttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg                                                                                                   <  2:1248346/150‑1 (MQ=255)
                                                ttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg                                                                                                   <  1:791584/150‑1 (MQ=255)
                                                  gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                                                                 >  2:1338849/1‑150 (MQ=255)
                                                   cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                                                                >  1:396135/1‑150 (MQ=255)
                                                   cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTACTCCTCCTCTTAGAGTGACt                                                                                                >  1:1393250/1‑150 (MQ=255)
                                                         cATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct                                                                                          >  1:628720/1‑150 (MQ=255)
                                                            aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                                                                       <  1:751515/150‑1 (MQ=255)
                                                              cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                                                                     <  2:906926/150‑1 (MQ=255)
                                                                gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                                                                   >  2:840255/1‑150 (MQ=255)
                                                                gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                                                                   >  1:1343712/1‑150 (MQ=255)
                                                                 tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                                                                                  <  1:1338849/150‑1 (MQ=255)
                                                                   gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                                                                                >  2:384853/1‑150 (MQ=255)
                                                                       ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                                                            <  1:647189/150‑1 (MQ=255)
                                                                       ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                                                                            <  2:1030086/150‑1 (MQ=255)
                                                                        tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                                                           <  1:762259/150‑1 (MQ=255)
                                                                         ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                                                                          <  1:391941/150‑1 (MQ=255)
                                                                           aaTTTTTAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                                                                        >  2:873965/1‑150 (MQ=255)
                                                                           aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                                                                        >  2:993908/1‑150 (MQ=255)
                                                                           aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                                                                        <  1:260094/150‑1 (MQ=255)
                                                                             tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                                                                      >  1:1056312/1‑150 (MQ=255)
                                                                                   aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCATGATAGAATTGCGGTATTGACTAt                                                                >  1:1053223/1‑150 (MQ=255)
                                                                                         ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                                                          >  1:911106/1‑150 (MQ=255)
                                                                                          tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                                                         <  2:52578/150‑1 (MQ=255)
                                                                                            aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCTCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc                                                       >  1:190521/1‑150 (MQ=255)
                                                                                              aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg                                                     <  1:185982/150‑1 (MQ=255)
                                                                                                aaaTCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt                                                   <  1:391829/150‑1 (MQ=255)
                                                                                                  aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                                                                                    >  1:86526/1‑115 (MQ=255)
                                                                                                  aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                                                                                    <  2:86526/115‑1 (MQ=255)
                                                                                                  aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc                                                 <  1:51793/150‑1 (MQ=255)
                                                                                                   tCAAGACAATAGCGATGTAAATCACTAAGTTTAAATTGCTCGCCAATCATGGCGCTCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg                                                <  1:504050/150‑1 (MQ=255)
                                                                                                    cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt                                               >  1:719196/1‑150 (MQ=255)
                                                                                                         cAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGttt                                          <  1:1217532/150‑1 (MQ=255)
                                                                                                          aaTAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTa                                         <  1:1232894/150‑1 (MQ=255)
                                                                                                           aTAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAATCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTAt                                        <  1:2494/150‑1 (MQ=255)
                                                                                                           aTAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTAt                                        >  2:475452/1‑150 (MQ=255)
                                                                                                            tAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTAtt                                       <  1:384853/150‑1 (MQ=255)
                                                                                                            tAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTAtt                                       >  1:1200890/1‑150 (MQ=255)
                                                                                                              gCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGt                                     <  2:632375/150‑1 (MQ=255)
                                                                                                                gATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTaa                                   <  2:610633/150‑1 (MQ=255)
                                                                                                                 aTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAat                                  <  1:580937/150‑1 (MQ=255)
                                                                                                                      aaGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATAtctc                             <  2:507640/150‑1 (MQ=255)
                                                                                                                       aGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcc                            >  1:444841/1‑150 (MQ=255)
                                                                                                                          aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct                         <  1:1016247/150‑1 (MQ=255)
                                                                                                                          aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct                         <  1:470441/150‑1 (MQ=255)
                                                                                                                           cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTa                        <  1:1305920/150‑1 (MQ=255)
                                                                                                                           cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTa                        <  1:118190/150‑1 (MQ=255)
                                                                                                                           cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTa                        <  1:840255/150‑1 (MQ=255)
                                                                                                                              aGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCt                     <  2:1230561/150‑1 (MQ=255)
                                                                                                                               gTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTg                    >  1:1137007/1‑150 (MQ=255)
                                                                                                                                     tgtGCGCGCCATGATTGGCGCGCGCAGGCTATAATGCCAGCCTAATTTCTCCTCCTCTTAGCGTGGCTATAGCTCCGGATAGAGTTGGGGTATTGACGATCAAGTGCGCTGCTTCGTTGTGTATTGTATGATCTCCTGTAGCTGCGGagg              >  2:1219396/3‑148 (MQ=255)
                                                                                                                                     tgtGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACg              >  2:276588/3‑150 (MQ=255)
                                                                                                                                     aTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACg              >  2:8275/1‑150 (MQ=255)
                                                                                                                                      gtGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGa             >  1:377510/2‑150 (MQ=255)
                                                                                                                                       tGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                                                                  >  1:37356/1‑96 (MQ=255)
                                                                                                                                       tGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                                                                  <  2:37356/96‑1 (MQ=255)
                                                                                                                                           cgcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTACTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACt        <  2:1393250/150‑1 (MQ=255)
                                                                                                                                             ccAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTa      >  2:956263/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATG  >  NC_000913/2470336‑2470630

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